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Abstract 

We give a self-contained introduction to the physics of ultracold atoms using functional integral techniques. Based 
on a consideration of the relevant length scales, we derive the universal effective low energy Hamiltonian describing 
ultracold alkali atoms. We then introduce the concept of the effective action, which generalizes the classical action 
principle to full quantum status and provides an intuitive and versatile tool for practical calculations. This framework 
is applied to weakly interacting degenerate bosons and fermions in the spatial continuum. In particular, we discuss 
the related BEC and BCS quantum condensation mechanisms. We then turn to the BCS-BEC crossover, which 
interpolates between both phenomena, and which is realized experimentally in the vicinity of a Feshbach resonance. 
For its description, we introduce the Functional Renormalization Group approach. After a general discussion of the 
method in the cold atoms context, we present a detailed and pedagogical application to the crossover problem. This 
not only provides the physical mechanism underlying this phenomenon. More generally, it also reveals how the 
renormalization group can be used as a tool to capture physics at all scales, from few-body scattering on microscopic 
scales, through the finite temperature phase diagram governed by many-body length scales, up to critical phenomena 
dictating long distance physics at the phase transition. 

The presentation aims to equip students at the beginning PhD level with knowledge on key physical phenomena and 
flexible tools for their description, and should enable to embark upon practical calculations in this field. 
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1. Introduction 

Cold atomic many-body systems make up a young 
and rapidly evolving area of modern physics. The field 
was born in 1995, where an almost pure and weakly in- 
teracting Bose-Einstein condensate (BEC) was created 
in alkaline atomic vapors [Q~|[2). This state of matter ex- 
hibits macroscopic phase coherence as a manifest fea- 
ture of many -body quantum mechanics. Soon after, fur- 
ther milestones were achieved - among them, the obser- 
vation of the quantum phase transition from a Mott in- 
sulator to a superfluid in optical lattices, to which atoms 
are confined due to their interaction with light (3]|4), and 
the implementation of the BCS-BEC crossover realizing 
strongly interacting fermion ensembles ||5] |6) . 

These achievements reflect both the high degree of 
control in the manipulation of the atomic constituents, 
as well as the remarkable tunability of the scales and in- 
teractions governing the system. From a practical point 
of view, this prepares the ground for efficient quan- 
tum simulations using cold atomic gases. A promi- 
nent example is provided by the determination of the 
ground state of the two-dimensional Fermi-Hubbard 
model, which remains a challenge to theory [7], but can 
be implemented in cold atomic gases with high accu- 
racy as a direct simulator to gain insight into its low 
temperature quantum physics, once sufficient cooling is 
achieved J8). With increasing control over the micro- 
scopic constituents in optical lattices, a long term goal is 
to go beyond such special task quantum simulators, and 
to build truly programmable, universal quantum simu- 
lation devices ll9T H2l . From a theoretical point of view, 
cold atomic samples offer a unique testbed for modem 
nonperturbative quantum field theoretical approaches: 
While a precise microscopic understanding is often not 
available in complex many-body quantum systems re- 
alized in condensed matter or high energy physics, the 
ability to experimentally probe a cold atomic sample at 



all scales, and to manipulate its microphysics in a con- 
trolled way, allows for a direct comparison of experi- 
ment and theory in a strongly interacting context. The 
transition from microscopic simplicity to macroscopic 
complexity, performed in a given theoretical approach, 
can be benchmarked in a direct way. 

Beyond these general considerations, cold atomic 
gases also offer a number of physical situations which 
do not have an immediate counterpart in other branches 
of physics. Without aiming at completeness, but for the 
sake of giving a flavor, we briefly point out two direc- 
tions which recently attracted interest. The first one 
harnesses quantum optical manipulation tools for con- 
structing microscopic Hamiltonians with exotic interac- 
tions. For instance, long range interactions such as l/r 3 , 
in part with strong spatial anisotropy, are obtained in the 
context of dipolar atoms and polar molecules lfl3UT5l or 
Rydberg dressed atoms fl6l . Moreover, it is possible to 
realize multicomponent interactions with high degree of 
symmetry such as SU(iV). This is achieved by control- 
ling and addressing internal states of atoms, which pos- 
sess a rich level structures, like earth alkaline-like atoms 
lfT7lfT9l . Reliably extracting the many-body physics 
of each of these systems poses its own specific chal- 
lenges to theory. A second direction is given by non- 
equilibrium physics with cold atoms. This comprises, 
on the one hand, the dynamics of closed systems, where 
key questions are related to the relaxation dynamics to- 
wards thermodynamic equilibrium Il20l424l . the behav- 
ior of a many -body system following a quench of micro- 
scopic parameters [25-28], or dynamical phenomena 
such as the propagation of local perturbations through 
the system [29 1 . On the other hand, this concerns open 
systems, where many-body ensembles are driven far 
from thermodynamic equilibrium by coupling to effec- 
tive external reservoirs - which may occur naturally, or 
via specific reservoir engineering. Such systems may 
exhibit stable non- equilibrium stationary states with in- 
triguing quantum mechanical properties and rich phase 
diagrams [30-33 1. In the non-equilibrium context, it is 
a first challenge for both experiment and theory to iden- 
tify situations which exhibit a sufficient degree of uni- 
versality, i.e. with a phenomenology which occurs in 
classes of systems and settings beyond a particular real- 
ization. Furthermore, this calls for the development of 
flexible theoretical tools - both numerical and analytical 
- to describe them efficiently. 

The aim of these lecture notes is to act as a door- 
opener to this exciting field. It should equip students 
at the beginning of their theoretical PhD studies with 
the knowledge of key physical many-body phenomena 
in ultracold quantum gases in the spatial continuum, as 
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well as with flexible functional techniques for their de- 
scription. In particular, we introduce the concept of the 
effective action, which generalizes the classical action 
principle to full status and is an intuitive and versatile 
tool for practical calculations. We describe the physics 
of weakly interacting degenerate bosons and fermions in 
this framework, and discuss the related BEC and BCS 
quantum condensation mechanisms. We also describe 
a Functional Renormalization Group (FRG) approach 
to ultracold atoms, for reviews see [7 34-37], and ex- 
tract the full finite temperature phase diagram for the 
BCS-BEC crossover within a simple approximation for 
the running effective action. This discussion not only 
provides the physical mechanism underlying this phe- 
nomenon, which interpolates between the above two 
cornerstones of quantum condensation phenomena with 
a particularly challenging strongly interacting regime 
in between. It also illustrates how the concept of the 
renormalization group can be used in practice beyond 
the realm of critical phenomena at very long distances. 
In fact, it provides a powerful tool to smoothly per- 
form the transition from micro- to macrophysics. This 
comprises the quantitative description of few-body scat- 
tering in the physical vacuum at short distances, the 
phase diagram and thermodynamics governed by many- 
body scales such as interparticle spacing and de Broglie 
wavelength, and finally also includes the description 
of critical behavior at large wavelength, in one unified 
framework. 

The lecture notes are organized as follows. In Sec. [2] 
we discuss the basic microscopic and many-body scales 
present in ultracold gases, and construct an effective 
Hamiltonian which universally describes the physics of 
bosonic and fermionic alkali gases at low energies, as 
appropriate for ultracold experiments. For the sake of 
a self-contained presentation, we also provide brief re- 
minders on thermodynamics and the quantum statistical 
mechanics of noninteracting bosons and fermions. In 
Sec. [3] we introduce the functional integral represen- 
tation of the quantum partition function (with technical 
details on the derivation provided in the appendix), and 
switch to a more intuitive object encoding the same in- 
formation, the effective action. We describe the key phe- 
nomenon of spontaneous symmetry breaking, and how 
to extract thermodynamic information from this object. 
We apply these concepts to weakly interacting bosons 
and fermions, with an emphasis on the condensation 
phenomenon and the nature of the low energy excita- 
tions in both cases. In Sec. |4] we introduce Wetterich's 
FRG framework for the effective action, which encodes 
the full information on the many-body problem in terms 
of an exact functional differential equation, and is ide- 



ally suited for the implementation of practical approx- 
imation schemes (truncations) beyond mean field the- 
ory. After a brief general discussion on the application 
of this framework to ultracold atoms systems, we apply 
it to the BCS-BEC crossover problem. The emphasis is 
on a detailed presentation of a simple truncation, which 
is able to produce the full finite temperature phase di- 
agram and already at this level demonstrates a number 
of improvements compared to extended mean field the- 
ories, due to a consistent inclusion of bosonic degrees 
of freedom. Since, on a technical level, this covers both 
the treatment of interacting fermionic and bosonic the- 
ories, it will enable the reader to embark upon practical 
calculations in this field. We provide a (subjective) list 
of challenges for the future. 

This work is based on lectures delivered by S. D. at 
the 49th Schladming Graduate School for Theoretical 
Physics. 

2. Basics of ultracold atomic physics 

2.1. Scales and interactions 

The physics of ultracold quantum gases is governed by 
the interplay of several scales. Tuning their relative size, 
it is possible to access different regions of the phase 
diagram of a given system, in this way exploring its 
physics. In this section, we show which scales are rele- 
vant in the context of alkali atoms. In particular, we will 
discuss the conditions under which we have an ultracold 
quantum gas. These model-independent considerations 
will also reveal why it is possible to formulate a sim- 
ple effective Hamiltonian, described by a few experi- 
mentally measurable parameters only, which governs all 
alkali (single valence electron) atoms. 

Given a homogeneous gas of atoms with density n in 
d spatial dimensions, we may write 

n = r d , (1) 

with t being the interparticle spacing. Indeed, consider 
a homogeneous system in a box of volume V. We divide 
this volume into cells of size l d each. Putting exactly 
one atom into each cell, it is possible to distribute N = 
V/{ d particles. Thus, we arrive at the density n = N/V = 

r d . 

Experiments on cold atoms are performed in either 
magnetic or optical traps (see [1, 38 , 39 1 and references 
therein). Therefore, the ground state of the many-body 
system will not be homogeneous. In particular, the den- 
sity depends on space. However, there are many cases 
where the picture of a locally homogeneous system is 
still valid and useful [40 1. In order to understand this, 
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we consider a time-independent external trapping po- 
tential of harmonic shape. We have 



2 ..2 



V ext (x) = 2 aJ o r 



(2) 



with r = and m being the mass of the atoms. The 
potential is characterized by the trapping frequency oiq. 
Equivalently, we may write V ext (x) = ^r(V/^ osc ) 2 with 
the oscillator length 



f 

1 (ISC 



h \ 



1/2 



(3) 



Thus, hcjQ and £ osc are the characteristic energy and 
length scales of the trap, respectively. We will later 
see that in a typical situation, f osc constitutes the by far 
largest length scale in the system. (It will, however, act 
as an infrared (IR) cutoff for very long wavelength fluc- 
tuations present e.g. at a critical point.) Accordingly, 
ha>o usually provides the smallest energy scale of the 
problem. 

If the physics under consideration takes place on 
much shorter distances than { osc , we can use this sep- 
aration of scales to work in a so-called local density 
approximation: Consider the density at points x t and 
x*2, respectively. We can then expand n{x\) — n{x2){\ + 
0(\x\ — X2\li 0S c))- Obviously, for both points being close 
to each other we can neglect the correction and assume 
the density to be locally constant. In particular, for 
large values of £ osc this may hold for subvolumes of the 
trapped cloud which contain many particles. The rules 
of thermodynamic equilibrium can then be applied to 
these small, homogeneous subvolumes. We will come 
back to this point in the section on thermodynamics of 
cold quantum gases. 

The statistical behavior of our trapped cloud is de- 
termined by the ratio between the interparticle spacing 
and the so-called thermal or de Broglie wavelength. To 
get an intuition for the latter quantity, consider a gas 
of atoms coupled to a heat bath of temperature T. The 
nonvanishing temperature induces a nonzero average ki- 
netic energy (p 2 )T/2m per spatial direction of the parti- 
cles. The de Broglie wavelength is the length scale asso- 

9 1 II 

ciated to this energy, according to At = h/(p ) T . More 
precisely, using p — hk — h/A (k and A the wavenum- 
ber and -length, respectively), we define A T as the de 
Broglie wavelength of a particle with kinetic energy 
p 2 /2m = nkftT. (The factor of ji is purely conventional 
but standard.) This leads to 



At — 



2nh 
mkftT 



2 \l/2 





Figure 1: Quantum degeneracy is reached when the thermal wave- 
length At is of the same order as the interparticle spacing C = n < d . 
In this regime, it is important for the statistics whether the particles are 
identical or not, leading to quantum many-body phenomena such as 
Bose condensation or a Fermi surface. We also indicate a typical order 
of magnitude for the scattering length a, which corresponds roughly 
to the radius of equivalent hard-core particles with contact interactions 
and cross section cr oc a 2 . 



Note that At ~ T~ 1 ^ 2 becomes large for decreasing tem- 
perature. The quantities i = l(n) and At = At(T) con- 
stitute the many-body length scales of the system due to 
nonzero density and temperature, respectively. 

Now we compare the length scales set by the interpar- 
ticle spacing and the de Broglie wavelength. Thinking 
of particles as being represented by wavepackets rather 
than pointlike objects, At determines the spread of these 
lumps. The ratio £/Ar is large if the wavepackets of 
the individual particles are widely separated and do not 
overlap. In this case the quantum nature of the parti- 
cles does not play a role. Indeed, we may follow the 
trajectory of an individual particle by subsequent im- 
ages, because position and momentum are determined 
simultaneously, i.e. the gas can be described classi- 
cally. However, for £/A T < 1, we are dealing with 
wavepackets which strongly overlap. The gas is then 
called quantum degenerate, or ultracold. Clearly, it is 
then no longer possible to distinguish the single atoms 
and their trajectories. In this case we rather have to deal 
with the whole many-body quantum system. The be- 
havior is then determined by quantum mechanics, with 
statistics resulting from the spin of the constituents; ul- 
tracold atoms allow for exploring both degenerate Bose 
and Fermi gases. 

The transition from the classical to the quantum de- 
generate regime occurs for nA d T 1, i.e. 



£/A T 1. 



(5) 



We visualize this situation in Fig. [T] The combination 

0) = n4 = (A T /() d (6) 

is called the phase space density. It indicates the num- 
ber of particles contained in a cube with linear extension 
set by the de Broglie wavelength. 



(4) 
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Interactions and effective Hamiltonian 

So far, our considerations did not depend on the inter- 
actions of the particles. The alkali atoms used in ul- 
tracold gas experiments are neutral and interact electro- 
magnetically through van der Waals forces. A typical 
interaction potential U(r) of two atoms separated by a 
distance r has a strongly repulsive part for small r. The 
physical origin of the latter is Pauli's principle which 
forbids the electron clouds of the two atoms to overlap. 
This repulsive part can typically be modeled by a term 
U(r) ~ 1/r 12 , but a hard core repulsion with infinite 
strength works as well. For larger distances, two atoms 
experience an attraction due to mutual polarization of 
the electron clouds. Each atom then acts as a small in- 
duced dipole, and they attract each other according to 
a van der Waals interaction U(r) ~ -1/r 6 . (We show 
the generic shape of the total interatomic potential, the 
Lennard-Jones potential, in Fig. [2]) We thus approxi- 
mate the microscopic interaction potential to be 





model potential with 




same scattering length 








>^ true interatomic 




/ potential U(r) 



Figure 2: The interatomic potential U(r) between two neutral atoms 
is of the Lennard-Jones type, with an attractive van der Waals tail 
~ 1/r 6 at large separations. From U(r) we can calculate the scattering 
length a, which is the only parameter relevant for low energy scatter- 
ing. The 5-like potential from Eq. j!4) , which is shown here in red, 
is an equally good description (and more handy for practical calcula- 
tions), as long as it has the same scattering length. The reason is that 
under ultracold conditions, the short distance details of U(r) are never 
resolved. 



°° (r < r ), 

-C 6 /r 6 (r > r ). 



(7) 



We can use this expression to provide a typical length 
scale, the van der Waals length, which characterizes the 
interactions. A typical length scale for zero (total) en- 
ergy scattering is obtained from equating kinetic energy 
of a particle with momentum p - H/€ v ^, and potential 
energy t/ v dw(A-dw), resulting in 



L 



vdW 



_ / mC 6 \ 

~\~w) 



1/4 



(8) 



For typical values of Ce, we find that ^ v dw = 
(50...200)ao (ciq - 5.3 x 10 2 nm the Bohr radius), which 
crucially is much smaller than the interparticle spacing 
and the thermal wavelength (cf. Tab. [TJ, 



£, Aj » ^vdw- 



(9) 



The many-body effects in an ultracold gas we are in- 
terested in thus never resolve physics beyond the van 
der Waals length. As a consequence, we will be able 
to specify an effective low energy Hamiltonian, valid on 
length scales > Ardw> as the microscopic starting point 
of our calculations. 

After indicating the rough scale associated to interac- 
tions, we now identify the relevant physical parameter 
which can be extracted from scattering experiments, the 
scattering length a [39|. This length scale character- 
izes two-body collisions and emerges universally as the 
sole parameter characterizing low energy collisions in 



potentials of sufficiently short range, such as 1/r 6 as we 
deal with here. To see this, let us consider low energy 
elastic scattering of two particles in a quantum mechan- 
ical framework. (As we explain below, we can assume 
only elastic two-body processes to be relevant; further 
note that our meaning of "low energies" is quantified by 
Eq. Q.) Restricting ourselves to three dimensions for 
concreteness, the relative wave function of two quan- 
tum particles colliding along the z-axis in a short range 
potential can be written as 



>/,„(£) = e 1 ^ + f(p,6) 



,ipr/h 



(10) 



The scattering amplitude f(p, 0) depends on the cen- 
ter of mass energy p 2 /m r (m x is the reduced mass) and 
the scattering angle 0. Solving the scattering problem 
for a particular potential U(r) consists in determining 
f(p, 0) or, equivalently, all partial wave scattering am- 
plitudes fi(p) in the expansion f{p,0) = J^. (2l + 
l)/;(p)P;(cos 0) with Legendre polynomials Pi. A non- 
vanishing relative angular momentum I of the scatter- 
ing particles introduces a centrifugal barrier term h 2 l(l + 
l)/2m r r 2 in the Schrodinger equation of relative motion. 
As a good estimate for the corresponding energy, we 
can replace r 2 — > ( 2 m and find that this barrier is far 
too high for particles with energies p 2 /2m r «: S 2 /^ dw . 
Therefore, only isotropic s-wave-scattering (/ = 0) oc- 
curs in ultracold alkali quantum gases. 

With p = hk, the low momentum expression of the 



Boettcher, Pawlowski, Diehl / Nuclear Physics B Proceedings Supplement 00 (2013) l- \73\ 



6 



s-wave scattering amplitude is given by 

/oO) = — WTr, — Z 

- l + r e k 2 /2 - lk + . 



(11) 



In this expansion, a is the scattering length anticipated 
above, and constitutes the most important parameter 
quantifying scattering in ultracold quantum gases in 
three dimensions. The coefficient r e is referred to as ef- 
fective range. It represents a correction which for the 
available k in ultracold gases is subleading, and thus 
we work with / -a. From Eq. (IOi we then have 
\)/{pt) a/r for large r and low momenta. 

The limitation to s-wave scattering has drastic conse- 
quences for ultracold gases of identical fermions. They 
are necessarily noninteracting. Collisions would only 
be possible in the p-wave or higher channels, but these 
cannot be reached due to the low energies. In order to 
have interactions between ultracold alkali fermions, we 
therefore always need at least two different species. 

From a low energy expansion of the s-wave scattering 
amplitude in one and two dimensions, respectively, it is 
possible to derive parameters similar to a which quan- 
tify scattering in reduced dimensionality. Such low di- 
mensional geometries can be designed in experiments 
by choosing a highly anisotropic harmonic potential 
with strong confinement in either one or two directions 

EQ. 

Let us briefly comment on the role of inelastic col- 
lisions. For collisions which do not change the spin of 
the particles, the most important inelastic mechanism is 
the formation of a molecule: If two atoms come close 
to each other, there may be energetically lower lying 
bound states and it is desirable for both atoms to build 
a molecule. However, without a third partner which al- 
lows for conservation of energy and momentum in this 
process, the excess energy from binding cannot go any- 
where. Therefore, in two-body processes molecule for- 
mation is ruled out. If a third atom is involved, we 
end up with a high kinetic energy of both the third 
atom and the molecule. These fast particles are then 
expelled from the trap. This three-body loss results in 
a finite lifetime of the gas. Due to diluteness and the 
contact interaction nature of ultracold atoms, such pro- 
cesses are suppressed and we find stable gases even at 
extremely low temperatures, where solidification would 
be expected. Increasing the density, we have to ensure 
the typical time scale of three-body recombination to be 
much larger than the experimental time of observation. 

Equipped with the length scale characterizing inter- 
actions, we give a concrete meaning to the notion of 
"weak" interactions by requiring the scattering length 



to be much smaller than the interparticle spacing in this 
case. This is equivalent to the gas parameter an 1 ^ be- 
ing small. The criterion for weak interactions 



an 



1/3 



« 1 



(12) 



is often referred to as diluteness condition. This in- 
terpretation is motivated by the fact that the scattering 
length provides the typical extent of a particle as far as 
its collisional properties are concerned. We indicate this 
in Fig. [T] In the dilute regime, it is possible to perform 
a perturbation theory in the gas parameter. 

For short range interaction potentials and low energy 
scattering, the s-wave scattering length can be calcu- 
lated within the Born approximation. It is then given 
by the Fourier transform of the interaction potential at 
zero wave vector (TJ, 



Anh 2 



s 



d 3 xU(r). 



(13) 



In particular, this formula can be applied to the 
Lennard-Jones potential for cold atoms introduced 
above. Importantly, from Eq. ( fl"3] l we learn (i) that 
value and sign of the scattering length may depend sen- 
sitively on the short distance physics of the interatomic 
potential and (ii) that we do not need to know these de- 
tails, since very different shapes of the interaction poten- 
tial will have the same scattering length, i.e. the same 
low energy scattering behavior. Quite remarkably, it is 
therefore possible to replace the microscopic Lennard- 
Jones potential by any other model potential producing 
the same scattering length (cf. Fig. [2]). For practical 
reasons, it is often convenient to work with completely 
local contact potentials 



U(r) = g A 6&. 



(14) 



This simple model potential needs regularization at 
short distances L = H/A, and a subsequent renormal- 
ization procedure. We remind to this fact with the index 
referring to an ultraviolet cutoff A, and refer to Sec. 4.2 
for a detailed presentation of the procedure. The cutoff- 
independent renormalized coupling constant g is related 
to the physically measured scattering length by the sim- 
ple formula 

m 

a - m* (15) 

The above considerations on ultracold atoms can be 
summarized in the effective Hamiltonian 

H = jT^OO ^~ + V«ttf)Wj& + Sa»(*) 2 ) , (16) 
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Scattering length 


Interparticle spacing 


de Broglie wavelength 


Trap size 


a/a 




At/Oq 




(0.05 ... 0.2) x 10 3 


(0.8 . . . 3) x 10 3 


(10...40)x 10 3 


(3 . . . 300) x 10 3 



Table 1: Standard scale hierarchy in ultracold quantum gases with typical values in units of the Bohr radius ciq = 0.53 X 10~'°m. The ratios of 
scales have the following physical meaning: ajt -SC 1 — weakly interacting or dilute; C/Aj « 1- ultracold. As long as t osc is the largest length 
scale, the local density approximation is valid (except for long distance physics in the vicinity of a critical point). 



t f t t 

interaction scale interparticle spacing de Broglie wavelength (Oscillator length of trap) 

i : i 

Feshbach resonance 

t or t 
optical lattice spacing 



Figure 3: Violations of the scale hierarchy, which do not invalidate the effective Hamiltonian, since all length scales are much larger than f vl jw 



where the operators aHx) and a(x) create and annihilate 
an atom at point x, respectively, and h(x) - a\x)a(x) is 
the local particle density operator. Note that the power 
of two in the interaction term ~ gh 2 stems from the fact 
that two particles have to meet at one point in order 
to interact. The trapping potential V ext (x) lifts the en- 
ergy of the particles at the point x and thus this term is 
proportional to h(x). This rather universal Hamiltonian 
provides an accurate description for all ultracold alkali 
atoms. 

It is a key feature of ultracold quantum gases that 
they are accurately described by effective microscopic 
Hamiltonians which depend only on a few system pa- 
rameters. The latter can be measured in experiments 
to a high precision, e.g. by spectroscopic methods or 
by measurement of the collisional cross sections |39], 
without the need to resolve the full interatomic poten- 
tials. This situation is very distinct from condensed mat- 
ter systems, where the underlying microscopic model is 
not known to such precision, and often has to be approx- 
imated by an educated guess. Moreover, realizations of 
ultracold quantum gases allow to change the system pa- 
rameters continuously and thus to understand their in- 
fluence on the many-body state. 

We summarize our discussion by indicating the stan- 
dard scale hierarchy in Tab. [Tj which is built up from 
the scattering length a, the interparticle distance (den- 
sity) (, the thermal wavelength (temperature) At and 
the oscillator length i osc . Moreover, the system has a 
natural UV cutoff A -1 <k £ V dw- Microscopic details on 



shorter length scales are irrelevant for our purposes be- 
cause none of the many-body length scales can resolve 
the underlying physics. 

It is both experimentally and theoretically appealing 
that ultracold atoms can be tuned such that they violate 
the scale hierarchy, allowing to reach strongly interact- 
ing regimes - crucially, without loosing the validity of 
the above discussion. One way is provided by Feshbach 
resonances of the scattering length [39|. Here, we can 
loosen the condition a <K t and explore new regimes 
of the many-body system which are not captured by 
mean field theory or perturbative expansions. Such res- 
onances are realized in cold atoms if a bound state is lo- 
cated close to the zero energy scattering threshold, and 
is tuned to resonance due to the variation of an external 
magnetic field. From this we infer that a Feshbach res- 
onance is a result of a specific fine-tuning of the micro- 
physics. The scattering length can then be parametrized 
according to 



where a^g, AB and B$ are background scattering length, 
width and position of the resonance, respectively. Ob- 
viously, approaching Bq, we can obtain an anomalously 
large scattering length, meaning that, by virtue of fine- 
tuned microphysics, it greatly exceeds the generic scale 
set by the van der Waals length, 

\a\ » { vdw . (18) 



Boettcher, Pawlowski, Diehl / Nuclear Physics B Proceedings Supplement 00 (2013) l- \73\ 



8 



We will discuss the issue of Feshbach resonances in 



more detail in Sec. 4.3 when introducing our effective 
model for the BCS-BEC crossover. 

Another way to reach an interaction dominated 
regime is by superimposing an optical lattice [3] 01. 
This is a standing wave of counterpropagating laser 
beams in each spatial direction, which provides a con- 
servative periodic potential landscape for the atoms via 
the AC Stark effect. Tuning the depth of the lattice wells 
via the laser intensity, we can withdraw the kinetic en- 
ergy more strongly then the interaction energy and thus 
arrive at a strongly correlated system. The lattice spac- 
ing is related to the wavelength of the light used for the 
optical lattice. By engineering neighboring sites close 
to each other, we can reach high densities ("fillings"). 
Each of these effects enhances the correlations in the 
system. 

Recall that the validity of the effective Hamiltonian in 
Eq. ( fT6] > is restricted to length scales sufficiently larger 
than £ v dw- Since the mentioned scale violations happen 
at larger scales, the faithful microscopic modeling is not 
touched. Therefore, the pointlike description of the in- 
teractions is also applicable in the dense and strong cou- 
pling regimes. 

For a more detailed presentation of low energy uni- 
versality in atomic few-body systems and from the 
viewpoint of quantum field theory, we refer to l42l . 

2.2. Thermodynamics 

In this section we review a few thermodynamic concepts 
which are of relevance for experiments with ultracold 
atoms. We derive general thermodynamic statements, 
which hold independently of the particular system un- 
der consideration. We will see that the phase diagram 
and the equation of state encode important, experimen- 
tally accessible information about a system and thus are 
desirable quantities to be computed from first principles. 
This also serves as one of the motivations to investigate 
cold atoms with the Functional Renormalization Group. 

For thermodynamics to be applicable, we require the 
internal processes of a many-body system to be such 
that the system is in equilibrium on the time scale of ob- 
servation. Strictly speaking, thermodynamic statements 
and, in particular, the theory of phase transitions are 
only valid in infinitely large systems. But this require- 
ment is less severe as it might seem at first sight, be- 
cause any thermodynamic relation can be expressed in 
terms of intensive quantities only, like particle density, 
entropy density, or magnetization per particle. Taking 
these densities to be local quantities, we can apply the 
laws of thermodynamics locally for small subsystems 
of finite volume and particle number. This procedure 



works perfectly at room temperature with large particle 
numbers N ~ 10 23 , and is still justified for trapped gases 
with typically N ~ 10 4 — 10 7 . Let us add that in addi- 
tion, it turns out that systems with low atom loss rate 
and long lifetime can indeed be assumed to be thermo- 
dynamically equilibrated over the period of observation. 
Such a system is provided by two-component fermions 
in the BCS-BEC crossover. 

We recall that the full thermodynamic information of 
a system is stored in the equation of state P(p, T), which 
can be expressed in terms of the pressure as a function of 
chemical potential and temperature. Using the Gibbs- 
Duhem relations dP = ndjj+sdT ande = Ts—P+fin we 
can calculate all other intensive thermodynamic quan- 
tities from the pressure. Here, n = N/V, s — S /V 
and e - E/V are the densities of particle number, en- 
tropy and energy, respectively. The chemical potential 
/i is a parameter which determines the particle number 
N(jS) for a given temperature. Eliminating p for the 
density n(p, T), the equation of state can also be formu- 
lated in terms of the free energy density /(«, T), which 
is the Legendre transform of the pressure according to 
f = fm-P. 

In order to understand the influence of a trap, we con- 
sider a cloud in a time-independent external potential 
V ex t(x) which varies on much larger length scales than 
the typical atomic ones (e.g. interparticle spacing and 
scattering length). Picking two neighboring small but 
yet macroscopic subvolumes V\ and V2 of the cloud, 
thermal and chemical processes between them will re- 
sult in the equality of their temperature and chemical 
potential. Since the subvolumes were arbitrary, we con- 
clude that temperature and full chemical potential are 
constant inside the trap. However, from the Gibbs- 
Duhem relation we infer that the full chemical poten- 
tial corresponds to the Gibbs free energy G — F + PV 
per particle, p = G/N. The latter is spatially inhomoge- 
neous due to the trap and we find p = ph om (n(x), T) + 
V ex t(x) = const. In this formula, phom(n, T) is the chem- 
ical potential obtained from a calculation in a homoge- 
neous setting, e.g. a box of volume V containing N par- 
ticles. 

We conclude that a system where the thermodynamic 
quantities are replaced according to 



P(p,T)^ P(p-VUAT) 



(19) 



behaves like a system trapped in a potential of large 
spatial extent. This prescription is called local den- 
sity approximation (LDA). The above derivation pro- 
vides an intuitive understanding why this procedure 
should give reasonable results. Of course, if we can- 
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Figure 4: Within LDA, the inner regions of the trapped cloud corre- 
spond to higher chemical potentials: p\oc(r) = p. — V cxl (r). We show 
here the density profile of a weakly interacting Bose gas in an external 
harmonic confinement. The inset shows a typical phase diagram in the 
(p, T)-plane, where the blue region represents the superfluid phase. 
We cross the critical line of the superfluid phase transition for fixed 
temperature T at a certain chemical potential pAT)- This corresponds 
to a kink in the density profile at a critical value n c (T). Note, however, 
that LDA breaks down in the outer regions of the cloud, where the gas 
is extremely dilute. 

not pick small, yet macroscopic subvolumes, the argu- 
ment breaks down. The applicability of LDA is there- 
fore limited to systems, where the trap { osc provides the 
largest length scale. This agrees with our earlier con- 
siderations. From a field theory perspective it is very 
promising that properties of homogeneous systems can 
be obtained from trapped gases and, indeed, there have 
already been beautiful measurements of the equation of 
state of the BCS-BEC crossover to a high precision us- 
ingLDA I43H451. 

The equation of state also contains information about 
possible phase transitions appearing in the many-body 
system. Phases consist of extended parameter regimes 
which can be distinguished from each other by macro- 
scopic observables. As an example, we consider the 
element iron. Despite the difference of solid, liquid 
and gaseous phase we can independently also distin- 
guish the ferromagnetically from the antiferromagnet- 
ically ordered phase, or furthermore the crystal struc- 
tures of y-Fe and a-Fe. Phase transitions manifest them- 
selves through kinks and jumps in the thermodynamic 
functions, typically in the higher derivatives of P(p, T). 
These root in non-analyticities contained in the parti- 
tion function. It is easy to see that true phase transitions 
need a continuum of degrees of freedom, i.e. occur only 
in the thermodynamic limit. Indeed, the partition func- 
tion is Z = Tre~P H = 2„e^ £ ", E„ the eigenenergies 
of the system. Each of the contributions is analytic. 
Non-analyticities can only be generated in the case of 
infinitely many states entering the sum. 



9 

More formally, we distinguish two phases by an or- 
der parameter po(jx, T), which depends on the thermo- 
dynamic variables. In different phases, it is either zero 
or nonzero, which gives rises to the phase diagram in 
the (p, T)-plane. For a fixed value of the chemical po- 
tential, we define the critical temperature T c {p) via the 
relation po(p-, T c {p)) = 0. Of course, we can also fix 
the density n to obtain the critical temperature T c (n) as 
a function of n. 

In the regime where LDA is applicable, the local 
chemical potential p\ oc (x) = p. - V ext (x) has its largest 
value at the minima of the trapping potential. Accord- 
ingly, an increase of the potential reduces p\ 0Q . For this 
reason, we can scan the phase diagram over a certain 
region from a density image in a harmonic potential 
V e xt(x) = f w o r2 ' see Fig- From our above consid- 
erations we conclude that the corresponding path in the 
(p, J)-plane is an isothermal line. In particular, we may 
cross the phase boundary when the local chemical po- 
tential reaches the critical value p c (T). For this reason, 
we can have a superfluid gas in the inner regions of the 
cloud, whereas the outer shell is in its normal phase. 
The lobes in the phase diagram of the Bose-Hubbard 
model lead to a wedding cake structure of the density 
profile 1461 . 

2.3. Noninteracting Bose and Fermi gases 

After these general remarks on thermodynamics we turn 
our attention to degenerate, noninteracting Bose and 
Fermi gases. The notions of Bose-Einstein condensa- 
tion and Fermi surfaces are introduced. They constitute 
the two cornerstones of quantum statistical phenomena 
and are crucial for understanding interacting gases. 

The state of a single particle can be addressed by its 
momentum p and spin-projection cr. The correspond- 
ing occupation numbers n^ are restricted to 0, 1 for 
fermions due to Pauli's principle, whereas they can have 
arbitrary integer values 0,1,2,... for bosons. As is 
known from statistical mechanics, we then find for the 
equation of state 

/J if 
^-^log(l+ e -« £ ^), (20) 

where s p = p 2 /2m and g is the spin degeneracy of the 
momentum states. We have g = 1 for spinless bosons 
considered here, and g — 2 for spin- 1/2 fermions. 
The upper (lower) sign in Eq. ( f20| ) holds for bosons 
(fermions). As we will see below, for bosons, this ex- 
pression is valid in the absence of a condensate only. 
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Free bosons and Bose-Einstein condensation 



p,(n, T) 



To understand the appearance of condensation as a 
purely quantum statistical effect, we consider an ideal 
gas of identical bosons. At zero temperature, we ex- 
pect all bosons to be in the single particle state with 
energy e — 0. In particular, this means that the occu- 
pation number iVo °f that state is extensive, No ~ V. We 
say that the zero mode e = is occupied macroscopi- 
cally. At low nonzero temperature, some particles will 
be thermally excited into the higher states. At very high 
temperatures, we approach the Boltzmann limit, where 
all occupation numbers are small (in particular, none of 
them is occupied macroscopically) and the distribution 
function n(s) is very broad. Therefore, there must be a 
critical temperature T c below which macroscopic occu- 
pation of the single particle ground state sets in. Since 
this particular behavior is due to quantum statistics and 
absent in a classical gas, we can estimate the critical 
temperature very roughly to occur for A Tc =t {, 



Starting from Eq. ( 20 1 for the pressure, we obtain the 
particle number in a three-dimensional box of volume 
V by virtue of a /i-derivative as to be given by 



N(T,V,fi)=Yj<^ = 



V 1 



Jo 



de Ve 

4T(3/2) Jo 



(21) 

For fixed temperature and volume, this formula has a 
maximum N mix at p = 0. However, if we decide to 
put more than N max particles into the box, the expres- 
sion necessarily becomes invalid. The critical tempera- 
ture T c (n) where this happens is determined by a critical 
phase space density 



= «4. = (fajt? = £(3/2) 2.612, 



(22) 



i.e. Ajjl = 0(1) as anticipated above. 

Since our starting point was physically sound, but we 
ended up with an unphysical result, we must have made 
an error. This led Einstein and Bose to treating the zero 
momentum mode separately B71 l48ll . Indeed, in Eq. 
( 21 1 we did not appropriately incorporate the states with 
s = 0: Replacing the quantized momenta of the finite 
system = 2nhnjL in the naive continuum limit 

we multiply their contribution with s = (or equiva- 
lently p 2 = 0). This corresponds to a vanishing occupa- 
tion of the single particle ground state, which constitutes 
a bad approximation, as is apparent from our above con- 
siderations. 




>d = 2 



Figure 5: We plot the chemical potential p(n, T) from Eq. (20) for 
an ideal Bose gas at fixed density n. In three spatial dimensions, the 
function hits zero at T c (ri) > 0. Since Eq. (20) cannot be applied for 
positive p, the chemical potential remains zero and condensation sets 
in. In contrast, the chemical potential in two dimensions is negative 
for T > and thus there are always enough thermally excited states 
and condensation is absent. 



Therefore, the situation at temperatures T < T c (ri) 
is as follows. Formula pi) with fi = describes the 
excited particles in the states with s > 0. The remain- 
ing Nq(T) = N — N ex {T) particles are condensed to the 
zero energy state, leading to its macroscopic occupa- 
tion. This resolves the puzzle from above. If we put 
more than A^ max particles into the box, they will add to 
the condensate. The particle number below T c is given 
by 

^(r,v) = <aJajf> + ^(3/2). (24) 

Obviously, No(T) = (flijflg) ~ V is extensive. The con- 
densate fraction Nq(T)/N is an order parameter for the 
Bose-Einstein condensation phase transition. From Eq. 
d22l we conclude 



Nq(T) 
N 



1 



3/2 



for T < T c 



(25) 



It vanishes continuously for T — > T c , which signals a 
second order phase transition. 

In Eq. ( |23| ) we used the three-dimensional density of 
states p(e) ~ Ve to show why condensation appears. 
In d spatial dimension, we have p(s) ~ e^ 2-1 and the 
ground state contribution is not multiplied by zero for 
d < 2. Indeed, a similar calculation shows that for 
one- and two-dimensional systems the particle number 
N(T, V,fi) does not have a maximum at nonzero temper- 
atures and thus Bose-Einstein condensation is absent. 
For d = 1 this also holds at zero temperature. In Fig. 
[5] we plot the chemical potential as a function of tem- 
perature. Whereas p(T c ) = for a nonzero T c in three 
dimensions, we find T c — for d — 2. Our finding for 
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noninteracting particles is a special case of the gener- 
ally valid Mermin-Wagner theorem ll49l . which states 
that there is no spontaneous breaking of a continuous 
symmetry in d < 2 (noncompact) dimensions. The in- 
gredients to this theorem are the locality of the under- 
lying Hamiltonian, and the universal relativistic long- 
wavelength form of the dispersion relation. The long- 
range order is then destroyed by fluctuations with very 
long wavelengths. However, in atomic gas experiments 
the trap provides the largest length scale i osc , such that 
these fluctuations are not present and condensation can 
be observed in lower-dimensional geometries. 



Free fermions and Fermi surface 

Whereas the appearance of a Bose-Einstein condensate 
is closely related to the fact that identical bosons can 
have arbitrarily large occupation numbers, the notion of 
a Fermi surface is a consequence of Pauli's principle for 
many-fermion systems. 

To get an intuition, we consider an ideal gas of N 
identical fermions. What will be the ground state of 
the quantum many-body system? (This state is realized 
at zero temperature.) Obviously, each of the particles 
seeks to minimize its energy. But since every single par- 
ticle state can only be occupied by at most one fermion, 
the ground state will be such that precisely the N ener- 
getically lowest lying states are occupied. Equivalently, 
due to rotation symmetry, all states with momenta in- 
side a sphere of radius pp in momentum space will be 
occupied. Restricting to three dimensions, we can count 
states by dividing the classical phase space into cells of 
volume h 3 . This yields 



N 



gV 4n 



"ft" 



(26) 



(27th) 3 3 

We call p P = Hkp the Fermi momentum and deduce 

k F (n) = (6n 2 n/g) 1/3 . (27) 
Up to a prefactor of order unity, kp(ri) equals the in- 



verse interparticle spacing (~ = n 1 . Eq. (27i can 
be used to express thermodynamic quantities as a func- 
tion of kp(n) instead of the density n. In this case, 
kp(n) is not bound to the presence of a Fermi surface. 
From pp we construct the Fermi energy and tempera- 
ture, sp = e Pf = p F /2m and Tp = sp/k^, respectively. 
Our simple picture of the many-body ground state is 



k B T 

< > 



£f = k B T F 



e/ef 



Figure 6: The Fermi-Dirac distribution at zero temperature (solid 
line) constitutes a step function. For T > (dashed line), broaden- 
ing appears in a region of width T around the Fermi edge located 
at e = £p = &b7*f. If the distance of the edge from the origin is much 
larger than the area of broadening, the distribution function still dis- 
plays the characteristic step-like behavior. Clearly, this picture is valid 
for the dimensionless parameter T/Tp being small. 



correct, because from Eq. (20 1 we have 



- — 
dp 

r->o 



r d 3 p 

J (2nK 



1 



H) 3 eP^t-A + 1 
0(e p -p) 



6n 2 h? 



(2mp) V2 . (28) 



On the other hand, from the zero temperature limit of 
the Fermi-Dirac distribution we infer that the highest 
energy present in the system is e F = p and thus we 
find pp = (2m//) 1 / 2 . As before, we finally arrive at 
kp = {6n 2 n/ g) 1 ^ 3 . The Fermi-Dirac distribution at zero 
temperature is shown in Fig. [6] 

What happens to this picture at nonzero temperature? 
The Fermi-Dirac distribution n F (s) = ( e {e -> 1)lk * T + 
is no longer a sharp step function but rather smears out 
around e - p. Nevertheless, the smeared out region 
is of order k^T, whereas the distance of the edge from 
e = is of order p =t e F = kp,Tp. Therefore, as long 
as T/Tp <k 1, the distribution function looks approxi- 
mately like a step function. We visualize this situation 
in Fig. [6] For T > 0, there are thermally excited par- 
ticles with energies close to the chemical potential. We 
conclude that the low lying excitations of a Fermi gas 
are not at zero momentum but rather at momenta close 
to the Fermi surface, which consists of the momenta 
\P\ = Pf- 
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3. Functional methods for interacting bosons and 
fermions 

3.1. Functional integral and effective action 

In this section, the quantum field theoretical formula- 
tion of interacting cold atoms is put forward. Starting 
from the functional integral representation of the quan- 
tum partition function Z, we introduce the effective ac- 
tion F |49|. The latter stores the same information as 
the partition function, however in a way that is more in- 
tuitive. In particular, it naturally provides the classical 
limit. The effects of both quantum and thermal fluctua- 
tions on physical observables can be derived from it in 
the few- and many-body context. Moreover, it allows 
for a transparent discussion of spontaneous symmetry 
breaking, and allows to leverage the power of symmetry 
considerations from the classical action over to the full 
quantum theory. We set h = k s = 1 . For the moment, 
we keep the nonrelativistic mass M in our formulation, 
but later we will set 2M = 1 in the same spirit as for the 
fundamental constants. 



Functional integral 



As we have shown in section 2.1 a system of ultracold 



atoms is accurately described by the effective Hamilto- 
nian 



„,(.?> I «(.?> + |S«n, (29) 



where a* (x) and a(x) are operators which create and an- 
nihilate an atom at position x, respectively. Depending 
on whether we consider bosons or fermions, these op- 
erators satisfy commutation or anti-commutation rela- 
tions. The density operator is given by h(x) = a 1 (x)a(x). 

The Hamiltonian in Eq. p9| ) defines a quantum field 
theory with operators a and iW on each point of space. 
Physical observables are derived from expectation val- 
ues of functions of these operators. However, the corre- 
sponding quantum field theory can also be formulated 
in terms of a functional integral. The latter does no 
longer depend on the notion of operators. In the context 
of quantum many-body systems, a possible derivation 
starts from the grand canonical partition function 



Z( M ,T) = Tr(e-K A -^), 



(30) 



where the trace is taken over Fock space. This trace can 
be represented in the basis of so-called coherent states, 



which are eigenstates of the annihilation operator a(x). 
We then obtain 



ZQx,T) 



j Dtp*D 



ipe 



-S [if", if] 



(31) 



The expression in Eq. ( 3 1 1 is called a functional 
or path integral. It contains the microscopic action 
S [(p* , <p] of a field theory, which in our case is nonrela- 
tivistic. It is related to the normal orderecQHamiltonian 
H = Hla* (x),a(x)\ according to 

S[y*,<p]= f dr f <p*(T,x)(d r - fi)<p(r,x) 
Jo Jj? 

i 

dTH[<p*(T,x),<p(T,x)]. (32) 



f 

Jo 



For the particular choice of the effective Hamiltonian in 
Eq. (|29|, we have 



+ ^(<P*(t,x) V (t,x)) 2 ). (33) 

The explicit construction of the functional integral 
representation of the partition function for a generic 



many-body Hamiltonian is carried out in Appendix A 
We summarize here the two main findings. 

1) Bosonic atoms are represented by complex fields 
<p(T, x), whereas fermions are described in terms of 
Grassmann valued fields i//(t, x). 

2) The non-commutativity of operators introduces the 
imaginary time t, which is restricted to the inter- 
val [0,/?]. Bosonic fields are /3-periodic in time, 
i.e. (p(J3, x) - ip(0, x). In contrast, fermionic fields 
satisfy ifr(/3, x) = -iff(0, x). 

The second property implies that the Fourier transform 
of the fields tp and ip in imaginary time direction reduces 
to a Fourier series with discrete Matsubara frequencies 



2nnT (bosons) 
oj„ = \ , n e Z. (34) 

\2n{n + l/2)T (fermions) 



We say that the imaginary time direction is compacti- 
fied to a torus of circumference j3. In fact, introducing 



' Normal ordering is a prerequisite to the construction of the coher- 
ent state functional integral, cf. |Appendix A| Starting from Eq. \29) , 
normal ordering introduces a shift /j — > p - g/2, which we simply 
absorb into a redefinition of the chemical potential. 
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Ising magnets 


Bosonic atoms 


lattice sites x, 


space-time points X = (t, x) 


magnetic moment m, ■ = ± 1 


complex field ip(X) 


magnetic field /i, 


external source j{X) 


mean field m, = (m,)/, 


cf>(X) = (<p(X))j 


partial derivative 


functional derivative 


summation over sites 2, 


functional integral fDtp'Dtp 


functional measure n, dm, 


= n* ^ { T X) 


partition function Z((A,j) 


generating functional Z[f, j] 


dlogz 
m >'= a*, 


</>(*) = wm 


{m imj ) h tmrnj = m . d l. 


(<P x (p Y )j <f> x <f>Y = gjjj* 


effective action r({m,() 




h> = ^- 

1 dm. 


- 64,* {X) 



Table 2: Correspondence of quantities for classical Ising magnets and 
cold atoms. The magnet case on a lattice is discussed in | Appendix B| 



a chemical potential fi and compactifying the time di- 
rection to < r < we can describe any euclidean 
quantum field theory at nonzero density and temper- 
ature. For example, this procedure can be applied to 
quantum chromodynamics. 



Generating functional and effective action 

Starting from the functional integral representation for 
the partition function Z, we now construct the corre- 
sponding effective action. The procedure outlined here 
focuses on the application to systems of ultracold atoms. 
However, additional insights into these concepts can be 
obtained from a comparison to classical Ising magnets 
on a discrete lattice. For this reason, we included a de- 



tailed discussion of the latter system in Appendix B 
There, we also perform the continuum limit and re- 
view the notions of functional differentiation and inte- 
gration. A dictionary for the translation between ul- 
tracold bosonic atoms and the classical Ising model is 
given in Tab. [2] For simplicity, we mostly restrict to 
the bosonic case in this section. There are only minor 
modifications for fermions, which are discussed at the 
end of the section. 

The bosonic functional integral in Eq. pT} allows 
for the definition of a probability measure on the set of 



fields ip. Given an observable 0({<p}) which depends on 
the field, we define 

(O) = 1 j D V *D<pO({<p})e- sl *'* ] . (35) 

Herein, the action S acts as a weight. For instance, from 
O = <p(X) or (9 = qf(X)(p(Y) we obtain the one- and (dis- 
connected) two-point correlation functions of the the- 
ory. More generally, we obtain averages of observables 
by introducing a complex source field j(X) according to 

W = / V D,,**W* (36, 



We call Z[j*,j] the generating functional and have 

51ogZ 



4>{X) = (<p(X))j 

(<p*(X)<p(Y))j 



sf(xy 

1 6 2 Z 



(37) 
(38) 



Z6j(X)6f(Y)' 
etc. The subscript j indicates that the external source is 



not yet set to zero. Generalizing Eqs. ( 37 i and ( 38 1, we 
find the representation of a general expectation value 



(39) 



We conclude that all correlation functions of interest can 
be obtained from the generating functional Z or 



W[j\j]=logZ[f,jl 



(40) 



The latter quantity is called the Schwinger functional. 
It constitutes the generating functional of connected n- 
point functions. For instance, we find for the connected 
two-point function 



W (2) [j](X,Y) = 



6 Z W 



6f(X)6j(Y) 
= ((p(X)ip*(Y))j - <f>(X)f(Y) 

= <^0^(r»j, c> 



(41) 



where the index refers to "connected". This object is 
also called (time ordered) Green's function or propaga- 
tor of the theory. Imposing the time-ordering for the 
propagators leads to time-ordered general correlation 
functions ([39]). 

The field expectation value carries a direct physical 
significance. For example, in a homogeneous situation, 
(f>(X) = 0o describes the Bose-Einstein condensate. It 
therefore seems desirable to implement it into the the- 
ory in a more direct way. In fact, by the aid of a Leg- 
endre transformation, this is possible and gives rise to 
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the effective action, see e.g. [49]. Similar to our con- 
siderations for the Ising magnets in Appendix B we 



introduce the latter as the generating functional which 
depends on the mean or classical field (f> defined by 



cf>(X) = {ip(X))j = 



6W 

SfVO' 



(42) 



Assume we have solved this equation. We can then 
construct the effective action according to the Legendre 
transformation 



= f (fx 

Jx 



*xjx + fx<t>x)-W[f,jl (43) 



where j and j* are defined implicitly by Eq. (42 1. Note 



that while the active variable for the partition function 
is the source, Z = Z[j], the active variable for the effec- 
tive action is the field expectation value, T = Y[(f>]. The 
effective action is thus parametrized directly in terms of 
a physical observable. Applying the chain rule for func- 
tional differentiation we find 



6F 



[f ,</>]= f(X). 



(44) 



For the derivation of this relation in the case of discrete 



variables, see Eq. (B.9i 



Technically speaking, the effective action is the gen- 
erating functional of one-particle irreducible (1PI) cor- 
relation functions. They can be obtained from F by 
taking successive functional derivatives with respect to 
<p(X) and (p*(X). Diagrammatically the 1PI correlation 
functions are given by all diagrams that cannot be split 
by cutting one (internal) line, hence the name. Physi- 
cally, such inhomogeneous mean fields can be obtained 
by applying external sources, namely by choosing j 
such that 5YI5(f) = j. Often, we are mainly interested 
in the situation of vanishing source. Then, given the ef- 
fective action, we have to solve the equations of motion 



6T 

mx) 



(45) 



to obtain the thermodynamic equilibrium state 0o - 
((p)j=o of the theory. The reference to the external field 
is no longer present and also not needed, because it is 
already included in r[(p*,<f>]. Typically, the solution </>q 
to Eq. (45 i is constant in space-time. This will be ex- 



plained in more detail below Eq. ( 56 1. However, in gen- 
eral there also might be inhomogeneous solutions (f>o to 
the nonlinear partial differential equation Eq. ( j45j ), such 
as instantons, solitons or vortices (see Sec. |3.2| i. 

Higher derivatives of the effective action with respect 
to the fields <p, <p*, denoted with r (n) for the nth deriva- 
tive, provide the one-particle-irreducible vertices. The 



second derivative of the effective action, 

6 2 T 



T (2) (X, Y) = 



8<p*{X)5<t>{Y)' 



(46) 



plays a special role, as it is the inverse propagator. This 
is easily proven by 



I 



r (2) (x,z)w (2) (z,r> 



6j*{Z) 5f(Y) + 6j(Z) 6<P*(Y)\ 



L[sf(X) 6HZ) S<f>*(X) 6j(Z) 
= 6(X - Y) , (47) 



where we have used (|42[>,(|44[> and the completeness re- 
lation of derivatives with respect to /' , /'*. 



In principle, Eq. ( 44 1 can be taken as a starting point 
to calculate the effective action in certain ap- 

proximations. However, the definition of F implies an 



exact identity, which is equivalent to Eq. (44 1, but more 
useful. Applying Eqs. (43 i, (44i and W = logZ, we 
arrive at 

e -WS _ e -J x (f"Wf)+w 

= J Dip* D^e _s[ ^* 1+ /xO'*(^)+(¥>*-^*)» 

= J Dtp* D(p e~ S [ * ! ** l+ ir ( I W<f-<t>)Hf*-4>'> [(*])_ (48) 

This equation is called the background field identity for 
the effective action. For (p = <po with tpo satisfying Eq. 
d45l, we recover 



ZQx,T) = e 



(49) 



i.e. the effective action then corresponds to the grand 
canonical potential. Furthermore, by performing a shift 



of the integration variable, we rewrite Eq. (48 i as 

D5<p*D6ye~ s % [#<V+^* •■§? M)_ 



(50) 



This functional integral representation of the effective 
action gives rise to the intuitive picture that the effec- 
tive action encodes the complete information on the eu- 
clidean field theory by means of summing over all pos- 
sible field configurations dip deviating from the classical 
one, 0. 
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We now show that in the classical limit, the effective 
action and the classical action coincide. Reintroducing 
Planck's constant h, we have F/h and S/H appearing in 
Eq. (50 1. The classical limit is obtained for h — > at 



fixed F. The integrand is then sharply peaked around the 
solution to the classical equations of motion 5S /8<p = 0. 
This results in F = S , which physically is the classical 
approximation. 

It is clear that the effective action lends itself ideally 
for semiclassical approximations, and also systematic 
improvements thereon. From Eqs. (48 1 or (50i, we 



can easily go one step beyond the classical approxima- 
tion, by expanding the exponent in the functional inte- 
gral around its minimum value <po determined by 

6S 6F SS 6T 

dtp of dip Of 

For the particularly simple case where (f> = <po, the linear 
derivatives cancel and we obtain 

e -r\m\ ^ e s\m\ J u^^-i f^'>^M{^) (52 ) 

with S <2) the second functional derivative of S with re- 
spect to ip, <p* . 



More generally, with the help of formula (B.19i, the 



Gaussian approximation to Eq. (50 1 can be evaluated at 
any field (p which ensures the path integral to be dom- 
inated by small fluctuations 5ip. This leads to the so- 
called one-loop formula 



T[f,f] *S[f*,f] + iTrlogS (2) [0*,# 



(53) 



In this order of approximation, the linear derivative 
terms cancel due to the tree level relation r - S . Note 
that the effective action equals the classical action also 
in the case of a free, noninteracting theory. Expanding 
the Tr log expression in powers of the field, we generate 
one-loop perturbation theory. We may therefore expect 



Eq. (53 1 to give good results in the perturbative regime 
of small coupling. 

Our considerations can easily be extended to 
fermions as well. We introduce independent Grassman- 
nian source terms t](X) and if (X) into the generating 
functional Z[j* , j,rf ,77], which couple linearly to the 
fields iff*(X) and i[f(X), respectively. The effective action 
is defined in the same manner as before via the Leg- 
endre transformation of log Z with respect to the mean 
fields. The ground state of the theory necessarily satis- 
fies (iff(X)), 1= {) = (if/*(X)) n= o = 0, since Pauli's princi- 
ple forbids macroscopic occupation of fermionic states. 
However, the generating functional T depends on non- 
vanishing fermionic "mean fields". Such fields \jj(X) can 



be constructed by applying a source rj(X) = 6F/Stff(X). 
They must not be regarded as physical objects, but 
rather as bookkeeping parameters used to generate the 
1PI correlation functions via Grassmannian functional 
differentiation. 



Eq. (53 1 is also valid for fermionic fields, but with an 



additional minus sign in front of the trace. For a mixed 
theory of both bosons and fermions we introduce the 
so-called supertrace, STr, which takes into account this 
sign for fermionic terms. Thus, we arrive at the one- 
loop formula 

r[0,<A] ^S[fM + isTrlogS (2) [«^]. (54) 

More detailed presentations on functional integrals 
can be found in I5t)ll5ll. 



3.2. Effective potential and spontaneous symmetry 
breaking 

In this section, we discuss how phase transitions and 
spontaneous symmetry breaking (SSB) find their natural 
description in terms of the effective potential U(p). The 
latter is the part of the effective action which does not 
contain derivatives of the field. It includes both quan- 
tum and thermal fluctuations, and typically changes its 
shape by tuning the system parameters like temperature, 
chemical potential or interaction strength. For parame- 
ter regions where the minimum of the effective potential 
is nonzero, small perturbations can drive the system into 
an equilibrium ground state which does not respect the 
symmetry of the underlying physical theory. The sym- 
metry is spontaneously broken. We exemplify this im- 
portant concept of many-body physics on systems with 
Z2- and U(l)-symmetry, respectively. 

An intuitive picture of SSB is provided by a simple 
daily life observation. Suppose a pencil is balanced on 
its tip to stand upright. Due to the cylindrical sym- 
metry, the pencil should stay in this position. Indeed, 
the underlying physics, here given by the gravitational 
force pointing downwards, does not prefer any direc- 
tion. However, if there is a small perturbation of this 
symmetry due to the environment, the pencil will imme- 
diately fall to the side and thereby minimize its energy. 
Even if the perturbation is removed now, the pencil will 
remain in the horizontal position. 



Thermodynamics from the effective action 

In order to study the properties of the thermodynamic 
equilibrium state, we consider a system of bosons. We 
assume the trapping potential V e!lt (x) to vanish and the 
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external source to be constant, j(X) = j. Hence, the set- 
ting is homogeneous in space-time. We learned in Eq. 
( |49| ) that the grand canonical partition function Zip, T) 
is related to the effective action according to 



r[0o] = -iogzoi,r). 



(55) 



Herein, the field expectation value <po(X) minimizes the 



effective action, as can be seen from Eqs. ( |43] l and (45 1 
The effective action F[0] has the structure 



T[0] = J (terms containing derivatives) + J U((f>(X)) . 

(56) 

If the part containing derivatives can be expanded in or- 
ders of the derivatives, it is necessarily non-negative for 
the sake of stability. Otherwise, F[</>] would not possess 
a minimum, because we could arbitrarily decrease its 
value by creating heavily oscillating fields. Then, (f>o is 
a constant field which additionally minimizes the effec- 
tive potential Ui<p) according to 



t/(0o) = min[t/(0)]. 



(57) 



Since the effective potential depends on both the exter- 
nal parameters /j and T, the same will be true for the 
field expectation value O = 4>oip, T). In the presence of 
a nonvanishing background source field, we also have 
an explicit dependence on j. Note that the above ar- 
gument does not exclude the existence of inhomoge- 
neous ground states as they cannot be expanded in terms 
of derivatives. The existence of such inhomogeneous 
ground states is common in low dimensions, in partic- 
ular in 1 + 1 dimensions, see e.g. [52 1 for the class of 
models under discussion here. 

Using Eq. (j49]), the effective potential at its minimum 
value is related to the pressure according to 



P(jx, T) = -U(4>o,n,T). 



(58) 



This constitutes the equation of state of the system. 
Often, we are mainly interested in the density nip,, T), 
which is found from dP = ndp + sdT. The relevant 
thermodynamic information contained in the effective 
potential can thus be summarized in the two equations 



dU 

dU 
dp 



i<po,p, T)-0 (gap equation), (59) 
((j)<j,H, T) - -n (equation of state). (60) 



These equations are generally valid and constitute the 
main building blocks for the evaluation of the phase 



diagram of the many-body problem. In particular, the 
above discussion is not limited to bosons, but can be 
applied to an arbitrary many-body system or quantum 
field theory, since the effective action approach is appli- 
cable to all of these system. For instance, as outlined 
earlier, the field <p(X) may as well describe the degrees 
of freedom in a Heisenberg ferromagnet with magnetic 
moments m, and m(x) on a lattice or in the continuum, 
respectively. However, we keep denoting the fields by <p 
and <f>, which may have to be replaced appropriately. 



Spontaneous symmetry breaking 

As a preparation for the more formal discussion of SSB, 
we first relate symmetries of the microscopic action to 
those of the effective action. To this end, we recall the 
definition of the effective action to be 



= jTXpe~ sw+ ^ MH,p ~^. (6i) 



Setting the external source j to zero, we see that any 
symmetry of the microscopic action which is respected 
by the functional measure, will also be a symmetry of 
the effective action. A nonvanishing source j(X), in- 
stead, typically leads to terms in the effective action 
which explicitly break the microscopic symmetry. This 
is accompanied by a nonzero expectation value <p{X), 
because j(X) either introduces a nonhomogeneity in 
space-time or at least singles out a direction in field 
space (p. 

Spontaneous breaking of a symmetry refers to a dif- 
ferent scenario. In this case, the external source van- 
ishes such that the effective action manifestly shares 
the symmetry of the microscopic action. Nevertheless, 
the ground state of the theory (or, more generally, the 
thermodynamic equilibrium state), may spontaneously 
break this symmetry due to a nonzero expectation value 
according to 

<h = * 0. (62) 

The symmetry is then broken because the field expec- 
tation value transforms nontrivially under the symme- 
try transformation. For a more detailed discussion on 
the interplay between the thermodynamic limit and the 
limit j — > 0, we refer to Appendix B 



We illustrate this discussion with examples. First, we 
consider classical Ising magnets on a lattice. The sym- 
metry transformation exerted on the Ising variables m, is 
a global reflection, m, — > -m, for all i. The Hamiltonian 
H[m] = -J Yji m i ,n i+i is reflection symmetric, meaning 
that 

H[m] = H[-m]. (63) 
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Since the functional measure J Yl t dm;<5(m?— 1) does not 
break this symmetry, we have for the effective action 



T[m] — T[— m]. 



(64) 



We call this a Z2-symmetry. 

Analogously, the microscopic action of cold atomic 



bosons given in Eq. (33 1 has a global U(l)-symmetry, 
meaning that it is invariant under the following global 
transformation of the fields, 



(*>*)' 



(65) 



with real parameter a. In the basis of real fields, <p - 
<pi + this corresponds to a rotation 



1 \ 




( 

cos a 


\ 

- sin a 




(66) 












sine 

V 


cos a 







in field space. Since the functional measure jDcp*Dip 
shares this symmetry, the effective action F[0*,0] pos- 
sesses the global U(l)-symmetry as well. 

By virtue of Noether's theorem, the global U(l)- 
symmetry in conjunction with a linearly appearing time 
derivative in the kinetic term of the microscopic ac- 
tion, leads to the conservation of total particle num- 
ber N = fJ<p*(Jt)<p(jt)). This is a characteristic fea- 
ture of nonrelativistic field theories. A brief review of 
Noether's theorem in the classical and quantum case is 



given in Appendix D 



The above mentioned properties of the effective ac- 
tion for vanishing external source have a profound con- 
sequence for the effective potential U. Indeed, from Eq. 
( [56] ) we deduce that the latter is not a function of (f>, <p* 
alone, but we rather have 



U = U(p), 



(67) 



where p is the most general combination of fields al- 
lowed by symmetry. For instance, we have 



P = 



|m 2 (Z2 - symmetry), 
1 0*0 (U(l) - symmetry). 



(68) 



We plot the effective potential U (p) for a second and 
first order phase transition in Figs. |7]and[8] respectively. 
The critical temperature T c (p) is defined such that the 
location of the minimum poipt, T) becomes zero - either 
continuously or discontinuously. In particular, for a sec- 
ond order phase transition we distinguish the following 
three cases: 

(i) po + 0, U'(po) = 0: phase with broken symmetry, 



T>T C 




T <T C 



Figure 7: The effective potential U(p) for vanishing external sources 
is a function of the symmetry invariant p. The latter is given by p = rrr 
or p = |0| 2 for magnets or ultracold bosons, respectively. Throughout 
a second order phase transition, the location of the minimum of the 
effective potential changes from po = to po > in a continuous 
manner. We have chosen here the temperature to be the control pa- 
rameter. However, since the effective potential depends on p, T, and 
the microscopic parameters of the theory (e.g. coupling constants), 
we may also drive the phase transition differently. 



T = T C 




Figure 8: In a first order phase transition, we have a jump in the order 
parameter po as we cross T = T c . From the plot of the effective po- 
tential we see how this discontinuous behavior can arise, although we 
smoothly vary the system parameters. Note that the effective poten- 
tial is actual a convex function, as the effective action originates from 
the Legendre transform of the convex Schwinger functional. The non- 
convex parts should therefore be replaced by straight lines according 
to the Maxwell construction, but this does not invalidate the overall 
picture of first order phase transitions. 



(ii) po = 0, U'(po) = 0: critical point, 

(iii) po = 0, U'(po) + 0: symmetric phase. 

In the broken phase, the location of the minimum 
Po of the effective potential does not necessarily com- 
pletely determine the thermodynamic equilibrium state 
of the theory. In the case of magnets, we have n& + 
and thus there is still the freedom to choose the sign 
of mo, which is a Z2-transformation. For the case of 
bosons, the condition \<Pq\ 2 + only fixes the amplitude 
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of the complex field <f>o = \(po\e l6 , whereas the phase 9 
can still be chosen arbitrarily. The possible nonequiva- 
lent choices are given by 8 e [0, 2n) U(l). We say that 
the condition on po singles out a manifold of possible 
ground states 0o> which in our examples is given by Z2 
and U(l), respectively. In the absence of explicit sym- 
metry breaking terms, the precise choice of the ground 
state in the degenerate manifold indeed happens spon- 
taneously - it is induced by fluctuations or perturba- 
tions due to the environment, which we can neither re- 
solve nor control ll53l . Nevertheless, this phenomenon 
is ubiquitously observed experimentally; for example, 
spontaneous phase symmetry breaking can be detected 
in interference experiments of initially disjunct conden- 
sates Il54ll . 

In Fig. [9| we plot the boson effective action in the tree 
level approximation T[tp] S [<f>] for a constant field in 
the complex 0-plane. The microscopic action S is given 



in Eq. ( 33 1. We write 



1 



T[<p = const.] = U(<f>) = -/#| 2 



(69) 



For obvious reasons, U((p) is also called the Mexican 
hat potential. Without loss of generality we assume the 
ground state tpo to be real, such that real and imaginary 
components of tp = d> l + \(p 2 are distinct directions in 
the complex plane. The ground state singles out the 
point (0o, 0). Now consider the field tp to be fluctuat- 
ing around this point. As usual, we write 



<p(r, x) = 00 + 8ip{T, x) 



(70) 



with (S(p) = 0. The fluctuations dtp are complex and 
can vary both amplitude and phase of (p. However, the 
fluctuations which increase the amplitude away from </>q 
have to climb up the hill and thus are energetically un- 
favorable, i.e. they are suppressed in the functional in- 
tegral by a term 



I 



(71) 



We call them radial or gapped excitations. In contrast, 
fluctuations of the phase are not hindered energetically, 
because they are along the well of the Mexican hat. 

The existence of a massless or gapless mode in a 
symmetry broken phase observed in the example above 
is a general phenomenon. In fact, it is an exact prop- 
erty of the full theory, as has been established by Gold- 
stone [55]. More precisely, Goldstone's theorem states 
that any spontaneous breaking of a continuous symme- 
try results in the appearance of gapless modes in the 




Figure 9: The Mexican hat potential from Eq. (69) only depends on 
the amplitude of the complex field <p = 0| + i(f>2- Thus, it reflects the 
U(l)-symmetry of the bosonic theory, which is invariant under phase 
rotations — * e m <f>. The ground state of the system will, however, 
spontaneously break this symmetry, e.g. by choosing 0o 6 R. 



excitation spectrum of the system. The proof of Gold- 
stone's theorem is very simple in the effective action 
framework. Since we are interested in a statement about 
the masses of the theory, i.e. properties of the system 
in the homogeneous limit of vanishing frequencies and 
momenta, we can restrict ourselves to the effective po- 
tential U. As we have seen above, U(p) only depends 
on the symmetry invariant p — \<p\ 2 . Consider the field 
equation = for the radial field. Since we assumed 
a homogeneous setting, this reduces to 

dU , ssb 

= — (to) = <hU (f>a) => U'ipo) = 0. (72) 
dtp 1 

Here, a prime denotes differentiation with respect to p. 
The mass of the Goldstone mode 6tp2 is then found to be 



d<f>j W 2 

(U'(p) + <fcU"(p))\ = U'(po) = 0. (73) 



* = U(p) ) 



We used the continuity of the symmetry by requiring 
p to depend smoothly on tf>i,2- We have carried out 
the proof for the symmetry U(l) =a 0(2). The above 
steps can be performed analogously for larger symmetry 
groups such as O(N), leading to N — I massless Gold- 
stone modes. In the above case, the vanishing of the 
mass term allows for strong fluctuations of the phase 
field in the phase of broken symmetry. In particular, 
they question the assumption of small fluctuations dtp in 
the functional integral, which lead to the one-loop for- 



mula given in Eq. ( 53 1. 
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3.3. Condensation of weakly interacting bosons 

We compute the properties of a weakly interacting Bose 
gas within the formalism introduced in the previous sec- 
tion. Mean field theory and the Gaussian approxima- 
tion, which is equivalent to the Bogoliubov theory in a 
second quantized formulation, already capture many ef- 
fects which are relevant for experiment. For reasons of 
stability, a Bose gas usually has to be dilute or weakly 
interacting such that these simple approximations work 
well. However, we show below that there are situations, 
where the Bose gas has to be treated by more sophisti- 
cated methods. 

We start from the microscopic action of the weakly 



interacting Bose gas given in Eq. ( 33 1, 



S [<p*, ip\ = J dr J J^p*(j, x) [d T - ^ - fij <p(j, x) 
+ V^nWiT, x)<p(j, x) + |(^*(t, $<p(T, f)) 2 } (74) 



with g > 0. We recall that <p is a complex field. As dis- 



cussed in Eq. (65 i, the microscopic action has a global 
U(l)-symmetry. 

The goal of this section is to calculate the correspond- 
ing effective action for this problem. Due to the inter- 
action term of fourth order in the field, the functional 
integral cannot be calculated in a straightforward man- 
ner and we rather have to rely on some approximative 
method. At sufficiently low temperatures (IjAj « 1), 
we can utilize the fact that the condensate is occupied 
by a macroscopic number of particles, <po ~ y/N. This 
scaling of the field expectation value had been found for 



noninteracting bosons in Sec. 2.3 for a sufficiently high 



spatial dimension. It is reasonable to assume that this 
scaling is also valid for the case of weak interactions. 
This identifies an ordering principle, which will justify 
our subsequent approximations. 

To understand the mechanism underlying an order- 
ing principle of this type, we study a simple toy model: 
Consider a smooth, real-valued function f(x) which has 
its minimum at Xq. Our goal is to evaluate the integral 



dxe- N ^\ 

DO 



(75) 



where N is large. Expanding the exponent around the 
minimum of /, we have f'(xo) = and /"(xq) > 0, 
leading to 



-Nf(x ) 



f 



Axe 



(76) 



However, performing a variable transformation x' — 
we find 



j _ e -Nf(xo) f°° dx ' c -kf"(x )(x'-x' o 



e - Nf<M yj27T/Nf"(x ), (77) 



where C is some constant. Here, means that the 
correction term vanishes in the limit N — > oo. If we 
are only interested in ^ log /, we can neglect the square 



root term in Eq. ( 77 1. Note that we did not require 



the minimum of f(x) at xq to be strongly expressed, 
because this is ensured by a sufficiently large N. This 
basic mechanism underlies all large N expansion strate- 
gies: after identification of a suitable parameter N, the 
remaining functional integral becomes Gaussian, just as 
in our one-dimensional toy model. 

In the context of weakly interacting bosons, we want 
to compute the effective action Y[<p* ,$]. We start from 



Eq. ( 50 1, which we write in condensed notation as 



-/ 



D6(pe 



(78) 



where <p(r, x) is a complex field vector. Expanding the 
exponent around <p, the linear terms cancel and we arrive 
at 



I 



D6<pe 



-si<p]- L 2 j XY s<p-s<- 2 n<f>]-s<p+ 



(79) 



Now, if the expansion point obeys = yfNcf)', where 
0' = 0(N°), we find S (2> [cf>] = NS' (2) [4>'], because the 
microscopic action only contains terms up to fourth or- 
der in the field. Accordingly, we have 



I 



D ^' c -iy^" 2 '-V+o(jv- 1/2 ) 
VA7 



(80) 



We will see below that we indeed have (f> ~ y/N for 
the expansion point. In comparison, the fluctuations 
scale with dtp ~ N° and their contribution is negligi- 
ble for large Af. The validity of the classical approxi- 
mation r ^ S , and subsequent improvements, thus re- 
lies on the existence of a macroscopically occupied con- 
densate. We emphasize that the term "classical" here 
refers to the absence of fluctuations, and not to the limit 
h — > 0. In fact, in the next paragraph we will discuss 
specific features which crucially build on a truly quan- 
tum mechanical feature: macroscopic phase coherence. 
In addition, we note that the approximation is not based 
on perturbation theory in the coupling constant g. The 
notion of weak interactions is, however, needed to jus- 
tify the scaling <f> ~ yfN derived from the noninteracting 
case. 
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Classical limit and Gross-Pitaevskii equation 



f/<PO 



From Eq. ( 52 1 we obtain an approximate expression for 
the effective action F, when evaluated for the solution 
of the classical equations of motion. It corresponds to a 
saddle-point approximation for the functional integral. 



With Eq. (51 1 we find the condition for vanishing exter- 
nal sources 

= s Or -a = ( dr ~ — ~ V + V ext (x)\ </>o(t, x) 
0<p*{T,x) \ 2M ) 



+ 8(<Po(t< %o(t, x)))<fo(T, x). 

(81) 

This equation can be analytically continued to real time 
t = -it and is then known as the Gross-Pitaevskii equa- 
tion [56, 57 1 . Restoring H, we find 



m—{t,x) 

at 



2vt2 



- H + V ext (x) <p (t, x) 



' 2M 

(ip* (t, x)<p (t, x))<po(t, X). 



(82) 



For vanishing coupling g, this formally is a Schrodinger 
equation for a single particle in an external potential. 
For this reason, the order parameter <po(t,x) is some- 
times called the macroscopic wave function and we can 
expect characteristic features from quantum mechanics 
to be found in weakly interacting Bose-Einstein con- 
densates. For instance, effects of phase coherence can 
be observed in a condensate, which as anticipated above 
is possible since the classical limit considered here de- 
rives from N — > oo and not h — > 0. For g + 0, the 
Gross-Pitaevskii equation is nonlinear and thus shows a 
richer spectrum of solutions than the Schrodinger equa- 
tion. 

One particular example for this interplay between 
nonlinearity and quantum mechanics is the existence of 
vortex solutions for the equations of motion with quan- 
tized phase. For this purpose, consider the situation of 
vanishing external potential, V ext (x) = 0. We may then 
look for a static solution with cylindrical symmetry ac- 
cording to (fo(t, x) = f(r)e lIe , where r — (x 2 + y 2 ) 1 ^ 2 
and is the polar angle. The winding number I must be 
an integer in order to guarantee the uniqueness of the 
macroscopic wave function as 9 wraps around the ori- 
gin. Plugging this ansatz into the equations of motion, 
we arrive at the ordinary differential equation 

0'-U r+f i-^)-» f+sf - ,83) 

The corresponding field configuration <po(jt) is called a 
vortex solution, with radial function / plotted in Fig. 




Figure 10: We locate the center of the vortex solution to Eq. J83J at 
position r = in the (x, y)-plane and plot the behavior of the ampli- 
tude fir). For distances which are small in comparison to the healing 
length = hl(2Mgpo) 1 ^ 2 , we have a strong deviation from the con- 
stant solution ipo = yjplg. For larger distance, however, the amplitude 
approaches this mean field value. We thus deduce the healing length 
to be the typical size of a vortex or, more generally, to be the char- 
acteristic length scale of possible inhomogeneities in the otherwise 
constant condensate amplitude. 



10 The solution has two qualitatively distinct regimes 
separated by the length scale related to the nonlinearity, 
the healing length 



= h/(2M gPo ) 



1/2 



(84) 



For length scales r » £h> the amplitude |^o(*)l ap- 
proaches the value ^nlg, which is a constant solution 
of Eq. (82 1. Instead, for r <k the condensate am- 



plitude must vanish due to the centrifugal barrier, which 
in turn roots in the quantization of the phase (the term 
~ l 2 /r 2 ). Quantized vortices are a hallmark of quan- 
tum condensation phenomena, and have been observed 
experimentally in bosonic |58| and fermionic |59| con- 
densed systems. The vortex core size £h gives direct 
information on the interactions in the many-body prob- 
lem. The action of a vortex field configuration is larger 
than that of a homogenous condensate. Therefore, it 
must be triggered externally, such as via rotating the 
trap, which imprints angular momentum onto the sys- 
tem. We refer to ll60l for an extensive discussion of 
these topics. 



Effective potential and spontaneous symmetry breaking 

The effective potential in the classical approximation 
F =! S is given by the Mexican hat potential discussed 
below Eq. 

Vext = 0, the solution to Eq. ( 8 1 
0o, which minimizes 



(69 1. Indeed, for the homogeneous case 
is given by a constant 



U{a» = -p\4>\ 2 + 



(85) 
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For p < 0, we have 0o = 0. For the cases of interest 
here, however, the chemical potential is positive and we 
have besides <*n = also the solution 



\(/>o(t, x)\ = po = pig G 



(86) 



which has a smaller value of the action S[(po] = F[0o]- 
We draw attention to the point that the amplitude -Jp/g 
was also found for the vortex solution far from the cen- 
ter of the vortex. The latter inhomogeneous solution 
necessarily has to approach this limiting value for the 
action to remain finite. 

Note that the phase of (po is not specified by the ac- 
tion principle 6S /dtp* - 0. Indeed, the global U(l)- 
symmetry of the microscopic action S is carried over 
to the equations of motion derived from it. Clearly, the 
actual ground state of the system must have a particu- 
lar value of the phase and without loss of generality we 
assume (po to be real, i.e. 



(87) 



As discussed above, this is a manifestation of sponta- 
neous symmetry breaking. 

From the effective potential U(<p), we can deduce the 



phase diagram and the equation of state using Eqs. (59 1 



and (60 1. We find in the classical approximation 



d U 

= ^(0o) = ("^ + #l0o| 2 )0o, (gap equation) 
o<p* 



(88) 

dU , 
n(p) = — — (0o) = m\ =Po, (equation of state). 

Ofi 

(89) 

The first equation yields the constant solutions from Eq. 
(86 1. From the equation of state n(p) = po we infer 
that all particles are condensed in the classical approxi- 
mation at zero temperature. In particular, the phase di- 
agram consist of the two regions, p > and p < 0, 
with and without particles, respectively. Furthermore, 
since the classical approximation we applied so far does 
not include any thermal fluctuations, we are restricted 
to zero temperature. To get a more physical picture of 
the phase structure and the thermodynamics, we have 



to improve formula ( 85 1 for the effective potential by 
including fluctuations. 



Quadratic fluctuations and excitation spectrum 

We now go one step beyond the classical limit and in- 
clude quadratic fluctuations in ip in the effective action. 



The treatment of quadratic (or Gaussian) fluctuations is 
often referred to as mean field theory, although some- 
times this term is also used for the classical approxima- 
tion r - S . The Gaussian corrections to the classical 
formulae for a weakly interacting Bose-Einstein con- 
densate have first been derived by Bogoliubov in the 
second quantized formalism [61 1. We are aiming at a 
discussion of this contribution with the help of the func- 
tional integral techniques developed so far. 

The correction to the effective action due to quadratic 
fluctuations is summarized in the one-loop formula ( |53] l, 

IW.fl =S[<f>*,</>] + ^TrlogS (2) [<f ,<^], (90) 

which is valid for small fluctuating fields. Here, <p(r,x) 
is an arbitrary complex field. After the trace has been 
evaluated, the ground state configuration O is found 
from the minimum of the full (one-loop) effective po- 
tential. However, to get a first idea of the underly- 
ing physics, we approximate the full ground state to be 
given by the classical solution tpo from Eqs. (81 1 and 



(86 1. In this case, Eq. (52 1 yields 



a -n<p<>] 



S -S[<p a ] J 



Dtp* Dipe 



^ jw^vk;,) 



(91) 



It is favorable to work in the real basis for the fluctuating 
field with transformation specified by (p = -q^ifi +1^2), 
and to switch to momentum space (using translation in- 
variance of the investigated situation). The quadratic 
fluctuations then take the form 



(<Pl-Q,tP2,-Q) 



S q + 2&O0 ~0)„ 



<Pl,Q 
<P2,Q 



(92) 

where e q — q 2 /2M, Q = {co„,q) and oj n is a bosonic 
Matsubara frequency, see Eq. ( [34] ). Moreover, we used 
Po = l^ol 2 = pig- 

Eq. (92 1 fully confirms our picture of fluctuations in 
the Mexican hat potential. Indeed, the real field 6ip\ con- 
stitutes the radial mode, which is gapped, i.e. it has a 
mass term 2gpo. This mass term suppresses the corre- 
sponding amplitude fluctuations of the field. The second 
kind of fluctuations, 6(p%, however, is not gapped and 
there can be arbitrary many excitations of this mode. It 
corresponds to phase fluctuations, which take place in 
the valley of the Mexican hat. It constitutes the Gold- 
stone mode associated to the spontaneous breaking of 
the U(l)-symmetry. 

Next we wish to calculate the excitation spectrum. 
To this end, we first note that the matrix appearing in 
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particles 




Figure 1 1 : The energy spectrum of Bogoliubov excitations consists of 
a phonon-like part for momenta q < q^, where q^ = h/^ is associated 
to the healing length t^. For larger momenta, the excitations show a 
quadratic momentum dependence, which is typical for nonrelativistic 
particles. 



Eq. ( |92| constitutes the full inverse propagator, or in- 
verse Green's function, G~ l {Q), within our simple ap- 
proximation (cf. the discussion of the effective action 
in Sec. |3.1| >. More precisely, it is defined according 
to r®(<2', QMo] = G l (Q)5(Q + Q') and thus coin- 
cides with the classical inverse propagator S (2) [^o] at 
this level of approximation. The excitation spectrum of 
the system, i.e. the dispersion relation, is obtained from 
the poles of the full propagator G(a>,q) after analytic 
continuation to real-time frequencies a>, which corre- 
spond to real times t. 

In order to find the right prescription for the ana- 
lytic continuation, we consider a field 4>{t,x) in eu- 
clidean time r. Its Fourier representation is given by 
0( T , jt) = J Q e i ^- ?+& '» T >0(g). The sign of the tempo- 
ral term in the exponent reflects the fact that time and 
space are treated equally in euclidean space. For real 
times f, we expect a reversed sign instead, describing 
wave propagation. Thus, using r = if, we find from 

co„ = ioj. (93) 

The poles of the full propagator at frequencies to - E q 
are then obtained via 

= det(G- V, = iE q , q)) = (e q + 2g Po )s q - E\. (94) 
This leads to the Bogoliubov excitation spectrum 

Eq = yjsqiSq + 2gf Q ). (95) 

The form of the spectrum of elementary excitations 
(95 1 has important physical consequences. For small 



momenta q' 
dispersion 



<s Mgpo, we have a linear and gapless 



c\q\, 



c = 



(96) 



which is characteristic for phonons with velocity of 
sound c, whereas the spectrum gets quadratic and thus 
particle-like for high momenta q 2 » Mgpo, 



11 
2M' 



(97) 



Thus, for long wavelength excitations (low momenta) 
the linear part of the spectrum dominates, whereas on 
short scales (high momenta) the particle nature of the 
system is still visible. The fact that the low-momentum 
degrees of freedom have a phonon-like dispersion hints 
at the typical collective behavior of many-body systems, 
where the effective quasiparticles do no longer coincide 
with the microscopic particle-like degrees of freedom. 

Below, we will see that the phenomenon of superflu- 
idity is intimately connected to this modification of the 
spectrum at low momenta. We can get more insight into 
the nature of these fluctuations by integrating out the 
massive modes tp\ in the functional integral Eq. <\9l\ . 
(This can be done by completing the square in the expo- 
nent). This procedure is useful if we focus on momenta 
with energies below the gap, e q <K 2gpo. It produces a 
"renormalized" low energy theory for the massless ex- 
citations if2, which reads 



rpo X 



4gpo 



■ c 2 q 2 )ip 2 , Q . 



(98) 



This low energy action reveals that the linear dispersion 
is due to the fluctuations of the phase. 

The characteristic momentum scale q^ which sepa- 
rates both regimes is given by the inverse healing length, 
qh = h/i;^, which we already encountered as the char- 
acteristic size of a vortex inside an otherwise homoge- 
neous condensate. The Bogoliubov dispersion and its 
limits are shown in Fig. 



11 



To complete the predictions available from Bogoli- 
ubov theory, we indicate the condensate depletion. 
Physically, this depletion is an observable effect of 
quantum fluctuations related to the presence of a finite 
interaction strength g, and results from scattering pro- 
cesses out of and back into the condensate. In partic- 
ular, it occurs in the absence of thermal fluctuations at 
T = 0. One finds 

1 f d?q ( Sq + gpq \ 
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V 



frame of reference of the 
moving object 



frame of reference of the 
ground state system 



Figure 12: Frames of reference for Landau's gedankenexperiment to 
identify superfluid flow. 



The calculation is rather lengthy, but instructive as we 
can learn from it how to deal with ultraviolet diver- 
gences in nonrelativistic quantum field theory. It is pre- 



sented in Appendix C 



Superfluidity and Landau criterion 



The excitation spectrum E q = sjs q (s q + 2gpo) allows 
for superfluidity. This phenomenon was first observed 
in liquid 4 He and manifests itself, for example, in the 
frictionless flow through small slits or capillaries. The 
onset of superfluid behavior constitutes a phase transi- 
tion, which is of second order for the weakly interacting 
Bose gas in three dimensions. The hydrodynamic de- 
scription of a superfluid system is altered by the sepa- 
ration of the macroscopic motion into normal fluid flow 
and superfluid flow, the latter being frictionless, irrota- 
tional and entropy conserving. The corresponding two- 
fluid hydrodynamic equations have been developed by 
Landau and coworkers (cf. [62 1). Above the critical 
temperature of superfluidity, which coincides with Bose 
condensation here, the two-fluid equations turn into the 
hydrodynamics of a one-component fluid. 

Condensation and superfluidity are related phenom- 
ena, but they do not necessarily coincide. Condensation 
refers to the macroscopic occupation of a single mode 
and thus is a statistical effect. In contrast, superfluidity 
refers to the response of a given system as we will see 
below. The fact that these are independent concepts be- 
comes particularly important in spatial dimension less 
than three, where the Mermin- Wagner theorem forbids 
the existence of a condensate. In contrast, superfluid- 
ity is present below a critical temperature. The related 
phase transition is known as the Berezinski-Kosterlitz- 
Thouless transition. For more information on this sub- 
ject, we refer to the literature (e.g. IBTl ). 



Landau established a beautiful criterion for the exis- 
tence of superfluidity, identified via the property of fric- 
tionless flow, from a simple yet general kinematic argu- 
ment. We consider an object uniformly moving in the 
liquid equilibrium state of a system at velocity v, and 
ask when it is favorable to create an excitation in that 
state, which leads to friction. To describe the situation, 
we introduce two frames of reference: S' is comoving 
with the object, 2 is the frame of reference for the liquid. 



(This is illustrated in Fig. 12 ) The general transforma- 
tion of energy and momentum under a Galilean boost 
with velocity v in these two frames is given by 



S: E, 



trP- 



S' : E' - E - p ■ v + \Mv 



p-Mv*, (99) 



where M is the total mass. We apply this transformation 
to the calculation of the total energy and momentum in 
the following two situations: 



(i) the ground state of the liquid 

I : £o, po = 5, 
I' : E' = E + \M#, f 



-Mv 



(100) 



(ii) the ground state of the liquid plus an excitation 
with momentum p and energy e p 



I : E ex - Eq + e t 



- r -p-v + \Mv 2 , 



(101) 



The creation of an excitation is unfavorable if E' ex -E' = 
e p - p ■ v > e p — \p\\v\ > 0. Thus, no excitations occur, 
and we have frictionless transport of the object, as long 
as the velocity is smaller than a critical velocity, 



\v\ < v c 



■ e P 
mm — . 

P \P\ 



(102) 



Moving the object with a larger velocity leads to fric- 
tion. Let us apply this criterion to two particularly im- 
portant systems. First, for a weakly interacting Bose 
gas, the low momentum linear dispersion is E p = c\p\. 
Thus, the system is superfluid with a critical velocity c. 
Second, in a free Bose gas, the scaling of energy with 
momentum is quadratic, E p = p 2 /2M. In consequence, 
a free Bose-Einstein condensate is not a superfluid; its 
critical velocity v c = 0. This is another example where 
the concepts of superfluidity and condensation do not 
coincide. 
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The Landau criterion can also be established from a 
field theoretical perspective. In order to study this prob- 
lem, we first need to analytically continue back to real 
times. The euclidean microscopic action S E and the 
real time action S are connected by the requirement that 
e~ SE = e lS appears in the path integral. Thus, starting 
from Eq. (74 1, we employ t = it to arrive at the real 
time microscopic action 

S[<p] = J dt jy(id, + ^+^- 8 -(<p*<p) 2 }. (103) 

Working in the Bogoliubov approximation, we expand 
this action to quadratic order in the fluctuations about 
the solution <po = ^fpQ = yfjlfg of the classical equa- 
tions of motion. Writing ip — <fo + A=(6(fii + id<p2), the 
corresponding real time action for the Bogoliubov exci- 
tations is found to be given by 



S bog [6ip u 5ip 2 ] = ^ J dt JgPo 
+ ~ J dt J(6(p u 6<f 



d t 



2M ) 



6tf! 

8<p 2 t 
(104) 



In view of identifying a critical velocity, we ask 
whether this action describing the excitations on top of 
the condensed ground state is stable under the transfor- 
mation of the field 



(p(t, x) — > <p'(t, x) = e 



Hp-x-Et) 



<p(t, X). 



(105) 



This describes the imprint of a plane wave with momen- 
tum and velocity p = Mv on the fields, and we com- 
ment on the temporal phase rotations below. In order 
to understand this transformation, we first decompose 
(pit, x) - (po + 6ip(t, x), where <po describes the condensed 
ground state, and S(p(t, x) small fluctuations around it. 
In particular, we see that under the transformation, the 
ground state picks up a position dependent phase 0o — > 
tp'Jt,x) = <p$e^ p ' x ~ Ef) and thus carries a supercurrent 
/= ^[V0 o *(f,f)^(f,f)-^(f,f)V^(f,f)] = <p* Q cf> 1 = 
<pQ<p{)V. We will now ask whether this supercurrent (su- 
perfluid flow) is persistent or stable, and thus the system 
a true superfluid ^To this end, we apply the transforma- 
tion ( |105| l to the microscopic action in Eq. ( 103 1, and 



"We note that this is not a Galilei transformation describing 
the change of frame of reference, given by f(t,f) — » if'(t,ic) = 
e^P '<p(t, 5t— vi) (the full microscopic action must be invariant un- 
der such transformations). It rather is a local gauge transformation. 



find 



s[y] -»sm + 



p 

2M 



y (io6) 



The temporal component of the transformation de- 
scribes an adjustment of the zero of energy, and, choos- 
ing E = p 2 /2M, we shift the latter back to its original 
value. 

In the momentum representation of the fields, 
<p(co, q) — J e^'^'^tpit, x), we find the quadratic part 



of the transformed Bogoliubov action in Eq. ( 104 1 to be 
given by 



+ 2gPo -i(w + v-q) 

i((L> + v-q) E CJ 



(107) 



The system is stable under the transformation in Eq. 



( 105 1, and therefore supports superfluid flow, if the ex- 



citation energies resulting from the matrix in Eq. ( 107 1 
are positive - in this case, the system is located at a (lo- 
cal) minimum in energy. Diagonalization shows that the 
effect of the above transformation is to shift the Bogoli- 
ubov excitation energies according to 



E' 



q>E q 



> 0. 



(108) 



102 1 for the 



We thus recover the Landau criterion Eq. 
critical velocity: Below the critical velocity for the per- 
turbation, the excitation energies E' for all modes are 
positive, and thus the equilibrium state (p'Jt, x) carrying 
superfluid flow is stable. Above v c , however, there exist 
unstable fluctuations which ultimately tend to destroy 
the superfluid flow, driving the system to a state which 
no longer is described by (p' (t, x). We note that the con- 
siderations also hold for the noninteracting case, where 
E q = q 2 /2M, resulting in a vanishing superfluid veloc- 
ity. 



Validity of Bogoliubov theory 

To close this section, we discuss the validity of Bogoli- 
ubov theory. In particular, our analysis will reveal why 
many experimental observations on cold trapped bosons 
are captured within this framework. 

We have seen in Eq. ( 80 » that the validity of Bo- 



goliubov theory is related to the ordering principle 
of a macroscopically occupied condensate, which al- 
lows for an approximate evaluation of the path inte- 
gral. Obviously, such a procedure breaks down if 
no condensate exists. This situation is found in two- 
dimensional systems at nonzero temperature and always 
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in one-dimensional systems. (See our discussion of 
the Mermin- Wagner theorem and its relevance for cold 
atoms in Sec. |2.3| ) In these lower-dimensional set- 
tings, one necessarily has to rely on nonperturbative ap- 
proaches. 

Even in three dimensions it is questionable whether 
an expansion in powers of the fluctuating field dtp is 
valid for low momenta. Indeed, from Eq. ( |92| for the 
classical inverse propagator we find for the occupation 
of the q-mode 



J CO 



(SiPq6(Pq) 



1 



1 

w 



(109) 



The high, diverging occupation of low momentum 
modes, allowing to roughly count 6<pg ~ \q\~^ 2 , ques- 
tions the validity of a simple ordering principle set by 
the macroscopic condensate occupation 0o ~ yN. 

In view of estimating the momentum scale where Bo- 
goliubov theory breaks down, we study the perturba- 
tive effects on the self-energy X for weakly interacting 
bosons at zero temperature (63). The full inverse prop- 
agator in the (<f>, 0*)-basis is given by 

G~\P) 

{ ia) n + s p -p + E n (P) Z an (P) 

(HO) 

We may regard Bogoliubov theory as the tree-level self- 
energies 



Jf (P) = 2gp , 



: 8Po- 



(111) 



The leading perturbative corrections are shown dia- 



grammatically in Fig. 13 The second diagram has an 



infrared divergence, which is logarithmic in d — 3 spa- 
tial dimensions and polynomial for d < 3. Indeed, the 
low momentum contribution to the corresponding loop 
integrals is given by 

EiV) ~ ?£t(P) ~ ~g 2 PQ f G 22 (Q)G 22 (P+Q) (112) 
Jq 



with 



G 22 (P) 



2gPo 



(l) 2 + c 2 ^ 1 



(113) 



and c = (gpo/M) 1 ^ 2 . For very low momenta we find 
^\P^0)~g 2 p M 2 c [ A d d+l q~ 4 

J,, <r 



g PoMp h . 

\p * 



log(A/p) (d = 3) 



(d<3) 



(114) 




P + Q 



t Q<fo 



Figure 13: First order perturbative corrections to the self-energy. 
Whereas the diagram to the left is infrared regular, the one to the right 
is divergent for low momenta P, see Eq. jl 14). 



where q = (a>/c,q) is a (d + l)-dimensional vector and 
q = \q\. In this regime we can use A = with p^ = 
(Mgpo) 1 ^ 2 - Mc being the momentum scale associated 
with the healing length, discussed earlier in this section. 

We expect perturbation theory to break down when 
both zeroth and first order corrections to the self-energy 
become of comparable size, i.e. 

|2<V np )| * \L m (P np )\. (115) 

From this relation we can deduce a characteristic 
momentum scale p np where the superfluid becomes 
strongly correlated and has to be described nonpertur- 
batively. We arrive at 



Pnp 

Ph 



exp(-l/gMp h ) (d = 3) 
gM (d = 2) 

(gM/ Ph yt 2 (d = i) 



(116) 



In three dimension we have g oc a. Together with po ~ n 
we find gph oc (a 3 «)^ 2 . Thus, in the dilute regime of 
small gas parameter, the nonperturbative physics hap- 
pen at exponentially small momenta. In two dimensions 
instead, where the coupling constant is dimensionless, 
the condition p np as p n is reached for g being of order 
unity. 

The finding of Eq. \\ \6) can also be expressed as 



Pnp 
Ph 



exp(-l/| 3 / 2 ) 



(4 = 3) 
(d<3) 



(117) 



The dimensionless quantity g constitutes the ratio of in- 
teraction versus kinetic energy in the nonrelativstic su- 
perfluid, 



gPO 



E m l/(M£ 2 ) 



gM P \; 2ld ~( Ph {) 2 . (lis) 



We used here again po ~ n for the weakly interacting 
condensate and I = rT l l d is the interparticle spacing. 
Accordingly, superfluids can be classified as 

(1) weakly correlated: We have g <g 1 an d the refore 
from Eqs. 



117i and (118l. Bo 



Pnp « Ph « 
goliubov theory is valid for a large part of the spec- 
trum, namely all modes with momenta \q\ > p D „. 
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Figure 14: The hierarchy of length scales is cut oft' by the trapping 
potential with characteristic size f osc . Fluctuations with larger wave- 
length do not appear in the system due to the finite extent. Since 
perturbation theory breaks down beyond that scale only, nonpertur- 
bative effects are difficult to observe in experiments and Bogoliubov 
theory is a sufficient approximation. The healing length fh. which is 
the characteristic size of a vortex, is smaller than £ osc . The observa- 
tions of vortices and vortex lattices in trapped systems confirm this 
picture I58II59I . 



(2) strongly correlated: For g ^ 1 we find p np — Ph — 
t~ x . Bogoliubov theory breaks down. 

In typical traps, the oscillator length ^ osc is smaller than 
the scale £ np = l/p np - Therefore, the external potential 
provides an infrared cutoff towards the strongly corre- 
lated regime. In the spatial continuum discussed here, 
the effects of fluctuations with very large wavelengths 
are not encountered in a finite trapping geometry. The 
observation of vortices, which are of size £h separating 
the particle- from the phonon-like regime, shows that 
£~ <K £h ^ ^osc- We note, however, that the strongly 
correlated regime can be reached by reducing the ki- 
netic energy with respect to the potential energy, thus 
decreasing g . This happens for bosons on a lattice close 
to the phase transition between a Mott insulator and a 
superfluid [64 65 1. The scale hierarchy for bosons in a 
harmonic trap is summarized in Fig. 



14 



We conclude that Bogoliubov theory provides a good 
description for many experimental settings. It may 
break down for special geometries like optical lattices 
and lower-dimensional traps, which then allow to study 
nonperturbative effects beyond the mean field approxi- 
mation in experiment. 

What happens beyond the scale £ np , where perturba- 
tion theory is plagued with infrared divergences? It has 
been recognized a long time ago that a phase-amplitude 
(hydrodynamic) description, useful for the regime of 
distances larger than £ np , does not have these problems 
ll66l . Based on this, an exact argument properly taking 
into account the coupling of phase and amplitude mode 
has been given that there must be two massless modes 
Il67l - l69l . This is in stark contrast to the naive expecta- 
tion from Bogoliubov theory, which predicts one mass- 
less (transverse) and one massive (longitudinal) mode 
as we have seen above. In a renormalization group lan- 
guage, this behavior is reflected in a fixed point of the 
scale dependent interaction coupling gk at zero, such 
that the longitudinal mass ~ gk^op vanishes ll70H72l . 
In addition, a symmetry enhancement from Galilean to 



Lorentz symmetry has been observed. This is reflected 
by the coefficient of the linear frequency term co tak- 
ing fixed point at zero, while the coefficient of the term 
u> 2 arrives at a finite value (72). This is possible since 
Galilean invariance is broken in a many-body system 
even at zero temperature due to the presence of a con- 
densate. The full crossover from the particle-like to the 
phonon-like region at momentum scale l/£h, and then 
from the phonon-like to the hydrodynamic regime at 
momentum scale 1 /£ np , has been followed continuously 
in unified approaches based on FRG techniques 1631172] - 

E7J. 

The vanishing of the longitudinal mass (or the diver- 
gence of the longitudinal static susceptibility x^M> — 4 — 
0) ~ l/gp) also has interesting physical consequences. 
For example, it has been shown that at zero tempera- 
ture and in spatial dimension d — 2 the (real frequency 
domain) spectral function possesses a critical contin- 
uum which starts directly above the spin wave pole at 
o = clgl 1551175117811791. Physically, the large spec- 
tral weight stems from the possibility of the (naively 
gapped) amplitude mode breaking into a pair of spin 
waves for any a> > c\q\. 

3.4. Superfluidity of weakly attractive fermions 
In this section we discuss the Bardeen-Cooper- 
Schrieffer (BCS) theory of weakly interacting two- 
component fermions applied in a cold atoms context. 
After a brief survey on the main statements of BCS 
theory we solve the theory in a Gaussian approxima- 
tion. On the technical side we perform a Hubbard- 
Stratonovich transformation in order to introduce effec- 
tive bosonic degrees of freedom into the fermionic the- 
ory. While the boson field here plays the role of an aux- 
iliary degree of freedom describing a Cooper pair con- 
densate, it will acquire a more direct meaning when con- 
sidering the BCS-BEC crossover in an FRG framework, 
where it will play the role of a molecular bound state in 
the BEC regime of this problem. Many of the formulae 
we derive in this section will be useful for the analysis 
of the strongly interacting case. The bosonic degrees of 
freedom also play an important role in view of structur- 
ing the phase diagram, and we will be able to apply our 
knowledge on phase transitions for Bose systems here. 
We will, in part, also adopt a "purely fermionic point of 
view", and show that the BCS instability is related to a 
divergence of the fermionic four-point vertex. 

BCS theory has been developed for superconduc- 
tors and thus is originally a theory of electrons, which 
are two-component fermions due to the spin 1/2. In 
ultracold quantum gas experiments, two-component 
fermions are realized by two distinct hyperfine states, 
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for example of either 6 Li or 40 K. A balanced mixture of 
the spin components is not fundamental in this context, 
but can easily be achieved via proper spin polarization. 

We describe the system by an effective Hamiltonian 
with local interactions. As in the case of bosons, for 
fermionic isotopes of alkali atoms, it is given by 



(119) 



6 = j d 3 xltf(H)(~ 

The creation and annihilation operators satisfy anti- 
commutation relations. The corresponding action reads 

S W, if/] = JT* dr j d 3 x L\t, x)(d r - ^ - nWr, *) 



(120) 



with independent Grassmann fields if/ = (if/i,if/2) r and 
if/* - (if/*, if/Vf- (A shift in p. from normal ordering, a 
prerequisite for the construction of the functional inte- 
gral, has been absorbed.) We formally defined the oper- 
ation 1 = (*)'. The microscopic coupling constant A will 
receive contributions from fluctuations. Including them 
leads to the renormalized coupling /Ir, which is then re- 
lated to the scattering length measured in experiments 
by relation (15 1, i.e. A R = 4-nh 2 a/M. 

We begin with some qualitative remarks. The expec- 
tation value of the fermion field vanishes due to Pauli's 
principle, (if/j) = for i = 1,2. However, there is no 
fundamental principle preventing a nonzero pairing cor- 
relation 

(if/^2) * 0. (121) 

As we will see below, this behavior is indeed found for 
low temperatures and attractive interactions. A descrip- 
tion of this phenomenon within the BCS framework is 
possible for 



a<0, \k ¥ a\ ^ \a/€\ <K 1. 



(122) 



It is a key feature of BCS theory that the small interac- 
tion scale cannot substantially modify the Fermi sphere 
which we encountered in the discussion of noninteract- 
ing fermions in Sec. 2.3 In particular, p = sp(n) re- 



121 1 is equivalent to 



mains valid. 

The nonvanishing correlation 
a condensation of bosonic quasiparticles in their low- 
est energy state. Indeed, suppose two fermions build a 
composite bosonic state. The relative energy of the part- 
ners will be minimized for a spin singlet and vanishing 
center of mass energy. Thus, the total momentum of the 



boson is zero and the momenta of the partners are oppo- 
site. Due to the fermionic origin of the excitations, the 
relative momenta lie on antipodal points of the Fermi 
surface. We arrive at a pairing 



(if/i(e F ,p)if/ 2 (e F ,-p)} + 0, 



(123) 



which is local in momentum space. This is a very small 
effect as we will quantify below, since it only occurs in 
the vicinity of the Fermi surface. The composite bosons 
just described are called Cooper pairs. Note that Eq. 



( 123 1 results in a ground state of the many -body system 
which breaks the global U(l)-symmetry of the action 
TM- 



Pairing field and Hubbard-Stratonovich 
transformation 



When looking at the action in Eq. ( |120| i, we may won- 
der whether the quartic Grassmann field term is iden- 
tically zero. Indeed, this would be the case for a one- 
component Fermi gas (which thus does not have local 
interactions). However, for our two-component spinors 
we have 



oaV) 2 = 



if/ 2 



2if/\if/* 2 ip x ifj 2 . 



(124) 



With the fully anti-symmetric tensor in two dimensions 
£;./ = Sji, £12 = 1, i-e. 



E 



1 
-1 



(125) 



we easily find 



(<aV) 2 



1 



0AV*)(iA'«A). 



(126) 



This rewriting is called a Fierz transformation. Note that 
(if/'eif/)* = -ifAsiff. 



From Eq. (B. 25 1 we deduce the simple identity 



J Dip* Dip exp j- J m 2 (ip* - tp. d )(ip - ^ b ) j = N. 

(127) 

Herein, ip a (X) and ipb(X) can be chosen arbitrarily, be- 
cause they can be eliminated by a shift of the corre- 
sponding integration measure. Note however that m 2 
has to be positive for the integral to converge. For this 
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reason, we wrote it in a suggestive manner as a square. 
N is related to the determinant of the identity operator, 
but irrelevant for our purposes since we later take the 
logarithm of the corresponding expression. The idea 
behind the Hubbard-Stratonovich transformation is to 
insert this unity into the path integral J Dfte ~ s and then 
choose the free parameters in such a way that the ac- 
tion gets more suitable for subsequent approximations. 
In our case, we want to eliminate the quartic fermionic 
interaction term ~ A(ft'ft) 2 ; it is traded for the complex 
field ip as additional bosonic degree of freedom. As it 
turns out, this is actually a very physical effect. 
We choose 

V* = fb = ^(ft'eft), (128) 

2m z lm L 

where the physical meaning of the constants m 2 and h 
will become clear below. This leads to 

N = J~D^*Dy> 
x exp j- J m 2 (ip* - 
= ^DyDtp 

x exp j- J m 2 (ip* + ^(i^eifr*))(<p - ^( | A' el A))| 



JL(ft>sft)*)(<p~^ 1[ (ft>sft)) 
lm L lm A 



D^*Dy>exp^ - | ynqftp 



+ h - [ipft^eft* - <PV«^) - ^(ft' t sft*)(ft t sft) 



(129) 



Obviously, inserting this into the coherent state path in- 
tegral for the fermions 



Dft*Dfte 



-Shir*,*] 



(130) 



is equivalent to having a theory with both fermions and 
bosons and action 

S [ft*, ft, ip*, ip] = £ Ll(d T - ^ - fx)ft 



■ m 2 ip*ip + — (tpif/'sft* — (fi'iffsij/) 



(131) 



Before proceeding, we remark on the strategy of 
the outlined procedure. The rewriting of the purely 
fermionic theory in terms of a theory of both fermions 



and bosons is exact and did not involve any approxi- 
mation. On the other hand, we obviously did not gain 
anything so far, because the functional integral still has 
to be evaluated. However, we need to recognize that we 
will essentially never be able to perform the full func- 
tional integral. Therefore, it is reasonable to reformu- 
late the theory in such a way that already the leading 
order captures the physically most relevant phenomena. 
Based on physical insight, this is achieved by intro- 
ducing the proper collective degrees of freedom. We 
found above that for attractive fermions, condensation 
of pairs should be the relevant mechanism. Our choice 
(pi — (fib ~ ft' eiff in Eq. ( 128 1 effectively substitutes tp for 



ft' eft — 2ft\ft2 in the action. Thus, tp is directly related 
to the pairing order parameter {ft\ft2>- A condensation 
of <p will then yield a nonvanishing pairing correlation. 
Therefore, we can hope that already the introduction of 
tp and a Gaussian approximation to the path integral can 
be sufficient to describe the theory. 

We emphasize, however, that the choice of the 
bosonic degree of freedom, is not unique from a math- 
ematical point of view. In fact, the "wrong" choice of 
tp a and (pt> does not lead to a satisfying result. For ex- 
ample, if we had chosen cp a and <p\, as hermitean bi- 
linears, such that <p ~ ft\fti, then this would describe 
well the features of a theory with particle-hole-pairing, 
but fail here poorly, because the instability occurs in the 
particle-particle-channel. 

For the case of attractive interactions, A < 0, the ac- 



tion in Eq. ( 131 1 can be simplified by choosing 



h 1 



(132) 



Thus, only the ratio h 2 /m 2 is a physical quantity. In 
particular, the sign of h is irrelevant. This can also be 



inferred from Eq. ( 128 1, because any phase of /; can be 
absorbed into a redefinition of tp K and <p\,. By a rescaling 
of the bosonic field according to <p — > hip, we arrive at 
the microscopic action 

S [ft*, ft, ip\ ip] = ^ L^(d T - ^ - p)ft 

+ + \ (<pft f sft* - ^ft'sft)^ , (133) 
which only depends on the single parameter A. 



One-loop effective potential 

We perform the calculation of the one-loop effective 
potential analogously to the case of weakly interacting 
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bosons in Appendix C Starting from the action in Eq. 
( |133| i, we introduce Nambu spinors 



^2 ) 



(134) 



The action can then be expressed as 



<p] =-- \ <p*<p+ \ 

* JX JX 



p <p 

^ P' j 



(135) 



which is manifestly quadratic in the fermions, with 

V 2 V 2 

P = d T u, P'=d T + +u. (136) 

2M H 2M P 



Using the generalization of Eq. (48 i for the effective ac 



tion r[ l t" , V F, <p*, 0], we can now evaluate the Gaussian 
integral over the fermionic fields and "F. This results 
in a theory F[0*, of interacting bosons. 

The BCS approximation consists of integrating out 
the fermions and neglecting the boson field fluctua- 
tions. Therefore, it is a mean field approximation for 
the bosonic degrees of freedom. We choose a constant 
background field 0, which is not yet evaluated at its 
equilibrium value. The second functional derivative of 



the action is found from Eq. ( 135 i to be given by 



5{Q + Q') 



UO n + S q -fi 



(137) 



The matrix appearing in this expression is the inverse 
classical propagator G V \ V (Q). For the effective poten- 
tial U - T//3V we obtain 

U(p = fc/>) = -jfcp - iTrlogS (2) [0*,0] 

= ~JP- JlogdetG^CQ) 

= -jp- f lo s( w « + ( £ ? _ V? + p) 

= - l -p-2T f^Vogcosh(^). 
X J (2nf B V 2T ' 



(138) 



We used 2„log(l + ^r^jj^p) - 2 log cosh x and 
dropped an overall constant in the last line, which is ir- 
relevant for the thermodynamics. We also introduced 



E q (p) = ^(e.-p^ + A 2 



(139) 



gap A 



1/2 



Figure 15: We plot the function E q = J(e q — flj 1 + A 2 from Eq. 
\139) for a nonzero gap A. For A = 0, the branch of the disper- 
sion relation hits zero at q = yffi and, therefore, it is even at very low 
energies possible to excite particles with momenta at the Fermi sur- 
face. For A > 0, these excitations are costly and at very low energies 
or temperatures, we cannot excite them. In the BCS case, the gap is 
exponentially small. 



with A 2 = p. 



Excitation spectrum 



The order parameter po(p, T) 
gap equation 

dU 



i is found from the 



= ^(0o) = 0o • U'(po). 



(140) 



The three types of solutions to this equation have been 



discussed in Sec. 3.2 Given a nonvanishing order pa- 
rameter po > in the phase with spontaneously broken 
symmetry, the fermion excitation spectrum found from 



detGy]y((o n = iE q ,q) = acquires a gap, 



+ A 2 . 



(141) 



Due to the appearance of the gap, even the excitation 
of fermions with momenta at the Fermi surface is sup- 
pressed. (Note that the fermion dispersion vanishes for 
s = e f in the symmetric regime.) The absence of single 
particle fermion excitations at sufficiently low energies 
has an important implication: The lowest-lying modes 
are bosonic phonons with linear dispersion. This results 
in superfluidity of the system according to Landau's cri- 
terion discussed in Sec. 13.31 



Ultraviolet renormalization 

The onset of superfluidity is determined by the condi- 
tions po = and U'ipo) = 0, which correspond to the 
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critical point. Using Eq. ( 138 i for the effective poten- 
tial, we conclude from the gap equation that this partic- 
ular point is given by 

f d 3 g 8E q 



1 

A 



(2nf dp 



>p= 



tanh(M^)) 
o V 27V ' 



^tanh(: 



27V 



-)■ (142) 



The integral is linearly divergent as the integrand tends 
to unity for q — > oo. This UV divergence is due to our 
simplifying model assumption that the interactions are 
pointlike and thus constant in momentum space for ar- 
bitrarily large momenta. In reality, the coupling A(q) 
is cut off smoothly at large momenta. However, as we 
have seen in Sec. 



2.1 



the details of the interatomic po- 
tential are not essential for low energy collisions, which 
are solely determined by the scattering length. We visu- 



alize this situation in Fig. 16 



We cure the divergence in Eq. ( 142 1 by a proper UV 
renormalization, where the goal is to trade the "bare" 
coupling constant A against a physically observable one, 
in this way eliminating the UV divergence. The proce- 
dure consists of two steps: (i) We regularize the mo- 
mentum integral by the introduction of a sharp UV 
(high momentum) cutoff A. A is then interpreted as a 
bare coupling, which depends explicitly on the cutoff 
A. (ii) Next we perform the renormalization. To this 
end, we observe that in experiments, we actually mea- 
sure a renormalized coupling /Ir at low energies, which 
necessarily includes the effects of quantum fluctuations. 
Therefore, the bare coupling is not accessible to us and 
might have a very large or very small value. Express- 
ing everything in terms of the renormalized coupling A R , 
the cutoff A will eventually drop out of the results. To 
compute /Ir, we consider the vacuum situation where 
T = p, = po = and the excitation of bosons is sup- 
pressed. The renormalized boson mass U'(po) - m 2 
should then be positive and equal to -1/Ar, see Eq. 
\\32\ (with h absorbed into the fields). In the vacuum 
limit, we deduce from the explicit form of the effective 
potential in Eq. ( 138 i the relation 



1 



1 

An 2 



f 

Jo 



dq 



(143) 



q 2 _ 1 MA 
s n A 2tx 2 

142b, we can send 



Inserting this expression into Eq 
A — > oo. 

This yields the renormalized gap equation which de 
termines the critical temperature T c (p, /Ir), 



1 2M r°° 

= r dq 

A R An 2 Jo ^ 



E q -p. 



tanh^- 



27V 



H 



(144) 



pointlike interaction 




"\ 




\ i sharp cutoff 
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■ 




i 
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Figure 16: The effective Hamiltonian in Eq. \\ \9\ assumes a pointlike 
interaction, which is valid for the scales kf, a\Aj l , we encounter 
in ultracold atom experiments. However, when calculating the effec- 
tive action in perturbation theory, we are confronted with divergences 
when integrating over all momenta. These singularities arise because 
of momenta if > A 2 . At these scales, the microscopic details of the 
interatomic potential can be resolved and we cannot rely on a pointlike 
approximation, see Fig. [2] We cure the problem by observing that the 
true coupling A(q) is derived from a more realistic potential and falls 
off smoothly in the UV. This cannot be described in the pointlike ap- 
proximation, but is taken into account by introducing a sharp cutoff at 
A. 



We note a difference to the usual condensed matter argu- 
mentation for BCS theory in solids. There, the validity 
of the approximation of pointlike attraction is restricted 
to energies in the vicinity of the Fermi surface, where 
the attraction is mediated by phonons. Therefore, the 
cutoff (Debye frequency) is much closer to the Fermi 
surface, cf. e.g. IBTl . Nevertheless, observables such as 
the critical temperature display the same functional de- 
pendencies, as they are dominated by effects very close 
to the Fermi surface as we will see next. 



Critical temperature and paired state 

We already mentioned that the weak interactions do not 
significantly modify the Fermi surface and we can as- 
sume p — ep. This will also be justified below by a 
direct computation of the equation of state. Rescaling 
the momenta in Eq. ( |144| > with the Fermi momentum, 
q = q/kp, and using sp = TV = k 2 /2M, we arrive at the 
condition determining the critical temperature T c 



1 

kpa 



2 r 

x Jo 



dq 



^tanh(^ll)-l). 
-1 \2TJTp) j 



(145) 

We introduced the scattering length for distinguishable 
fermions according to Ar = 4na/M < 0. For T c /Tp — > 
0, the integral develops a logarithmic singularity at the 
Fermi surface, where q — 1 . Thus, we obtain a solution 
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of Eq. ( 145 i for arbitrarily weak interactions a. The 
corresponding BCS critical temperature is 



To 



»r-2 



- e 2hh 



(146) 



Here, y = 0.577 is Euler's constant such that the overall 
prefactor becomes 0.61. Note that T c /T$ is exponen- 
tially small. 

For T < T c , a gap po > develops, which cures the 
logarithmic divergence. In particular, we find from the 
gap equation at zero temperature 



p (T = 0) 



n 
e~r 



0.57. 



(147) 



In the spirit of Landau's theory of second order 
phase transitions, we expand the effective potential (i.e. 
the grand canonical potential) around its minimum po. 
Since the effective potential depends on the renormal- 
ized quantities including all fluctuations, we have 

U(p,p, T) = -P(ji, T) + m 2 R (p - po) + ^(p - po) 2 . 

(148) 

In the symmetric phase, we have po = and, accord- 
ingly, rrPp R = > 0. The order parameter po acquires 
a nonzero value when either R becomes zero or Ar 
diverges. The fact that the "mass" m? R of the bosons 
becomes zero when going from the symmetric to the 
broken phase will later be important in our renormaliza- 
tion group study. In particular, we will employ a trun- 
cation of the effective potential which closely resembles 
the shape of Eq. ( |148[ ). 

Let us finally comment on the relation to a "purely 
fermionic" approach which does not introduce bosonic 



degrees of freedom. To this end, we compare Eqs. ( 143 i 



and ( |144[ ). We have interpreted the left hand side of 
Eq. ( |143[ ) as the inverse vacuum interaction parameter, 
which includes the effects of vacuum quantum fluctu- 
ations - in an effective vacuum action, it would appear 
as the renormalized four-fermion (two-body) interaction 
term ^(^V) 2 - I n full analogy, we can interpret the left 



hand side of Eq. ( 14-4 1 as the inverse "many-body" in- 
teraction parameter l/A m \,(T,fi), which in addition to 
the vacuum quantum fluctuations also includes many- 
body fluctuation effects related to the many-body scales 
r 1 = k F ~ y/Ji,A T ~ y/f. (Note that A mh (T = 0,p = 
0) = /1r.) The zero of the inverse many-body interac- 
tion strength, i.e. the divergence of the many-body in- 
teraction vertex, signals the onset of new physics, more 
precisely an instability towards a state which needs a 
qualitative modification of the theoretical approach for 



its proper description (such as a nontrivial minimum in 
its effective potential). If this phase is understood qual- 
itatively, and supposed to be described quantitatively, 
the approach via Hubbard-Stratonovich decoupling is 
preferable. On the other hand, the purely fermionic ap- 
proach is less biased and may serve as a guide to iden- 
tify the proper decoupling channel. We emphasize that 
the possibility of reducing the question of stability to a 
single parameter is remarkable in a many-body context: 
The full four-fermion vertex of our problem is, in prin- 
ciple, a tensor Aq u ^q 4 depending on four 4-momenta. 
This highly nontrivial reduction of complexity in the 
vicinity of the Fermi surface has been shown by Shankar 
|, and is reviewed in [81 1. 



Equation of state 

Besides the phase structure, we also obtain the ther- 
modynamics from the one-loop effective potential. Eq. 



( 138 1 implies for the density 



—< 
dp 



C d J o ,E Q , 
n(p.T) = (/.„) - -2 J ^jtanh( ; q 



^2T> 



= 2[[^L^—- f (149) 
U (27T)W + l J {2nf2}- y ' 

In the second line, we have split up the integral into 
a physical contribution and an artefact from the func- 
tional integral. It accounts for the relation between 
field expectation values and operator expectation values 
(6i//*6if/) = (fl'a) - i for each momentum mode. (We 



refer to Appendix C and Eq. (C.18 i for a more detailed 
discussion.) We thus arrive at 



n(p,T) 



= 2 fix 



1 



•)3 e E q IT + ! ' 



(150) 



From Eq. ( 150 1 we find our above considerations to 
be consistent: The Fermi surface is clearly expressed for 
small temperatures and an at best exponentially small 
gap or order parameter identified above, cf. Eqs. ( 146 1 



and ( 147 i, and we have p ep. Despite its smallness, 
the nonzero gap of the dispersion relation around the 
Fermi energy results in an important qualitative effect, 
namely the superfluidity of the system. 



Validity and experimental (ir)relevance of BCS 
superfluidity 

We close this section by an estimate of the temperature 
regime which has to be reached in order to observe BCS 
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superfluidity in cold atom experiments. We found in Eq. 
(146 1 that for a given density (i.e. Fermi momentum 



kf(n)) we have 



r F (n) 



0.61e 



-jr/2\kpa\ 



(151) 



Due to \k F a\ « 1, we have an exponential suppression 
of the critical temperature. To experimentally reach this 
exponentially low temperature regime is difficult. 

We compare to a bosonic system. For particle number 
density n and boson density n B , we can formally define 
a Fermi temperature via = (67T 2 n/gjs) 2 ^ /2M. To- 
gether with the critical temperature for Bose condensa- 
tion found from g B £(3/2) = n B /^. = n B (2n/M B T c ) 3/2 , 
we get 



4nM I n B /n 



T ¥ (n) 



2/3 



M B \ 3^(3/2)/ 
j0.692 (n B — n ,M B = M) 
10.218 (n B =n/2,M B =2M) 



(152) 



where we assumed boson and fermion degeneracy to 
be g B — 1 an d gF = 2, respectively. The second case 
corresponds to a gas of composite bosons made of two 



fermions each. Since the ratio in Eq. ( 152 1 is of order 
unity, it is easier (yet nontrivial) to access bosonic su- 
perfluidity in ultracold quantum gases. In addition, due 
to Pauli blocking, the cooling of degenerate fermions is 
more challenging than for bosons. 
From Eq. 



(151 



we see that there is an exponential 
increase in T c for rising k?a, i.e. towards strong inter- 
actions. We may ask how this trend continues, and if 
it can help to achieve fermion pair condensation. More 
generally, one can imagine to tune the interaction pa- 



rameter kpa entering Eq. ( 145 1 from negative values 
through a resonance to positive ones. Such a knob in- 
deed exists in cold atom experiment thanks to Feshbach 
resonances. Not unexpectedly, the critical temperature 
doe not rise indefinitely when cranking up the interac- 
tions; in addition, for large positive values of kpa we 
will smoothly approach the limit of weakly interacting 
pointlike bosons. This is the BCS-BEC crossover. How- 
ever, the one-loop formula for the effective action used 
so far does not capture the crossover problem, because 
it relies on weak coupling and small amplitudes of the 
field. A more complete description is needed. 

4. Strong correlations and the Functional Renor- 
malization Group 

A full, quantitative grip on the strongly-correlated 
physics discussed in the previous chapter requires the 



use of nonperturbative techniques. Various methods 
have been applied to the physics of ultracold gases, 
ranging from quantum Monte-Carlo and diagrammatic 
Monte-Carlo Simulations to functional methods and 
resummation schemes such as Dyson-Schwinger and 
Kadanoff-Baym equations as well as the Functional 
Renormalization Group (FRG). The latter allows to ac- 
cess the whole phase diagram of ultracold gases in a 
unified approach, and, in particular, is also applicable 
in strongly correlated or strongly coupled regimes. Its 
setup and applications are described in the present chap- 
ter. This work does not aim at a fully detailed introduc- 
tion to the FRG. This has been done in various topical 
and general reviews and lecture notes, and we refer the 
interested reader to introductory and advanced general 
reviews [ 82-90 1, and to more topical reviews on nuclear 
and atomic physics Il7l [3414371 [914 - 19311 . non-equilibrium 
RG 194H961 . on gauge theories, QCD and QCD effec- 
tive models, I97ti99l . and quantum gravity lfT0W02"ll . 

4.1. Flow equation 

We begin with an introduction of the basic concepts 
of the FRG. It is based on the continuum version of 
Kadanoff's block-spinning transformations on the lat- 
tice 11031 . and has been formulated for the continuum 
by Wilson 1 1041 11051 . Its modern functional form for 
the effective action used in the present work has been 
put forward by Wetterich in 1 106]. 

For the description of ultracold atom experiments, the 
action S derived from the Hamiltonian in Eq. ( |T6] > is a 
microscopic starting point. It is related to an ultraviolet 
momentum scale A » Ciw' The relevant physics, how- 
ever, takes place at momentum scales k far smaller than 
A, and the respective quantum and thermal fluctuations 
have to be included. In the FRG framework, these fluc- 
tuations are included successively at a given momentum 
scale k starting at A with the action Fa = S , leading to 
an effective action r*. The latter already includes all 
quantum and thermal fluctuations above the momentum 
scale k. It can be interpreted as a microscopic action for 
the physics below the scale k in the very same way S has 
been introduced as the microscopic action of ultracold 
gases. After the inclusion of all fluctuations we arrive at 
the full effective action T, 



Ijfc=A - S, 

r k=0 = r. 



(153) 
(154) 



The effective action F^ interpolates smoothly between 
the microscopic (or initial effective) action Ta and the 
full effective action F = F^o- An infinitesimal change 
of the effective action with the scale k is described by 
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a flow equation for d k T k , which depends on the correla- 
tion functions of the theory at the scale k as well as the 
specific way the infrared modes with momenta smaller 
than k are suppressed. Such an RG-step has similarities 
to a coarse graining where details on short distances are 
continuously washed out, the difference being that the 
effective action F^ still keeps the information about the 
fluctuations between A and k. At the end of the process, 
for k — > 0, we include fluctuations with large wave- 
length. These are the problematic modes which cause 
infrared divergences in other approaches. Due to the 
stepwise inclusion of fluctuations, the renormalization 
group procedure is not plagued by such divergences. In 
conclusion, a given initial effective action Ta and the 



flow equation Eq. ( 165 1 define the full quantum theory 



analogously to the setting with classical action and the 
path integral. 

In this section, we derive the flow equation for T k and 
discuss its practical solution. To that end, we specify a 
suppression of low momentum fluctuations to, q 2 < k 2 . 
This is most easily achieved via a mass-like infrared 
modification of the dispersion relation, while the ultra- 
violet modes should remain unchanged: we add a reg- 
ulator or cutoff term AS k [(p] to the microscopic action 
S [(p] which is quadratic in the fields, 



S[<p]^S[<p]+AS k [<pl 



(155) 



The field ip is general and may be a collection of fields. 
For concreteness, we will use a notation analogous to 
ultracold bosons and write <p{X) = (tp(X)). We have 



AS*M= f <p\Q)R k (Q)<p(Q)- 

Jq 



(156) 



The requirement of the suppression of low momentum 
modes entails that R k (Q — > 0) + 0. In turn, for large mo- 
menta (in comparison to k), the regulator has to vanish, 
Rk(Q —* °°) —> 0. These properties can be summarized 
in the conditions 



lim R k (Q)=k 2 



lim R k (Q) = 0. 



(157) 



For the sake of simplicity, we have restricted ourselves 



in Eq. ( 157 1 to regulators that only depend on q and have 
a standard normalization R k (0) = k 2 in the infrared. The 
extension to general regulators is straightforward. 



If we interpret the action in Eq. ( 155 1 as the micro 



scopic action of a theory, it has a trivial infrared sector: 
The fields are gapped with gap k 2 . The generating func- 
tional of this theory is given by 



Z k [j] = \ Dipe 



-S-AS k +fj-tp 



(158) 



From Eq. (157i, we infer that Z[j] = Z k= o[j] is the 
full generating functional of the theory introduced in 
Eq. ( fJT) . For k — > A, the regulator term dominates the 
path integral as all physical scales are far smaller and we 
are left with a trivial Gaussian integral. Moreover, for a 
given k, the correlation functions {<p{Q\) ■ ■ ■ f(Q n )) tend 
towards the full correlation functions for Q 2 » k 2 for 
all i = \, ...,n. In turn, for Q 2 <sc k 2 , the correlation 
functions are trivial, as the fields are gapped. 

For explicit computations, it is more convenient to 
deal with the effective action r*, which is obtained via a 
modified Legendre transform according to 



logZ k [j] - AS k [<p] ■ 



(159) 



where j = j k [<p] satisfies (6 log Z k / 6 j)[j] = (p. We 
have already shown that the effective action has the 
simple physical interpretation of a free energy in a 
given background (p. Diagrammatically, it generates 
all one-particle irreducible diagrams. As in the case 
without regulator term, F k satisfies a functional integro- 



differential equation similar to Eqs. (48 I, (50i. Apply- 
ing the definitions of Z k and F k we find 



expf-r^]} 
= J Dip expj- 



■Sl<p + <p]-AS k [<p] + £ ^j-V} 

(160) 



where we have used the condensed notation introduced 
in Eq. (|78"]l, as well as 



m 



6(T k + AS k ) 



(161) 



following from the definition of the Legendre transform 



(159 1. Eq. (160 1 makes the suppression of the fluctua- 
tions even more apparent. Note first that the action S 
in the exponent depends on the sum + <p, whereas the 
cutoff term only depends on the fluctuation <p. Hence, 
for large cutoff scales k — > A, the functional integral in 



Eq. ( 160 1 gets Gaussian and the effective action tends 
towards the microscopic action, F^ A — > S . For k — > 0, 
the regulator vanishes, R k — > 0, and we are left with Eq. 

For a successive integration of momentum modes we 
need to know the flow d k r k . Applying the ^-derivative 



to Eq. (159 1 leads to 



d k r k [<f>] = -d k \.logZ k [j] - d k AS k [<f>]. 



(162) 



The notation signals that j is ^-dependent, but the 
terms proportional to d k j cancel. We have d k AS k [(f>] = 
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f Q dkRk(Q)<KQW(Q). The generating functional Z k 
only depends on k via the cutoff term AS*. Taking the 



^-derivative of Eq. ( 158 I, we can compute <9*|,logZ* to 
arrive at 

= j 9 k R k (Q) [(mwiQ)) k - </>(Q)<p*(Q)] ■ 

(163) 

Herein, we have restricted ourselves to bosonic fields 
(p. In the case of fermions, a global minus sign occurs 
due to the Grassmann nature of the fermions. The ex- 
pression in the square bracket in Eq. ( 163 i is the full, 



field-dependent propagator, which reads in terms of the 
effective action 



{<P(Q)<P*(Q)) - <f>(Q)f(Q) = 



l 



r< 2) + ** 



(2,0- (164) 



In Eq. (TT64J, we have used the property of Legendre 
transforms that the second derivatives of a functional 
and its Legendre transform are inversely related. In the 
present case, we note that the Legendre transform of 
log Z* is F* + AS *, as defined in Eq. 
are led schematically to 6 2 logZ k /5j 2 
which we have used in Eq. ( [164) . 
The momentum integral in Eq 



159 1. Hence, 
IT® + R k ) -- 



we 
1. 



163 i can be conve- 



niently written in terms of a trace. Including also the 
possibility of internal indices and different species of 
fields, we are led to the final expression for the flow 
equation for r*, 



d,T k 



1 



-STr 



1 



d,R k 



(165) 



the Wetterich equation. The supertrace includes the mo- 
mentum integration, and a summation over internal in- 
dices and field species, see also Eq. ( |54) , In Eq. \\65\ 
we have introduced the RG-time t = log k/ky with some 
reference scale ko, typically being either the ultraviolet 
scale, ko = A, or some physical scale. For a given quan- 
tity O k , the logarithmic scale derivative d,O k = kd k O k 
has the same properties under RG-scaling as the quan- 
tity itself. It also is convenient as one usually integrates 
the flow over several orders of magnitude in the mo- 
mentum scale k. Henceforth we shall use the standard 
choice t = log k/A. 

Above, we have argued that regulators with the prop- 



erties ( 157 1 lead to a suppression of the infrared physics 
of the theory. Moreover, since the finite initial effec- 
tive action Ta at the initial scale A already includes all 
fluctuations of momentum modes with momenta larger 
than A, no ultraviolet divergences should be present. 
These properties have to be reflected in the flow equa- 
tion ( |165) : It has to be both infrared and ultraviolet 




Figure 17: We plot a typical cutoff function R k (Q) = R k (3 2 ), which 
only depends on the spatial momentum. The function is nonzero for 
i} 2 < k 2 and thus provides an infrared cutoff for the propagators. For 
large momenta, it falls off rapidly, thus becoming inactive in the UV. 
The scale derivative R k (Q) = kd k R k (Q) is sharply peaked at 3 2 k k 2 . 
For this reason, the loop integral on the right hand side of the flow 
equation is dominated by these modes. This provides the mechanism 
how momentum shells are successively integrated out in the FRG 
framework. 



finite. Here, we show this explicitly for the case of 
bosonic fields. For low momenta, the regulator adds a 

(2) 

positive mass to F k in the denom inator. The typical 
size of this mass is k 2 , in Eq. ( 157 1 we have normalized 
it to k 2 . For the sake of simplicity, consider a classi- 



cal dispersion F 



(2) 



q 2 (with 2M 



1) for small 



momenta which tends to zero for vanishing momentum. 
Schematically, we have for small momenta 



1 



1 



rf(0 + fl*(0 iu + F + k 2 ' 



(166) 



which is finite for Q — > 0. For fermions, the infrared 
singularities arise close to the Fermi surface. Accord- 
ingly, the propagators have to be regularized there. In 
summary, this implies infrared safe flows. 

In turn, for large momenta, the scale-derivative 
d,R k (Q) vanishes due to Eq. ( 157) . If this happens suf- 
ficiently fast, 



Urn 4 2 d,R k (Q) 



0. 



(167) 



the momentum integral in Eq. ( 165 1 is finite. In the fol- 
lowing, we shall show results for regulators that satisfy 
Eq. 



167 1. We also remark th at m ass-like regulators, 
i.e. R k = k 2 , do not satisfy Eq. ( 167 1 and hence require 
UV renormalization. The related flows are functional 



Callan-Symanzik equations as first derived in 11071 . 
They are sometimes used due to computational simplic- 
ity, see e.g. [ 35 1 . The generic shape of a cutoff is shown 
in Fig. 



17 



It is apparent from the derivation that F* [(/>] depends 
on the shape of the regulator. This regulator-dependence 
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Figure 18: The flow of connects the microscopic action to the effec- 
tive action in the theory space of all possible action functionals. The 
latter is of infinite dimension since the effective action is character- 
ized by an infinite set of couplings (or correlation functions). This is 
indicated here schematically by the couplings gi,g2, g3 and {#,}. For 
different choices of the regulator R. , the trajectories in theory space 
differ, as is indicated in the figure. At k = 0, however, the particular 
paths merge again and eventually arrive at the full effective action. 



disappears for k — > 0, hence physical observables are in- 
dependent of the choice of R^, but the trajectory F^ from 
k — A to k — depends on R^, see. Fig. 



18 This leaves 



us with some freedom for the choice of the regulator. 
Indeed, its choice can be optimized to the approxima- 
tion under investigation ll8~7l 1 10814 1101 . In general, such 
a choice is further guided by computational simplicity, 
as in complicated systems the computational costs can 
be high. Typical choices are functions Rk(Q) which de- 
cay exponentially or even vanish identically for high Q. 
A slight complication for nonrelativistic system is pro- 
vided by the fact that frequencies and spatial momenta 
appear differently. The Galilei symmetric combination 
is given by ia> + q 2 , in contrast to the 0(4)-symmetric 
combination (q^) 2 = cl> 2 + q 2 for relativistic systems. At 
nonvanishing temperature, Galilei symmetry is broken. 
In the Matsubara formulation used in the present work, 
the coupling to the heat bath leads to periodicity in the 
imaginary time r with period Therefore, we may also 
choose a regulator which only depends on frequency or 
momentum space. Moreover, we may sacrifice Galilei 
symmetry in order to obtain simpler expressions for the 
flow equation. 

Here, we briefly discuss some common choices and 
their pros and cons. This should give the reader the 
chance to embark on first computations on their own 
while being aware of the advantages and limitations of 
the choices. A detailed discussion, however, is beyond 
the scope of the present work. 

For instance, exponential cutoffs for ultracold bosons 



are given by 



RkiQ) = 



(168) 



The power of n can be chosen such that the cutoff falls 
off sufficiently fast for high Q. We remark that despite 
its exponential nature, one has to choose n > 1 in order 
to have flows which are local in momentum space, see 
01111 . This is mirrored in the property that for n — 1 
the cutoff insertion d,Rit{q 2 ) is not peaked at about q 2 » 
k 2 but is a monotonously decaying function. Only for 
larger n do we get peaked cutoff insertions. We also 
remark that rapid decay of R^ as a function of a> poses 
problems for computing thermodynamic quantities, as 
they lead to oscillations in the flow. These oscillations 
mirror the property that the periodicity in /3, even though 
present for smooth regulators, is lost for a sharp cutoff 
in to. 

A particularly useful cutoff is the three-dimensional 
Litim cutoff ifTMl fTTIl . 



RkiQ) = (k 2 



■ q 2 )6(k 2 



->2\ 

q ), 



(169) 



which effectively reduces the momentum integral to 
q 2 < k 2 and replaces q 2 by k 2 . It facilitates the analytic 
derivation of flow equations for correlation functions in 
the derivative expansion, and hence leads to important 
computational simplifications. Its analytic property also 
allows an easy access to the structure and interrelation 
of the flows (and hence the correlation functions). These 
properties make it the standard choice within (lower or- 
ders of) the derivative expansion. Moreover, in three- 
dimensional theories, the cutoff in Eq. ( 169 1 provides 



an optimal choice 11871 11081 within the lowest order of 
the derivative expansion scheme as will be employed 
below. In (3 + l)-dimensional theories, it still shares 
some of the stability features it has in three dimensions 
HI 12L but loses locality in momentum space necessary 
for full quantitative precision 111 111 . 

A manifestly Galilei symmetric regulator is provided 
by Hm, 

RkiQ) = t^I' (170) 

i +C (^r 

where n determines the algebraic decay for large mo- 
menta and c is a prefactor of order unity. Eq. ( 170 » 



can be extended to more general rational functions in 
the Galilei invariant \co + q 2 . Its key advantages are its 
Galilei invariance as well as its analytic structure. The 
latter allows to continue the results to real time, and 
hence may give access to transport properties or more 
generally dynamics of ultracold gases. Similar choices 
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in relativistic theories can be used for computing decay 
properties [ 1 1 3 1 . 

Note that the flow equation ( |165| l has a one-loop 
structure, which can be traced back to the quadratic 
form of the regulator in Eq. ( 156 1. Indeed, we may 



rewrite Eq. ( 165 1 as 



d t T k [<Pl = ^STr 3, log (rf[0] 



(171) 



where the derivative d, only acts on the ^-dependence of 
the regulator, i.e. we have 



d t = d \ T ay 



(172) 



In Eq. ( 171 1, we identify the one-loop formula for the 
effective action ( |54| on the right hand side. Therein, we 
have to substitute 5 (2) with the full two-point function 



Eq. ( |171| i is a very convenient form of Eq. ( 165 1 for 



deriving the flows for correlation functions, e.g. <5,r^"\ 
It also allows for an easy access to the fluctuation- 
dependence of specific correlation functions without 
performing any calculation. First we note that the flow 
equation for any correlation function (e.g. a coupling 
constant) can be obtained by writing down all one-loop 
diagrams which contribute to the expression, replace the 
propagators and vertices by full ones, 



1 

5(2) 



1 



r? ] + R t 



,<») 



(173) 



approximation schemes to the full effective action Ta, 
which include the physics at hand already at a low order 
of the approximation. The systematics of a given ap- 
proximation scheme and the control of the related sys- 
tematic error is of chief importance when it comes to 
the discussion of the reliability of the results. A discus- 
sion of this interesting point is beyond the scope of the 
present work, and is tightly linked to the discussion of 
optimal choices of regulators mentioned above. 

Here, we briefly discuss the most important approxi- 
mation schemes which cover (in variations) all approx- 
imation schemes used in the literature. The most im- 
portant scheme, which is partially behind all approxi- 
mations used, is the vertex expansion about a specific 
background (/>, schematically written as 

1 f " 

nw = 2_i - *t\4>Wu...,x N ) Y] [<p(Xt) - 0(Xi)] . 

(174) 

The information about the effective action is encoded in 
the vertices r[ n) . The related flow can be derived from 
that of the effective action as 



5" 



6(KXi)...6<f>(X„) 



d t T k m . 



(175) 

On the right hand side, we have to take the nth derivative 



of the one-loop diagram in Eq. ( 165 i, 



S" 



1 



<5 < *(X 1 )...^(X„)2 



STr 



1 



d,R k 



(176) 



and then take the 5,-derivative. Remarkably, this ren- 
ders the one-loop expression an exact (flow) equation. 
We emphasize that this only holds true for additive IR 
regularizations of the one-loop formula for the effective 
action (54) , see 11 1 1411 . Note also that one also has to 
take into account perturbative one-loop diagrams where 
the vertices involved vanish classically, = 0. Still 
this one-loop structure is very useful: If the loop expan- 
sion of a given correlation function does not exhibit a 
one-loop diagram, this correlation function is not sen- 
sitive to quantum fluctuations. This either happens due 
to internal symmetries or the pole structure of the dia- 
grams. The latter is specific to nonrelativistic theories 
and, as we will see below, in the case of ultracold atoms 
it leads to strong simplifications for vacuum scattering 



properties (cf. Appendix E i 



The Wetterich equation (165i is an equation for a 



functional and thus may be evaluated for any (possi- 
bly inhomogeneous) mean field <p(X). It is a functional 
integro-differential equation and its full solution is, in 
most theories, beyond reach. Instead, one has to use 



which produces all possible one-loop diagrams with 
cutoff insertions. Evidently, the diagrams for the flow 
of TT depend on IT 1 ' with m < n + 2. Hence, within 
the vertex expansion described above, we arrive at an 
infinite hierarchy of equations, because the flow equa- 
tion for r (,l) requires input from r (n+I) and F ( " +2) . The 
flow of the latter two quantities depends again on higher 
correlation functions and eventually the system never 
closes. We should not be surprised about this, as the 
effective action necessarily contains infinitely many in- 
dependent terms, and we have just rewritten the func- 
tional integro-differential equation in terms of infinitely 
many partial integro-differential equations. In most in- 
teresting cases, it will not be possible to derive a closed 
expression for the functional T[<p], Practically, one trun- 
cates the hierarchy of flow equations at a given order 
n, i.e. approximates r* m>/,) ~ o, and solves the re- 
stricted, finite set of partial integro-differential equa- 
tions for r!! m (Qi, Q m ). Examples for this scheme 
can be found in e.g. 1711631 l75¥77l fTBHTTSl The self- 
consistency of this approximation can be checked by 
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computing the flow d t TT l>n ' as a function of r?" . This 
provides some error control. 

A further important approximation scheme is the 
derivative expansion. Formulated in momentum space, 
it is an expansion of the vertices in powers of the mo- 
menta (derivatives). Its nth order relates to the nth- 
order in \cl> + q 2 . In contrast to the vertex expansion, 
all vertices are present already at the lowest order of the 
derivative expansion. Here, we exemplify this expan- 
sion for the case of the effective action r* of a Bose gas. 
An often used ansatz for this theory is provided by 

T k [f, 0] = ^{f(Z k d T - A k V 2 )<f> + U k (f<f>j) . (177) 

Herein, U(<p*(j)) is the full effective potential. It is a gen- 
eral function of p — <p*(f). Accordingly, we have t/ ( "- ) + 
and thus vertices rjj" ) to all orders in n. For Z k — A k = 1, 



the ansatz in Eq. ( 177 i has the same momentum de- 



pendence as the classical action and is the lowest (or 
zeroth) order in the derivative expansion. For flowing 
Z k ,A k , one goes beyond the lowest order. Note however 
that it is not the full first order in the derivative expan- 
sion, as this requires field-dependent Z k (p),A k (p). The 
derivative expansion and in particular the above ansatz 



( 177 1 assumes lower orders of the differential operators 



to be more relevant than the higher ones. In the pres- 
ence of a mass gap m gap , this is expected to be valid in 
the infrared, because 



\q/m 



i -a 



gap I 



0. (178) 



Hence, within the derivative expansion, we make an ex- 
pansion of the effective action about the low energy ef- 
fective action. From the more technical point of view, 
we project the flow of onto a subspace of function- 
al, which mimic the shape of the microscopic action. 
The latter is given by Za = 1, A\ = 1, and the effective 
potential Ua(p) = -yip + \p 1 - These values constitute 
the initial conditions for the flow equation. During the 
flow, dressed quantities Z k , A k and U k (p) emerge. Be- 
sides the ansatz in Eq. ( 177 1, we also have to specify a 



projection description which deter mine s the flow equa- 
tions Z k , A k and U k (p) from Eq. ( 165 i. Examples for 
the full lowest order derivative expansion in bosonic as 
well as mixed fermionic-bosonic theories can be found 
in e.g. j84ll99lfTT9l. 

Most applications to ultracold atoms discussed in the 
present work are done within low orders of the deriva- 
tive expansion within an additional field expansion of 
the effective potential U k up to the nth order of the 
fields. Of course, such an approximation can also be 
interpreted as the nth order of the vertex expansion with 



an additional expansion in powers of momenta and fre- 
quencies. Indeed, as has been mentioned before, any 
approximation scheme used in the literature can be seen 
as combination, deformation or further approximation 
of the vertex expansion and the derivative expansion. 
In any case, when using such an approximation, we re- 
strict the space of functionals. For this reason, although 
we started from an exact flow equation, we may accu- 
mulate errors. In particular, given the exact flow equa- 
tion, every regulator satisfying the mentioned properties 
should give the same result. But since we never inte- 
grate the full flow, we may end up at different "effective 
actions" r,t = o if we used two different regulators. This 
regulator dependence can be applied to partially check 
this source of uncertainty. The approximate indepen- 
dence of the results at vanishing cutoff, k = 0, guaran- 
tees the self-consistency of the approximation. In turn, 
the independence of r*=o of the chosen regulator R k or 
the chosen trajectory in theory space can be utilized for 
devising regulators that are best-suited (optimal) for the 
given order of a given approximation scheme at hand, 
see 11871 [TOfTT0l[T2UHT23l . 

In summary, the Functional Renormalization Group 
approach for the effective action constitutes a fully non- 
perturbative approach to quantum field theory. In fact, 



the functional differential equation (165i may be seen 



as an alternative but equivalent formulation to the func- 
tional integrals (36 1 or ( 50 » in Sec. 3.1 In partic- 



ular, it is neither restricted to small couplings nor to 
small amplitudes. For this reason, it can be applied 
to many strongly coupled and/or strongly correlated 
systems such as superconductors, superfluids, quantum 
chromodynamics, quantum gravity, or - in our case - 
the unitary Fermi gas. 

4.2. The many-body problem in ultracold atoms from 

the FRG point of view 
We are now ready to proceed to the application of the 
FRG to the physics of ultracold atoms. Particular em- 
phasis is put on the difference between the vacuum and 
the many-body limit of the system, which can both be 
accessed in experiments and thereby allow for a high 
precision comparison between experiment and theory. 
For this reason, cold quantum gases also provide an 
ideal testing ground for different approximation and 
truncation schemes within the FRG approach. Our con- 
siderations build the basis for the detailed analysis of 
the BCS-BEC crossover of two-component fermions in 
the subsequent section, but are generally applicable to 
all ultracold atom settings. In addition, it will provide 
an RG point of view on UV renormalization discussed 
in Sec. 
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The microscopic action — S of an ultracold Bose 
jas at the initial ultraviolet scale A is given by 



S[<p*,<p] 



V2_ 

2M 



(179) 



As an effective action, this expression is only valid in 
the UV, which is given here by the van der Waals length. 
Going to the low energy scales realized in experiments 
with ultracold atoms, fluctuations are included and the 
parameter g^ gets replaced by a dressed quantity g. The 
corrections to the UV value arise from quantum and 
thermal fluctuations. 

Assume we want to measure the scattering length. 
Say, we perform a scattering experiment between two 
atoms such that there is no influence from other parti- 
cles. Hence, there are no effects which are associated to 
nonzero temperatures or densities in this setting. There- 
fore, we call such processes to take place in the vacuum, 
because they could also be observed if we had nothing 
but the two colliding partners. In practice, this situation 
is achieved in ultracold experiments for T — > 0, but with 
vanishing phase space density a> = (Ajj€) d <sc 1 (cf. the 



discussion around Eq. (22 1). 

The scattering length a of two identical bosons is re- 
lated to the coupling constant according to 



g = Snh 2 a/M . 



(180) 



It is important to realize that the dressed or renormalized 
coupling constant g appears in this equation. Although 
the scattering takes place in the absence of statistical 
many-body effects, we still have quantum fluctuations, 
which are always present. For this reason, the measured 
coupling constant g does not coincide with the bare cou- 
pling constant g^ of the microscopic action. 

In fact, the experimental relevance of the bare pa- 
rameters in Eq. (179 1 is at best indirect. Our micro- 



scopic action is not a realistic description for atoms at 
high energies, but rather a simple model with the same 
low energy physics as a more elaborate description. Put 
differently, we utilize here that only the renormalized 
parameters are important for observations in cold atom 
experiments. 

It is a key property of ultracold quantum gases that 
the interaction parameters of the atoms can be measured 
without reference to the many -body system. In our case, 
we know the value of the scattering length a and can use 
this as an input for many -body predictions. For instance, 
we may express the equation of state at zero temper- 
ature as a function of a and the chemical potential fi, 



P = P(ji, a), and verify experimentally the predicted de- 
pendence on both fi and a. This optimal situation is not 
met genetically in condensed matter systems like e.g. 
solids, where the parameters of the model Hamiltonian 
are not known and have to be adjusted according to the 
observed many -body physics. In addition, a systematic 
investigation of interaction effects is often not possible 
in solid state physics, because the parameters are fixed 
by the sample and cannot be tuned arbitrarily. 

The "dictionary" between bare and dressed parame- 
ters of the microscopic action in vacuum is sometimes 
referred to as the UV renormalization. We already en- 
countered such a procedure in the treatment of weakly 
interacting bosons and fermions in Sec. [5] Within the 
FRG approach, the notion of UV renormalization of mi- 
croscopic parameters arises naturally. 

The microscopic action in Eq. ( 179 1 enters the flow 
equation ( 165 i for the effective action. Depending on 



how we choose the initial parameters fi and T of the 
microscopic action, we will arrive for k — either in 
vacuum or at a many-body system. Herein, the vacuum 
effective action is defined through a diluting procedure 
n, T — > 0, such that the system is kept above criticality, 
T/T c (ri) > 1. In this way, condensation is excluded and 
we end up in the physical vacuum of a few scattering 
particles. An equivalent way to express the same idea is 
to take the limit T — > 0, and at the same time sending 
the phase space density u = (Ar/€) d — > (n = t~ d , A T ~ 
l/yff, cf. Eqs. |l}and Q). In other words, the system 
gets dilute faster than it gets cold. Indeed, the vacuum 
scattering experiments work precisely in this ultracold 
limit at low phase space density. These considerations 
are summarized with the definition 



lr-»0, E)=nAi-*0' 



(181) 



Hence, in this case, F^ is in the symmetric phase for all 
k and tends towards that at n — T — for k = 0. 

This allows us to set-up the solution of the many- 
body problem with the FRG in a two-step procedure: 

(1) Solving the flow in vacuum 

Choosing fi and T such that r* = o = r vac , we start 
from a given set of bare parameters and compute 
the resulting renormalized ones via a successive in- 
clusion of quantum fluctuations. The renormalized 
coupling can be measured and tuned in practice, 
and thus allows for a direct matching of theory and 
experiment. 

(2) Solving the flow for the many-body system 

We set the chemical potential and the temperature 
such that we arrive at a desired n > and T > 
for k = 0. 
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r A = s 




Figure 19: For given /j and T, the effective action V^-q = V will ei- 
ther describe a few-body or a many-body system. We denote these 
two cases by r vac and r m b in the plot, respectively. To compute ob- 
servables for ultracold quantum gases we need both: First, we solve 
the vacuum case, which provides us with the renormalized parameters 
of the microscopic action. Then, in a second step, we switch on the 



many-body scales k m \, 



(inverse interparticle spacing and 



de Broglie wavelength), which influences the trajectory of F/, once 
k ~ k m \,. 



Let us now follow the integration of the flow equa- 
tion at finite density and temperature. For large k with 
T /k 2 ,fi/k 2 <K 1, the flow agrees with the vacuum flow 
in the symmetric regime up to subleading contributions. 
Once the flow parameter k reaches the many-body mo- 
mentum scales fi l ^ 2 ,T 1 ^ 2 , the flow deviates from the 
vacuum trajectory and eventually arrives at the effec- 
tive action T of the many-body system. Having started 
with a given set of bare couplings, we know from step 
(1) what the corresponding renormalized couplings are. 
Thus, we can express the many-body observables in 
terms of n, T and physical microscopic couplings. The 



situation is illustrated in Fig. 19 



A simple example for this two-step procedure is pro- 
vided by the microscopic action for bosons given in Eq. 
(179 1. Herein, g A is the only bare parameter which 
receives substantial renormalization. In order to com- 
pute the renormalized coupling, we follow the flow for 
T = n — and g A > down to k — 0. The renormal- 
ized coupling g - g(gA) can then be read-off from the 
vacuum effective action F^o. This completes step (1). 
Now, we again solve the flow equation but with n > 
and T > 0. For instance, we may compute the equa- 
tion of state n(jj., T, g A ). From step (1) and Eq. 



180 1 we 



can replace gA by the scattering length and obtain the 
equation of state in the experimentally accessible form 
n(jj, T,a). 



4.3. BCS-BEC crossover and unitary Fermi gas 
In this section, we investigate the BCS-BEC crossover 
of two-component ultracold fermions. Applying the 
FRG in a simple truncation, we capture both the weakly 
and the strongly interacting regime. In particular, we 
show how the BCS and BEC ground states are linked 
by the unitary Fermi gas. Furthermore, we are able to 
resolve physics at all scales, from the few-body physics 
over the many-body sector down to very long wave- 
lengths, relevant for critical behavior. One of the goals 
of this section is to show how the FRG provides a uni- 
fied approach which captures all effects of the system 
within the same description. 



Feshbach resonances and microscopic model 

We have discussed the two cornerstones of quantum 
condensation phenomena in the weak coupling regimes: 
On the one hand, attractive interactions lead to super- 
fluidity of two-component fermions via the formation 
of Cooper pairs. The momenta of two fermions consti- 
tuting such a pair are located on opposite points of the 
Fermi surface. This locality in momentum space im- 
plies that the spacing between them may be large in po- 
sition space. On the other hand, we discussed Bose con- 
densation of weakly repulsive bosons, which are micro- 
scopically pointlike objects, localized in position space. 
Tightly bound pairs of two fermions could effectively 
realize such bosons. 

There exists an experimental knob to connect these 
two scenarios. It is provided by the Feshbach reso- 
nance [39], which allows to change the scattering length 
through the variation of an external magnetic field B. 
We write 

a(B) = a b A l-^—^-j, (182) 

where flb g , AB and Bq are parameters which can be 
determined experimentally. This formula is a decent 
parametrization in a range of order AB around B « Bq. 
In particular, at B = Bq, the scattering length changes 
sign and becomes anomalously large, \a\ » Adw- Re- 
call that this does not invalidate our effective Hamilto- 
nian and the fact that scattering can be assumed to be 



pointlike, see the discussion at the end of Sec 2.1 



Sufficiently stable ultracold quantum gases of two- 
component fermions are built from either 6 Li or 40 K, 
which are alkali atoms. Their internal structure is rel- 
evant in order to have fermionic s-wave interactions at 
all, but also manifests itself in the appearance of Fesh- 
bach resonances. To understand this, we consider a 
single alkali atom. We can approximate the system to 
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consist of an atomic core and a valence electron. The 
ground state of the system is given by the electron be- 
ing in the s-orbital. Correspondingly, the orbital angular 
momentum of the valence electron vanishes and thus a 
fine structure does not appear. However, the electron 
spin S couples to the spin of the nucleus /. The re- 
sulting quantum number F introduces a (tiny) hyperfine 
splitting of the ground state. Since S = 1/2, the value of 
F is given by F = I ± 1/2. In addition, every hyperfine 
state has a (2F + l)-fold degeneracy m F = -F, . . . , F. 
Thus, alkali atoms in their electronic ground state can be 
distinguished according to their hyperfine state \F, mp). 

Now suppose that two atoms in different hyperfine 
states scatter off each other. Due to the internal struc- 
ture of the colliding partners, we call this a multichan- 
nel scattering. The two-body system of atoms will be in 
a superposition of the singlet and the triplet state. De- 
pending on the species of atoms, the former will have a 
higher or a lower energy than the latter, while the first 
option is more generic. Moreover, due to the hyper- 
fine coupling there will in general be a mixing between 
both states. For our purpose it is enough to restrict our 
considerations to two relevant channels, an open and 
a closed channel, which have different magnetic mo- 
ments. We normalize the potential such that two atoms 
in the open channel at infinite distance have zero energy; 
this sets the scattering threshold. The closed channel is 
separated from the open one by a large energy gap AE 
(cf. Fig. 20 1. It thus cannot be accessed by atoms in 
the lower channel. The relevant feature of the closed 
channel is a bound state lying close to the open channel 
scattering threshold. It is evident that this situation is 
not particularly generic since typical bound state level 
spacings are much larger than typical collision energies 
in ultracold gases, and thus requires specific, fine tuned 
conditions. 

Because of the difference Ap. in magnetic moment, 
open and closed channel couple differently to an exter- 
nal magnetic field B. For this reason, the difference in 
energies between both channels can be tuned according 
to AE — > AE + Ap ■ B. Consider a particular bound state 
from the closed channel. Its energetic distance from the 
scattering threshold E = is called the detuning 



v(B) = An ■ {B - B„) • 



(183) 



Due to second order processes, where two colliding 
atoms virtually enter the closed channel and then leave 
it again, a bound state with small v affects the scatter- 
ing properties of the alkali atoms. In particular, chang- 
ing the magnetic field such that v — > 0, both channels 
become resonant and we obtain a strongly interacting 




open channel 



Figure 20: Interatomic potential U between two fermions in distinct 
hyperfine states separated by a distance r. The closed channel consists 
of bound states and low energy scattering can only take place in the 
open channel. However, by changing the external magnetic field B, 
we can drive one of the bound states close to the scattering threshold 
U (r — > oo) = 0. The energy distance related to this particular bound 
state is denoted as v(B) = Ap ■ (B - Bo). The resulting scattering 
length a = a(B) is parametrized according to Eq. {182). For B x 
Bq, it becomes anomalously large, and thus it can largely exceed the 
interparticle spacing: \a\ 3> £ v( w => Kifepa) I $ !• 



system. The corresponding field dependent scattering 



length a(B) in Eq. ( 182 1 constitutes the most important 
application of such a Feshbach resonance lfTll2"l ll24l in 
ultracold atomic physics. 

We now incorporate the physics of a Feshbach reso- 
nance for two-component fermions on the level of the 



microscopic action. Recall from Eq. ( 120 1 that the ac 



tion for a fermionic theory with pointlike interactions is 
given by 

S t [f, iff] = §{f(d T - ^)<A + ^%V) 2 ) , (184) 



with Grassmann fields iff = (if/\,tf/2). Eq. (184i consti- 
tutes a single-channel model. In order to emphasize the 
structure of this and the following expressions, we have 
dropped the chemical potential for the moment. It will 
be reintroduced below. 

The closed channel is included explicitly in terms of 
a microscopic bosonic field <p, which constitutes a com- 
posite degree of freedom resulting from the intercon- 
version of two fermions into a closed channel molecule. 
The action for this boson field is modeled as 



S<pW,<P\ 



4M 



<P- 



(185) 



The most important term is the detuning v, which acts 
as a mass term for the bosons. In addition, we allow 
for a Galilean invariant kinetic term, where the prefac- 
tor of 1 /AM is related to the mass 2M of the composite 
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object. As we will see in a moment, the microscopic ki- 
netic term is, however, unimportant for the case of broad 
Feshbach resonances which are studied here, and could 
be equally well omitted. The full microscopic action 
from which we will extract the physics of the BCS-BEC 
crossover is then given by 11251 

s [p, «a, v \ <p\ = J ty\a T - + ^%V) 2 

+ ^*' <5t ~ 4M +V r~ %>i^2 ~ )• (186) 

As anticipated above, the Yukawa-type cubic coupling 
~ h(p*t//i4^2 (called Feshbach coupling in the cold atom 
context) allows for the interconversion of two fermions 
of opposite spin into one molecule. 

The parameters Ab S , v and h of the microscopic action 
can be measured in experiment. We show here that they 



correspond to the three parameters in Eq. ( 182 1 for the 
scattering length a(B) across a Feshbach resonance in 
the broad resonance limit. For this purpose, we consider 
the functional integral Z = J D<pDi//e~ s W , ' ip \ For fixed tfr 
and i}r m , we can perform the Gaussian integral in <p* and 
(p. This is equivalent to the saddle-point approximation 
about the solution of the mean field equations of motion 



6S 
6^ 



d T - V 2 /4M + v 



ifnfo. (187) 



As the action is quadratic in tp, the saddle-point approx- 
imation is exact. If we formally insert this into the par- 
tition function and integrate out the bosonic fields, we 
arrive at the action 



Sl ' J/ •' Pl=S " + i^ 2 d 7 -^/4M + r 



W2, (188) 



with from Eq. (184i. We emphasize that the 



procedure described here corresponds to reversing a 
Hubbard-Stratonovich transformation explained in Sec. 
|3.4| for a slightly more complicated inverse boson prop- 
agator; this is possible due to the fact that in Eq. ( |185| l 
we work with a quadratic bosonic action. We now take 
the broad resonance limit, where h, v — > 00 with h 2 /v 
kept fixed. Then, we can neglect the derivatives corre- 
sponding to frequency and momentum dependence of 
the effective four-fermion vertex. More precisely, we 
scale v ~ h 2 for h — > 00, while leaving the derivative 
coefficients of order unity. We then obtain the action 



Sty*,*] =SW,ifr] 



\_h 2 
2 



v Jx 



(<AVr 



(189) 



Apparently, this coincides with a purely fermionic ac- 
tion with an effective coupling 



h 2 

^eif - As • 

V 



(190) 



We conclude that in the broad resonance limit the two- 
channel and the single-channel model become equiva- 
lent Il6l ll26l[l27l . The single channel model, however, 
acquires an additional effective contribution to the cou- 
pling constant. The Feshbach resonances in 6 Li and 40 K 
are broad. The narrow resonance limit is conceptually 
interesting, as it can be solved exactly [ 126 1. Moreover, 
recently, examples of narrow resonances have been 
studied experimentally II 1281 . In an RG language, it cor- 
responds to a Gaussian fixed point, while the broad res- 
onances are governed by an interacting (Wilson-Fisher) 
fixed point I129L While the macroscopic physics de- 
pends on microscopic details of the closed channel, the 
broad resonance fixed point is characterized by a large 
degree of universality, i.e. a pronounced independence 
on the microscopic details of the closed channel, as is 
plausible from the above scaling [129, 1 30 1 . 



Assuming the couplings in Eq. ( 190 1 to be the renor- 
malized ones, we can relate them to the scattering length 
according to A — 47ra/M0Wefind 



M ( h 2 



(191) 



Comparing this to Eqs. ( |182| i and ( |183| l, we find that 
indeed v - Afi ■ (B - Bq) corresponds to the detuning 
from resonance. The four-fermion coupling Ab g is re- 
lated to the background scattering length in the usual 
manner via a bg = MAh g /4n. With AB = h 2 Ap./A bg , the 
Yukawa/Feshbach coupling h is seen to determine the 
width of the resonance. 

For magnetic fields close to Bo, the scattering length 
becomes anomalously large and the background scatter- 
ing length can be neglected. In what follows, we assume 
Ob g = throughout the whole crossover. 

As anticipated above, only the value of the scatter- 
ing length plays the role of a relevant parameter for 
the crossover in the broad resonance limit. Given the 
density n of atoms, we build the dimensionless param- 
eter (kpd) . Since the interparticle spacing is given by 
( =s k F l , we find the following scheme for the crossover: 

1) kpa — > -00 : weakly interacting fermions, 



3 Note the difference in convention to identical bosons, where A 
8na/M. 
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2) \(kpa) l \ < 1: strong interactions, dense regime, 



3) Mr 



oo: weakly interacting molecules. 



The regions a < and a > are called BCS and BEC 
side of the crossover, respectively. 

Region 2) is often referred to as the unitary Fermi 
gas. The origin of this term is the following. The cross 
section of two-body scattering in the s-wave channel 
is given by o";=o = 47r|// = o| 2 . For the perturbative re- 
gions with p\a\ <K 1, where p is the relative momen- 
tum of scattering particles, we then find <x = Ana 2 . 
Naively extrapolating this to the resonance \a\~ l — > 0, 
this would imply a divergent cross section, which is ex- 
cluded from the fact that the scattering matrix is unitary. 
Recalling however Eq. (fTTh, we find in the latter limit 
that <t a An\{— I /a - ip)\^ — > 4-n/p 2 . Since scatter- 
ing is meaningful only for nonzero relative momenta, 
the expression on the right hand side constitutes the up- 
per limit on possible s-wave scattering; in the unitary 
Fermi gas, the typical scale for the scattering momen- 
tum are kp, yT. This effect has been observed in II 1 3 1 H . 
We note that exactly at the unitary point aT l = 0, the 
scale associated to interactions drops out and the only 
remaining scales are interparticle spacing and temper- 
ature. This hints at highly universal properties in this 
regime II 1321 (distinct from the broad resonance univer- 
sality described above). However, at this point and in 
its vicinity where l^pflT 1 <s 1, the gas parameter Eq. 
( p"2| ) is large and cannot be used to as a control param- 
eter for systematic expansions. In this regime, where 
the interaction length scale greatly exceeds the interpar- 
ticle spacing, we deal with a strongly coupled and dense 
quantum system. 

In the perturbative regime of small a, we have a sec- 
ond order phase transition towards the superfluid phase. 
We will see that a second order transition is found for 
all values of kpa. The value of the critical temperature 
T C /T F as a function of k?a is particularly interesting. 
Deep in the BEC regime, we expect the noninteract- 



ing formula (22 1 to hold, with a shift due to the small, 



but nonvanishing diluteness parameter kpa. On the BCS 
side, Eq. ( |152| > will turn out to be insufficient, because 
it is lowered by a factor of approximately two due to 
particle-hole fluctuations. A great challenge in many- 
body theory is the calculation of T c /Tp at unitarity from 
first principles. 



Ansatzfor the effective action 

After these preliminaries we turn our attention to an 
FRG study of the crossover. Since this approach is not 



limited to weak coupling, it can be applied to the lim- 
iting BCS and BEC regimes as well as to the unitary 
Fermi gas. Our particular interest lies in the transition 
from micro- to macrophysics. 

We restrict our analysis to the three dimensional case. 
However, performing the same calculations with gen- 
eral dimension d of the loop integral, we can also ana- 
lyze lower-dimensional systems. The scattering proper- 
ties have to be adjusted appropriately. We expect quan- 
tum effects to be more pronounced in reduced dimen- 
sionality, because long-range effects are more signifi- 
cant there. This is reflected in more severe infrared sin- 
gularities in the respective loop integrals. We also re- 
strict to the balanced case of equal densities of spin up 
and down atoms. For the imbalanced case, a first order 
phase transition is expected and thus its implementation 
is more demanding. The latter corresponds to unequal 
chemical potentials p\ + pi for the fermions. We set 
2M = 1 with the fermion mass M. 



The chemical potential which enters Eq. (186 1 can 



be found from the following consideration. Since every 
bosonic molecule consists of two fermionic atoms t//, 
we have for the total number of atoms N - N^, + 2N$. 
For this reason, what appears in the grand canonical par- 
tition function is the combination H - pN - H -p(N^, + 
2N$). On the other hand, we can also regard the system 
to be composed of two species with individual chemi- 
cal potential, i.e. H - p^N^, - p^N^ enters the partition 
function. Of course, both expression have to be equal, 
which yields p^ — p and p,p = 2p. The full microscopic 
action r A = S thus reads 



S [ip, ifr] = J {tf(d T - V 2 - p)if, + ^A(iAV) 2 



+ <p*(z A d T - ^A A V 2 + v A - 2p)y. 



2 



((p^si/f-yifSsifr*)). (192) 



The bosonic mass is related to the detuning from the 
Feshbach resonance as explained in the last section. 

Already a simple approximation for the effective ac- 
tion allows for a qualitative study of the full crossover. 
Since the key to the BCS-BEC crossover consists in the 
formation of a bound state, the possibility of describing 
a bosonic molecule must be contained in any reasonable 
approximation. We thus employ the following ansatz 
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for the effective action, 



+ f(z k d T - \A k V 2 )cp + U k (p) 



(193) 



As explained in Sec. 4. 1 this approximation, up to the 
wave function renormalizations Z k ,A k , constitutes the 
lowest order in the derivative expansion. 

The ansatz in Eq. (193i forces the effective ac- 
tion into a particular form. The general fluctuation- 
dependence of the vertices H is replaced by the flow of 
Z k , A k , h k and U k (p). This simplified picture, however, 
already encodes the most relevant physics of the theory. 
For instance, the fermion propagator is considered clas- 
sical in this approximation. The inclusion of the full 
fermion dynamics yields higher quantitative precision, 
but is not required to obtain a qualitative picture of the 
phase transition, which is driven by the particle-particle 
loop contribution to the boson propagator as we will see 
below. It is in this sense that a successive improvement 
of the truncation of the effective action can lead to new 
physical insights. 

The scale dependent effective potential U k (p) can 
only depend on the U(l)-invariant quantity p - 
<p*d>. This follows from our considerations in Sec. 



3.2 The U(l)-symmetry of U k is not violated dur- 
ing the renormalization group flow, if the regulator 
AS k respects this symmetry. As we employ AS k = 

Jg <P*(QWk(Q)<p(Q) + 4^(0^(0^(0, which is 

manifestly U(l)-symmetric, this is the case. It is a par- 
ticular strength of the FRG approach that symmetries of 
the theory are conserved if the truncation and the regu- 
lator are chosen appropriately. The initial condition for 



U k (p) can be deduced from Eq. ( 192 1 and is given by 



Ua(p) = (va - 2fi)p . 



(194) 



Note that the microscopic potential does not contain a 
term ~ u^^P 2 with dimer-dimer coupling Ua\ in accord 
with our discussion of the Feshbach model ( |186| l. De- 
spite being zero at k — A, this coupling will be gener- 
ated during the RG flow within our truncation. 



In Eq. ( 193 i, we have neglected a possible running 
of Apjc. Certainly, for A^ t \ < 0, we can eliminate the 
four-fermion vertex in the ultraviolet by choosing h A 
appropriately, see Eq. (132i. However, during the RG 



flow, the coupling A^k is generated again. By neglecting 
the flow of this coupling, we simplify the flow equations 
but miss an important screening effect in the many-body 
problem. 



Another interesting feature is the irrelevance of the 
initial conditions for Z\, A A and h A in the broad reso- 
nance limit. First of all, one may argue that the coupling 
A A is actually not present, because it can be absorbed 
into the definition of 0^] We apply this redefinition pro- 
cedure at all scales k by introducing the renormalized 
field 

= A' /2 0- (195) 



As a consequence, the effective average action is given 
by 

+ ~4> t {s k d T - l -v 1 )^ + u k (p) 

-^(^V^-^V)), (196) 



where we have introduced the renormalized couplings 
S k - Z k /A k and/i^ = h 2 /A k . The designation 'renormal- 
ized quantities' is common but somewhat unfortunate, 
since the actual renormalization procedure takes place 
in the flow with k. Note that Eq. ( |196[ ) has a canoni- 
cal kinetic term for the bosons without prefactor of the 
V 2 -term at all scales. 

As mentioned below Eq. 



194 1, the flow starts at 



k = A in the symmetric phase. In particular, at this 
stage of the flow, k is much larger than the many-body 
scales k-p and T 1 ^ 2 , which thus cannot be resolved. Ad- 
ditionally, bosonic contributions to the right hand side 
of the flow equation vanish, because they are propor- 
tional to the condensate po. Solving the remaining 
flow equations, we then find the dimensionless combi- 
nation ki/k to be rapidly attracted to the (partial) fixed 



point value 6n H1291 . This means that after a few RG 
steps, the Yukawa coupling shows the scaling behavior 
h k ~ y/6n 2 k in the ultraviolet. The approach to the fixed 
point is faster for a larger initial value of h\. It is only a 
partial fixed point, because the scaling solution becomes 
invalid when the scales provided by chemical potential 
and temperature enter the flow. We visualize this behav- 



ior in Fig. 21 



After h 2 k has approached its partial fixed point, the 
value of Sk is attracted to unity, irrespective of the ini- 
tial choice of S A . Therefore, the choice of Z A and A A 
is irrelevant for a sufficiently large h A . This is the case 
for a broad resonance. For definiteness we start with 



4 One of the couplings in the inverses propagator can always be 
absorbed, because the field equation can be premultiplied with an ar- 
bitrary nonzero number. Here we choose A k . 
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the propagator by a set of ordinary differential equa- 
tions. (For such a procedure to work well, the trunca- 
tion of Fk has to be physically well motivated.) The 
corresponding coupled equations are easily solved with 
standard numerical methods - although we are calculat- 
ing the effective action of a strongly interacting quantum 
field theory. We could also expand Uk(p) to a certain 
order in p and then consider only ordinary differential 
equations for the coefficients. However, for the moment 
we keep the full flowing potential, because we will re- 
cover several structures which are already familiar from 
the previous sections on functional methods. 



Figure 21: Flow of the Yukawa coupling = h/A^ 2 for different 
initial values ft/^/A^ 2 = 7, 4, 2 as a function of t = log(i/A). The 
curves are attracted to the scaling solution lit = (67r 2 &)'' 2 shown as 
an orange dashed line. Memory of the initial value is then lost. The 
approach to the corresponding partial fixed point of lr /kis faster for 
a larger initial value. In particular, for h\ = (6ff 2 A)'' 2 = 7.7A 1 ' 2 
we start the flow in the scaling regime. For this plot we have chosen 
T = and (kpay 1 = 0.3, although the situation is generic and also 
valid at nonzero temperature and on the BCS side of the crossover. 
For the chosen parameters, the mass becomes zero for t x -6.5 and 
the flow leaves the scaling regime at this point. 



the fixed point value h A - in the UV and set 

A A = Z A = 1 . Even if we had set the latter two quan- 
tities to zero, the term S = 1 would immediately 
be generated. This provides an RG argument for the 
equivalence of the models with and without dynami- 
cal bosons <p. In particular, after introducing the mass 
2M = 1, we find the boson kinetic term for S — > 1 to 
approach the form 



(197) 



which reveals the boson mass to be twice the fermion 
mass, = 2M. 

Choosing the initial conditions such that we start at 
the fixed point value of hk, we assume the flow of the 
Yukawa coupling to vanish: 



hk — h\ — h . 



(198) 



This is valid in the symmetric phase, but neglects a con- 
tribution proportional to po in the broken phase. The 
latter, however, is expected to be subleading, and this is 
confirmed numerically (35). Note that h\ = h 2 /Ak has a 
scale-dependence due to the running of Ak- 



in Eq. ( 193 1, the inverse boson propagator is assumed 
to be of the simple form iZ^w + ^Akq 2 . Thus, instead of 
resolving the full functional dependence of the propa- 
gator, we only follow the flow of and A*. This corre- 
sponds to replacing the partial differential equation for 



Flow of the effective potential 

Now we are aiming at a computation of the effective 
potential U(p) — Uk=o(p). For this purpose, we start 
from the Wetterich equation 



d k T k [(/,,i/f] = ^Tr 



d k R, 



bosonic contribution 



-Tr 



dkR</,,k 



r (2) 

V i//,k 



(199) 



fermionic contribution 



Here, we have already performed the summation over 
the field indices with a minus sign for the fermions. The 
field content is quite simple - bosons and fermions if/. 
So far, 4>(X) and i//(X) are arbitrary fields parametrizing 
the effective action. We need to keep them since we 
may want to take functional derivatives of Eq. ( 199 1. 



After we have derived all flow equations T k n) of inter- 
est, we insert physical values (possible solutions of the 
field equations) </> = const, and i// = 0, because the 
ground state of the theory will necessarily have a van- 
ishing fermion expectation value and the bosonic field 
will be constant. 

Without loss of generality, we assume a non-negative 
real value <p = yfp for the expectation value of the com- 
plex Bose field. We emphasize that we do not explicitly 
break the symmetry yet, because p could be zero. We 
rather assume that // there is a symmetry breaking, then 
the vacuum expectation value is real. The advantage of 
this is that by decomposing 

<P(X) = 4> + -U<tyi (*) + i<%(*)) , (200) 
V2 

into real fields 6<f>i, we can distinguish the radial mode 
from the massless Goldstone mode via 5<p\ and 6<f>2, re- 
spectively. At the end of the calculation, we determine 
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the ground state value po by minimizing Uk=o(p) with 
respect to p. If po = 0, we arrived in the symmetric 
phase above T c , whereas for po + we have a spon- 
taneous breaking of the U(l)-symmetry and are below 
T c . 

The inverse propagators in the real basis have the 
form lf\Q', Q) = G _1 (Q)tf(Q' + Q), which we encoun- 
tered above and are given by 



\A k f + U' k (p) + 2pU' k '(p) -Z k cj n 



Z k co„ 



\A k c? + U' k {p) ) 
(201) 

(202) 



-h(p*s ia>„ - (q - p) 
ito n + q 2 - p h(ps 



Compare this to the similar expressions in Eqs. (|C.8 



and ( 137 i. Here, e is the antisymmetric tensor in two di- 



mensions (cf. Eq. ( 125 1), and the unit matrix in two di- 



RkAQ) = 





) 



o 







mensions on the off-diagonal of Eq. (202 1 is suppressed. 
The matrices correspond to the orderings {4>\,<t>i) and 
(tf/\,i[f2, iAp 1A2 ) °f tne f" e ld variables, respectively. 



The flow equation ( 199 1 contains the regularized ex- 
pression (Tf } +R k )~ 1 . We use Litim cutoffs here. With 
this choice, we effectively replace q 2 /2 — > k 2 for the 
bosons. The fermions are regularized around the Fermi 
surface such that +\q 2 -p\ is replaced by +k 2 . Explicitly, 
we have 



'4x4 

(203) 



with regulator functions 

r = A k (k 2 - f/2)6(k 2 - <f/2) , (204) 
f> = (sgn(^ 2 - p)k 2 -f+ p)6(k 2 - \f - p\) . (205) 

Due to overall theta-functions arising from d k R k , the q- 
integration is restricted to regions where q 2 is replaced 
by k 2 such that the integration over spatial momenta be- 
comes trivial. 

Then, the full flow equation for the effective potential 
is given by 

d k U k (p) = \ T Y;f ^^{G^ k (Q)d k R k 4Q)) 

-\ T ^J ^{ G ^Q)dkRkm))- 



(206) 



Besides the simple integration over spatial momenta, 
the Matsubara summations can be carried out explic- 
itly, because the cl>„ are not cut off at all. The result- 
ing flow equation is a partial differential equation for 
U k (p) = U(k,p). We find 



d k U k {p) = 



y/2k b 
3n 2 Si 



1 + W\ 

1 + W2 



I +W2 



1 + W 



1 J 



'1 1 

X I - + - 

[ 2 e V(i+wi)(i+»'2)A: 2 /Si.r _ i 



2k 3 



3n 2 ^|T+ 



W>3 



1 



1 



/T+mk 2 /T + i 



(207) 



where m = U' k (p)/k 2 , w 2 = (U' k (p)+2pU' k '(p))/k 2 , w 3 = 
h 2 p/k 4 and t{p) = (p + k 2 ) 3/2 6(p + k 2 ) -(p- k 2 f l2 6{p - 
k 2 ). At zero temperature, we set the Bose and Fermi 
distribution functions to zero. 



The flow of the effective potential in Eq. ( 207 1 ex- 
plicitly depends on S k — Z k /A k . Therefore, in order to 
close the expression, we also need flow equations for the 
wave function renormalization Z k and the gradient coef- 
ficient A k . These are obtained from the flow of the in- 
verse boson propagator G^ 1 via suitable projection pre- 
scriptions. For example, given the flow of the boson 
two-point function if^Q', Q) = 5{Q' + fi)G^(Q), we 
obtain the flow equation for Z k according to 



3 



^(G-) I2(W ,0) 



(208) 



w=0 



The flow equation for A k is derived analogously. We 
can then again evaluate the Matsubara summations and 
^-integrations. These calculations are a little intricate, 
but standard and straightforward in principle. 



Three building blocks 



Approximative solutions to the flow equation ( 207 1 for 



the effective potential have to fulfill three key require- 
ments. These are independent of the particular trunca- 
tion and summarize our earlier considerations on gen- 
eral properties of the effective potential. The compu- 
tation of U k (p) has to account for the following three 
building blocks: 

1) Phase diagram po(p,T) 

The phase structure of the system under consider- 
ation is found from the minimum of the effective 
potential U k= o(p). Therefore, we have to allow for 
spontaneous symmetry breaking in the construc- 
tion of the flow equations. 
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2) Equation of state n(p, T) 

We obtain the equation of state from the full effec- 
tive potential via 



PQx,T) = -U k=0 (ji, T,po). 



(209) 



In most cases, we will be satisfied by having an 
expression for the density n(p, T) = (dP/dp)r- In 
particular, by inverting the relation n = n(p, T) for 
p = p(n, T), we can eliminate the chemical poten- 
tial for the density. 

3) Vacuum flow of the couplings 

As discussed earlier in this section, the micro- 
scopic couplings appearing in the action S are not 
the ones we will measure in experiment. For this 
reason, we equipped them with a subscript A. In 
the BCS-BEC crossover, we obtain the equation of 
state and the phase diagram as a function of va and 
h. In order to find out how these two parameters 
are connected to the scattering length a, we have 
to solve the vacuum flow equations. After that, we 
can express our observables as a function of k F a, 
which can then be compared to experiment. 

We use the first and second points as guidelines for the 
construction of a suitable truncation of Ukip). 

In order to include the effect of spontaneous symme- 
try breaking, we remark that 







dU, 



(210) 



has the three types of solutions discussed in Sec. 3.2 



however, at the scale k. The observable phase diagram 
is found from the corresponding value of po^ m the limit 
k — > 0. Since we expect a second order phase transition 
here, an expansion 



U k (p) = m^ip - p ,i) + -tt(p ~ Pojc) > 



(211) 



is sufficient to get the phase structure, cf. Eq. ( 148 1. In- 



serting this ansatz into Eq. ( 207 1 we get a set of coupled 



ordinary differential equations for m 2 k , u^ k and po,/fc- 
The inclusion of higher terms yields further quantitative 
improvement, but can be neglected in a first attempt. 
By construction, the effective potential has a minimum 
at po k . We divide the flow into two regimes: 



(i) symmetric regime: m 2 . k > 0, po^ - 

(ii) broken regime: m 2 . k = 0, po,jt + 0. 



0, 



There is no condensate for k — A, because the micro- 
scopic potential is given by U\(p) = (v A - 2p)p. We 



start in the symmetric regime and follow the flow of 
m^ k . If this quantity hits zero, we switch to the flow 
equations for the broken regime. The mass is then fixed 
to zero, tn^ k — 0, and the flow of the minimum pQ k is 
found from 

= d k (U' k (p 0M )) = d k U' k \ Po (p 0ik ) + U k (p 0tk )d k p Qtk , 

(212) 

to be given by 



dkPo,k = - 



dkU' k \ Po (po, k ) 



u 4>,k 



(213) 



The quantity in the numerator can be obtained from a 
derivative of the flow equation for U k . Note that we 
follow the physics from small to large scales by inves- 
tigating the scale-dependent position of the minimum 
po,t. If we arrive at po > for k — 0, we are in the 
phase of spontaneous symmetry breaking. However, it 
may even occur that in the broken regime of the flow 
the minimum becomes zero and we have to switch flow 
equations again. Only the (k - 0)-value is of physical 
relevance. We will discuss the related issue of precon- 
densation below. 

The expansion in Eq. (211 1 can be extended to in- 



clude pressure and density. For this purpose we write 
Uk(p) = -Pk-n k -{p - p) 

+ m\ k {p - po,*) + -7jr(p- Po,kf , (214) 

where p' is an artificially introduced off-shell chemical 
potential which replaces p in all previous calculations. 
It is used to generate the flow equation of n k by taking 
a p'-derivative of Eq. (207 i for the effective potential. 



Once we have the corresponding flow equation, we set 
p' = p. In fact, this is the same procedure as for obtain- 
ing slope and curvature of the effective potential {m 2 ^ k 
and u$jc) by taking derivatives with respect to p and set- 
ting p = pojc afterwards. Again, we may include higher 
orders in p - po,* and p' - p to obtain quantitative im- 
provements. 



UV renormalization and fine-tuning of the mass 

The solution of the flow equations derived above results 
in «(p, T) and po(p, T) as functions of Va and h. For a 
comparison with experiment we have to relate this set 
of parameters to the scattering length a. Another neces- 
sity for this UV renormalization is that the values of the 
microscopic action are strongly cutoff dependent. In- 
deed, we have a dependence m 2 k ~k for the running 
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mass term for k — > oo and thus, if we slightly change A, 
we drastically change rrp\ A . Since our UV cutoff is only 
roughly determined to be A » C<jw> tne precise value 
of m 2 A cannot be of any physical relevance. 

Let us briefly recapitulate the general reasoning de- 
scribed earlier and then apply it to the present case. The 
input for the flow equations is given by yU, T and the 
microscopic couplings. As k is lowered from A to 0, 
we include quantum and thermal fluctuations. For high 
momentum scales k, only quantum fluctuations are im- 
portant and a few-particle system will behave similar to 
a many-body system. This situations is altered at the 
scales kp and A T l which are present at nonzero densi- 
ties and temperatures. When k reaches these thresholds, 
many-body effects set in and the flow of observables is 
different from the pure vacuum flow. The FRG idea of 
UV renormalization is to follow the flow twice: First 
we solve the flow equations for /i = yu v and T — 
such that the many-body scales are never reached. This 
tells us what scattering length a we use. Then, in a sec- 
ond attempt, we use the same initial conditions but vary 
/J. = fi v + fi m b and T > to get the functions nip., T, a) 
and po([i, T, a). 

In the ultraviolet, the effective potential is given by 



U A (p) = Oa - 2/i)p . 



(215) 



For conceptual clarity we split up the chemical potential 
into a vacuum and many-body part according to 



(216) 



We have /i v < and // m b > 0. This decomposition is, 
of course, artificial, because the chemical potential is 
simply a parameter to change the density. However, we 
will see that // v is related to the binding energy of the 
molecules and only appears with nonzero value on the 
BEC-side, whereas /z m (, determines the density in the 
crossover and provides the scale kp. 

The chemical potential appears in the inverse fermion 
propagator = (d T - V 2 - //), which is shown below 
not to be renormalized in vacuum. A vanishing Fermi 
distribution function n$ = J^P^ ] (Q) thus requires p. < 
0. For the vacuum problem we set p m \, = 0. Moreover, 
we have T = 0, and since the flow will always be in the 
symmetric phase, we effectively have p — 0. The latter 
fact is related to the absence of condensation. 

The flow of the mass m 2 k = U' k (po,k) is found from 



Eq. (207 1 to be 



= gj^s (k 2 + ^f l2 0(k 2 + Hv) ■ (217) 



The right hand side is solely due to the fermionic con- 
tribution. This equation is solved by 



with 



£_ 

(VT 2 



Va^(i-^) 

-f v^arctan( A ^ v ) 



0"v = 0), 

/i v < 0, 

k < Kl 



(218) 



(219) 



Since m?(& = 0) and \fi Y \ act as gaps for the excitation of 
bosons and fermions, respectively, we can split up the 
physical vacuum at k = into three sectors according to 



(i) atom sector (a < 0): rn 2 ^ k=(l > 0, fi v — 0, 
<M=0 



(ii) resonance (a 1 = 0): ni 2 , =0 = fi v — 0, 



(iii) dimer sector (a > 0): m 2 ^^ — 0, /v v < 0. 

The three sectors are distinguished by the corresponding 
type of particles which interact in the vacuum and, in ad- 
dition, they are related to a certain scattering length. The 
latter statement is proven below. Note the formal anal- 
ogy to the classification of the thermodynamic phase di- 
agram, cf. the discussion below Eq. ( |68j ). The positive 
mass term (i) corresponds to the symmetric phase, (ii) 
relates to the critical point. In (iii), a scale fj. v is gener- 
ated, similarly to po in the many-body problem. Despite 
the fact that no symmetry is broken, the onset of the 
molecular bound state shares features of a phase transi- 
tion in vacuum. Indeed, there is a (spatial) continuum 
of degrees of freedom. 

The condition on /i v in (i) - (iii) is easily satisfied by 
setting the chemical potential in the microscopic action 
to a certain value. For the condition on m 2 k=0 , we use 



our solution (218 i of the flow equation to find a suitable 
choice of va. Obviously, given yu v = 0, we will have 
m 2 k=0 = precisely if va = ^ A. In the same fashion, 

> 0, i.e. the atom 



yu v = and v A > g^A yields m 2 ^ 
sector of individual fermions. Note that indeed m 2 



4>,A 



A strongly depends on the cutoff. If // v < 0, we take the 



second line of Eq. (|218| and solve m 2 k=0 = for va. 



This can easily be done numerically. With this choice 
of va we are in the molecule phase. 

If we now switch on fx m b, we could already tell from 
the choice of va and p. w whether we are on the BCS- 
side, on resonance, or on the BEC-side. However, we 
still need to relate Eq. ( |219| l to the precise value of 
the observable fermion scattering length a. From Eq. 
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( 190 1 we know that /l^, e ff,A = _ ^ a /va is valid in the UV. 



A similar equation holds for the renormalized coupling 
/V.eff = 87ra - It is given by /l^ ;e ff = -h 2 /P,p(cj,q = 0,fi v ), 
where is the boson propagator analytically contin- 
ued to real time frequencies u> = -\io„. The expression 
has to be evaluated for the on-shell condition of fermion 
scattering. We do not dive into the details here. One can 
show that in our case we always have 

h 2 

8na = A^eB = 5 ; — . (220) 



<p,k 



=0 (jUv = 0) 



From Eq. ( 220 1 we obtain the following dictionary to 
translate (h^, v\) — > a. The fermion scattering length a 
is always given by 



a(B) 



h 1 



v(B) = v A 



r A. (221) 



Herein, h 
Eq 



8ttv(B)' ' v ' 6n 2 
h\ and v(B) is the physical detuning from 



( 183 i. It corresponds to the energy difference be- 
tween the closed channel bound state and the open chan- 
nel scattering threshold. We write the explicit depen- 
dence of a and v on the magnetic field B to emphasize 
that this parameter can be tuned experimentally. If we 
choose v(B) > 0, the expression in the denominator of 
Eq. ( |221| i is positive and we are in the atom sector, or 
- at finite density - on the BCS side. We then have 
M = i-i m b- If v(B) < 0, the scattering length gets positive 
and we are on the BEC side. Part of the chemical po- 
tential yu = ju v + yu m b is then determined by the condition 
that molecules in the vacuum are gapless, propagating 
degrees of freedom. This contribution p v is found as a 
solution of 

h 2 







l tf>,k=0 



v A - 2fi v 



(^2+^(1 




/At \ 

2 A 2 ^ 



yU v arctan( 



VA 2 + fly 



(222) 



for given va. 

We observe that only on the BEC-side the chemical 
potential partially contributes to the vacuum flow. The 
physical interpretation of p v is very intuitive. It is half 
the binding energy of a dimer, 

£b 
2 

The binding energy is a negative quantity and thus the 
total chemical potential fi = ^ + fx m b might be negative. 



/A 



(BEC side) . 



(223) 



In fact, in the limit 
< |222} 



we find from Eq. 



Ma 2 



(BEC-side) . 



(224) 



This is the well known universal relation for a molec- 
ular bound state from quantum mechanics. If a gets 
small, sb becomes large and negative. Note that if we 
do not neglect higher orders in -J—p v /A — > 0, we get 



deviations from Eq. ( 224 1, which are related to micro- 



scopic details of the interaction potential and thus non- 
universal. 

We briefly comment here on the dimer-dimer scat- 
tering length. The system of two-component fermions 
can be described as a gas of bosons on the BEC side. 
This Bose gas has a scattering length between its con- 
stituents. The exact value is known from the quantum 
mechanical calculation [133] to be 



a 



0.6. 



(225) 



A mean field calculation omitting bosonic fluctuations 
gives a^/a = 2. Including fluctuations of the molecules 
yields a lower value. Within our truncation we find 
as I a — 0.72. 



Crossover at zero temperature 

The analysis of the vacuum problem allows to interpret 
the results of the flow equation for given initial con- 
ditions in terms of the scattering length a. Tuning the 
chemical potential fj., we can now create a nonzero den- 
sity. Many-body effects become important once the flow 
parameter k reaches k m <o. The earlier considerations on 
the initial conditions in the UV regime remain valid. 
Making contact to the discussion at the beginning of this 
section, we take the limit T — > 0, but now at fixed phase 
space density a>. 

From our building blocks stated above we conclude 
that the solution of the flow equation for the effective 
potential consists in essence in the determination of 



n(p) (equation of state) , 
A(p) = h 2 p (p) (gap) . 



(226) 



We introduced here the gap, which has precisely the 
same meaning as in the discussion of the BCS theory. 

We employ the truncations of the effective action and 
effective potential given in Eqs. ( 193 i and (214 1 for 
A<mb = fi - /M/ > 0. The resulting flow equations for 
Z t , A 



h m 



<l>,k " 



Po,yt> u <p,h Pk ar, d n k can easily be solved 
numerically. The general picture is that the fermions 
break the symmetry in the early stages of the flow en- 
tering the many-body regime k w k m b, i.e. the fermionic 
contributions on the right hand side of the flow equation 



d207b are responsible for a decrease of the mass m 2 



<p.k ' 
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(tea) -1 




(fc F a) 



Figure 22: We plot the chemical potential (equation of state) as a 
function of (kfCiY 1 for T = and T = T c , respectively. The Fermi 
momentum kp = (3;r 2 «)'' 3 corresponds to the density n. The binding 
energy of the molecular bound state on the BEC side for positive a 
has been subtracted. The po sitive quantity p - eb/2 corresponds to 
p m h in the notation of Eq. |216| . For (kpay 1 jS _ 2, we observe 
£b = and p = s^(n). This is equivalent to the equation of state for 
ideal fermions. We find an explanation in the exponentially small gap 
A > 0, which opens up in the dispersion E q (cf. Fig. |15| ) Due to its 
tiny value, it does not influence the equation of state significantly. For 
zero temperature, we can read off the Bertsch parameter f = p/£f at 
unitarity. 



which ultimately hits zero at k > 0. At this stage of the 
flow, we switch flow equations and follow the flow of 
po t k instead of m 2 k . The bosons were quite uninvolved 
in the symmetric regime flow. Now, the bosons become 
very active and try to restore the symmetry, and thus 
slow down the running of po,/t- We find that for k —> 
the values of all quantities saturate and we can read off 
the result at k = 0. (These observations also hold at 
nonzero temperature. But, depending on yu and T, the 
mass term may not be zero for k — and we end up in 
the symmetric phase.) 

The solution to the flow equation for the equation 
of state and the gap parameter as a function of the 
crossover parameter k^a are shown in Figs. (22i and 



( 23 1. These plots already yield a full qualitative under- 
standing of the BCS-BEC crossover at zero tempera- 
ture. If we wish to extend our truncation, we simply 
have to derive more flow equations and integrate them 
numerically. 

In the BCS limit (a < 0, \k?a\ «: 1), we find /i -> 
sp{n). Thus we arrive at a weakly interacting Fermi gas 
with a clearly expressed Fermi surface. We find good 
agreement with the perturbative one-loop result from 
BCS theory 



1 



Figure 23: The gap parameter A 2 = /i 2 po, shown here for 7 = 0, con- 
stitutes the order parameter of the superfluid phase transition in the 
BCS-BEC crossover. It is nonzero for all values of (kpa) ~' at zero 
temperature. We find excellent agreement of our truncation with the 
BCS result from Eq. ^228) . It is well-known that there is a correc- 
tion to the BCS result due to screening effects, which has first been 
calculated by Gorkov and Melik-Barkhudarov. We show the corre- 
sponding behavior of the gap Aoorkov = (2/e) 7 ' 3 e _ ' r ' 2 '' rFa ', which can 
be captured with the FRG from a more elaborate truncation including 
particle-hole fluctuations. 



p/e F 




Figure 24: The fermion chemical potential p includes the binding en- 
ergy of the molecule on the BEC side of the crossover. It eventually 
becomes negative for (fed) » 0.5 and the behavior of the system is 
then dominated by few-body physics. The chemical potential in vac- 
uum (n = 0), which is equivalent to half the binding energy of the 
bound state, is shown in red. We observe that indeed both functions 
merge for ffepa) -1 > 1. The smooth behavior at nonzero density ter- 
minates in a sharp second order phase transition in vacuum, which is 
related to the formation of a molecular bound state. 



150 1. For small 



with E q = y/isq - yu) 2 + A 2 , see Eq. 
values of \k F a\, the gap parameter agrees with the BCS 
result 

AfiCS _ ^_ e -n/(2\k f a\) _ ^28) 



•)3 fR.IT + 1 



(227) 



Since only fermions around the Fermi surface contribute 
to the pairing, condensation is weakly expressed. We 

find strong deviations from the BCS result for k^a 1 . 

In the BEC limit (a > Q,kpa <sc 1) we find ju/sf — > 
-oo. Comparing this result to the solution of the vacuum 
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limit, we observe that this regime is largely determined 
by few-body physics. We have 



p(n) = SB/2+fi mb (ri) 



1 



(229) 



The chemical potential approaches half the binding en- 
ergy in this limit and the many-body scale £f(«) £b 
drops out. In particular, we have a zero crossing of the 
chemical potential at (k F ay l w 0.5. This is demon- 



strated in Fig. 24 



The function p = p(kpa) behaves smoothly when go- 
ing from the BCS to the BEC side, which is the manifes- 
tation of a crossover. In contrast, the vacuum chemical 
potential p v = y 6(aA) has a discontinuous derivative at 
a' 1 = 0. 

In the unitary limit cC x = 0, the ratio p/sp is a univer- 
sal number, which is also called the Bertsch parameter 

We find here £ = 0.55. The experimental value is 
given by £ = 0.376(5) P51 . which shows that our sim- 
ple truncation can capture the qualitative effects but fails 



for quantitative precision. In Fig. 22 we display in ad- 
dition to the zero temperature case the behavior of p m \, 
at T — T c and find a larger value of p/sp at unitarity. 
This is also found in experiments [45] but, again, the 
chemical potential obtained here exceeds the measured 
value. 



Finite temperature phase diagram 

It is a particular strength of the FRG that calculations 
at nonzero temperature are conceptually and technically 
as straightforward as the corresponding computations at 
zero temperature. We show the finite temperature phase 
diagram of the crossover in Figs. ( |25] l and ( |26] > in terms 
of (p,T,a) and (n,T,a), respectively. The curves are 
obtained with the FRG from the basic truncation given 
in Eq. 



(193i 



The plots show regions, where the sys- 
tem is either in the normal (symmetric) phase or in 
the superfluid phase with spontaneous breaking of the 
global U(l)-symmetry. The superfluid phase transition 
is found to be of second order throughout the whole 
crossover. This justifies our truncation of the effective 



potential in Eq. (21 1 1 a posteriori. The temperature de- 



pendence of the gap is shown in Fig. 27 

Note that the phase diagram is equivalent to a plot of 
the critical temperature T c . The value of T c for strong 
coupling is particularly interesting. Within our trunca- 
tion we find T c /p = 0.44 and T c /Tp = 0.28 at unitar- 
ity. This has to be compared to the experimental values 
TJp = 0.32(3) El and TJT F = 0.167(13) 03. As 
for the Bertsch parameter, quantitative precision is not 
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Figure 25: Phase diagram of the three-dimensional BCS-BEC 
crossover in the plane spanned by temperature, chemical potential 
and scattering length. Here, /j ra (, = fi - sb /2 is positive because the 
binding energy has already been subtracted. Thus, y// m Bfl < and 
■Jlimh > correspond to the BCS and BEC sides of the crossover, 
respectively. The critical temperature separates the superfluid from 
the normal phase. In addition, there is a precondensation temperature 
below which a nonvanishing field expectation value po,* appears at in- 
termediate stages k of the RG flow. In the BCS limit, T c and T„ c are 
practically indistinguishable. 
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Figure 26: Combing the equation of state at T > and the phase 
diagram from Fig. [53] we obtain the phase diagram in terms of T c /Tj: 
and kpa. We indicate the limiting cases on the BCS and BEC sides 
by dashed lines. The BCS formula \146) is found here to match for 
(k^a) < — 1. The condensation temperature of ideal bosonic dimers 
Tc.id/TV = 0.218 is approached on the BEC side with a correction 
proportional to the diluteness parameter k^a, see Eq. ^231). 



found within this basic truncation. Part of the error in 
T c /Tp consists in determining 7f(h) = (37r 2 «) 2 ^ 3 at uni- 
tarity, i.e. the equation of state. 

We find remarkable agreement of the value for the 
critical temperature in the limiting cases of small cou- 
pling. On the BCS side, formula ( 146 1 is found to be 
valid. In particular, the exponential vanishing of T c 
with the correct exponent can be verified in a logarith- 
mic plot, which is not shown here. For small positive 
a = ay, we expect the system to be described by a gas 
of weakly interacting bosons with scattering length a^. 
Due to interactions effects, the critical temperature de- 
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Figure 27: Temperature dependence of the gap A = h po, which con- 
stitutes the order parameter for the superfluid phase transition in the 
BCS-BEC crossover. The dotted, solid and dashed lines correspond to 
BCS side, unitarity and BEC side, respectively. We rescaled the gap 
by its zero temperature value to obtain better visibility of the expo- 
nentially small gap in the BCS phase. Note that A(T) vanishes contin- 
uously for T — > r c , thus revealing the phase transition to be of second 
order. It is a particular strength of the FRG that it properly accounts 
for the fluctuations at criticality, which diminish the order parame- 
ter and eventually drive it to zero. Neglecting these fluctuations, the 
phase transition may (wrongly) appear to be of first order. 



viates from the ideal gas result T^d/Tp = 0.218 (cf. Eq. 



( 152 1) according to 



TV - T, 



1/3 



(for small a > 0) (230) 



with a dimensionless constant k. Rewriting this expres- 
sion as 



Tc - T c 



T cM (6n 2 y^ a 
we can apply our result a^/a — 0.72 to find k — 1.7. 



(231) 



In Fig. 25 we show the precondensation tempera- 



ture r pc . For T < Tpc, there is a momentum scale 
^sbb where a nonzero field expectation value po,t ap- 
pears during the flow. The superfluid phase corresponds 
to those temperatures, where po,jfc=o = Po is nonzero for 
k — and thus constitutes the order parameter of the 
phase transition. For intermediate temperatures T such 
that T c < T < 7p C , the field expectation value does not 
survive in the infrared and we arrive in the symmetric 
phase. This precondensation region can be viewed as a 
state of the system where we have correlated domains 
of size ~ k~ d . We visualize the flow of po,k in Fig. 28 



Let us pause here for a moment to discuss the rela- 
tion of the truncation presented in this section to other 
approximations. For this purpose, we start from a 
more basic ansatz and then successively build in ad- 
ditional effects until arriving at the present truncation 
scheme. Quite remarkably, the qualitative physics of the 
crossover problem at zero temperature can be described 



Figure 28: The precondensation temperature T pc is defined as the 
highest temperature where a nonzero value of po,t appears at k > 0. 
This minimum of Uk is not necessarily nonzero for k = 0. For this rea- 
son, there is a region in the phase diagram where precondensation oc- 
curs, but the system is not yet superfluid. The distinction between T c 
and T pc is illustrated in this plot. The flow parameter is parametrized 
according to t = log(fc/A). We have po,j = for T > T„ c . In this plot, 



we keep y// ra (,a fixed, which corresponds to (kva) 
above (below) T c , respectively. 



0.74 (0.72) 



in terms of just a BCS gap equation ( 145 I, together with 



a self-consistent treatment of the equation of state ( 149 1 
with fermionic contribution alone. While in the BCS 
regime the equation of state is well approximated (up 
to exponentially small corrections) by the Fermi sphere 
contribution (i.e. n - (2Mp) 3/2 /(37r 2 ) with = e F ), 
ramping the inverse scattering length to large positive 
values leads to n — > 20*0, which describes a (renormal- 
ized) condensate without depletion. This is the original 
approximation of Eagels |134| and Leggett [135|. The 
BEC regime here is described in terms of just a con- 
densate, but clearly there are no propagating bosonic 
degrees of freedom. The next stage of approximation 
improves on this point, and includes fermionic fluctua- 
tions which build up a bosonic propagator 1 136.1 11371 . 
This then allows to qualitatively describe the crossover 
at finite temperatures, including the phase border to 
the symmetric phase. This approximation yields a full, 
qualitatively correct finite temperature phase diagram, 
which often is referred to as the extended mean field 
or Nozieres-Schmitt-Rink approximation. In the FRG 
language, it corresponds to keeping only the fermionic 
contributions to the flow equations of Uk, Z^, and A^, 
respectively, and no bosonic feedback. This allows to 
describe an effective theory of pointlike bosons in the 
BEC regime, which is characterized by a nonrelativistic 
mass 2M and an effective scattering length = 2a. 

A consistent generalization of the extended mean 
field theory by semi-analytical means is not straight- 
forward. It has been addressed in the frame of func- 
tional field-theoretical techniques, in particular via e- 
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expansion lfT38Ul43l . 1 /^-expansion Q301 [144|, t- 
matrix approaches 111451 41501. Dyson-Schwinger equa- 
tions lfTl6l [T27l [1511 . 2-Particle Irreducible methods 
1 152 1, and renormalization group flow equations |129 
I153H1571 . An immediate challenge can be inferred 
from the fact that the inclusion of fluctuations related 
to the bosonic sector leads to infrared divergent inte- 
grals already in perturbation theory. We have met an 
example of this situation in our analysis of extensions 
of Bogoliubov's theory for weakly interacting bosons 
at zero temperature; at finite temperature, those diver- 
gences become even more serve. In addition, a notori- 
ous problem known for bosonic theories is the order of 
the finite temperature phase transition, which typically 
is wrongly found to be of first order. 

Already the simple truncation scheme Ml 541 advo- 
cated above, provides the means to consistently and sys- 
tematically deal with these problems. The flow equa- 
tion is both UV and IR finite and thus it is not plagued 
by the problems mentioned above. Physically, this al- 
lows to include the order parameter fluctuations in the 
BCS regime, and the thermal and quantum fluctuations 
of the effective molecular bound state degree of free- 
dom in a single framework. We have described a few 
important improvements resulting from this approach 
on various scales: On short distances, bosonic vacuum 
fluctuations renormalize the effective bosonic scattering 
length already close to the exact value obtained from 
the solution of the four-body Schrodinger equation. On 
thermodynamic scales, we observe a shift in the criti- 
cal temperature, predicted for pointlike bosons from ef- 
fective field theory and Monte Carlo simulations. At 
long distances, we correctly capture the second order 
phase transition throughout the whole crossover, char- 
acterized by an anomalous dimension close to the best 
available estimate from the e-expansion and numerical 
approaches. This discussion is summarized in Tab. [3] 



Efimov physics at resonance 



Physical scale 


MF 


FRG 


Other 


Microphysics, — 


2 


0.72 


0.6 (Ref. 11331) 


Thermodynamics, k 





1.7 


1.3 (1158H160I) 


Critical behavior, 77 





0.05 


0.038 (Ref. 111611) 



Table 3: Quantitative precision on all scales: The BCS-BEC crossover 
shows a separation of scales, which can be benchmarked by several 
key observables. We display here a selection of representatives for 
each sector and its value in extended mean field theory (MF), from 
our FRG truncation introduced above, and other methods. The mi- 
crophysics are sensitive to the scattering of composite bosons with 
effective molecular scattering length A nonperturbative thermo- 
dynamic effect on the BEC side consists in the shift of the critical 
temperature relative to the ideal gas value, &T c BEc/T c jj = Ka$n^ 3 , 
The critical exponent n known from the theory of critical phenomena 
has to be compared to the three-dimensional 0(2) universality class, 
because this is the symmetry of the complex order parameter (pQ. 



These procedures remain tractable since massive dia- 
grammatic simplifications occur in the vacuum limit 
( |181[ ), leading to a closed hierarchy of flow equations of 
the A^-body sector. (In other words, not unexpectedly, 
the nonrelativistic A^-body problem can be solved with- 
out knowledge of the (N + l)-body problem.) Here we 
give a flavor of this physics only, discussing the Efimov 
effect of three resonantly interacting particles, which is 
explained by a limit cycle behavior of the RG flow. For 
some further aspects of vacuum physics, we refer to 



Ap- 



421 



pendix E In particular, we point to the reviews l37l 
addressing Efimov physics from a field theoretical per- 
spective. 

Efimov found the effect named after him from a con- 
sideration of the Schrodinger equation for three reso- 
nantly interacting identical bosons. He showed that this 
problem can be mapped to the scattering in an inverse 
square potential at short distances 1166111671 . The latter 
potential has a discrete spectrum of bound states which 
form a geometric series, i.e. two neighboring bound 
states E„ and E„ + \ satisfy 



At resonance, the gas of two-component fermions is 
dense and strongly interacting. In particular, three-body 
processes become relevant. Here, we would like to 
give a short excursion which demonstrates that the few- 
body sector in the strongly interacting regime aT x — > 
holds interesting physics in its own right, the quali- 
tatively most interesting part occurring for bosons or 
three-component fermions. The corresponding physics 
of few-particle systems in vacuum can be investigated 
within an FRG framework using vertex expansions 
111621 11631 or refermionization techniques [164 



f 



(232) 



where sq 1.00624 is know from quantum mechan- 
ics. This scaling of the infinite tower of bound states 
is known as the Efimov effect. With the help of Fesh- 
bach resonances, this effect can be studied using ultra- 
cold quantum gases [ 1681 41811 . 

In order to describe this problem in a field theoreti- 
cal FRG framework, we need to include an atom-dimer 
interaction. The full solution for the momentum depen- 
dent atom-dimer vertex can be found within a vertex ex- 
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pansion of the FRG 11631 . The dimer can consist of 
either two fermions or two bosons, see the examples be- 
low. Using kinematic simplifications and projection to 
zero angular momentum partial waves, the flow equa- 
tion for the atom-dimer vertex 73^ can be reduced to a 
quadratic matrix differential equation 11621 11631 . We 
will not present details on these approaches here, but 
rather explain the basic phenomenon underlying the Efi- 
mov effect: a renormalization group limit cycle. It is a 
strength of the FRG that this qualitative effect can be 
identified in a minimalistic truncation in terms of a sin- 
gle atom-dimer interaction parameter [ 163 1. 

To this end we consider the following types of 
three-particle systems: (i) identical bosons, (ii) two- 
component fermions (i.e. two hyperfine states), and (iii) 
three-component fermions (i.e. three hyperfine states). 
Systems (i) and (iii) are assumed to have a SU(2) and 
SU(3) spin symmetry, respectively. For the latter case, 
this is not generic, because the resonances of the three 
mutual scattering lengths a^iB), 013(B) and 023(B), 
as a function of the magnetic field B, are in general 
not identical - there is no fundamental SU(3) symme- 
try in cold atoms. However, fine-tuning of parame- 
ters can realize such a situation to sufficient accuracy 



II 1741 1178L For theoretical work on this problem, see 
Refs. lfT6^fT65lfT82UT84l . 

We can now give a simple RG argument to identify 
which of the above three-particle systems allows for the 
Efimov effect |163|. We approximate the atom-dimer- 
vertex matrix by a single entry independent of momen- 
tum, i.e. a single coupling constant A^. The corre- 
sponding flow equation for the dimensionless^Jcoupling 
A3 k — M kk 2 has the structure 



d t Xz = aA\ + BA3 + y, 



(233) 



where the constants depend on the choice (i) - (iii). The 
solution to this equation is given by 



hit) 



tanh( VZ>r/2) (D > 0) 
tan( VfDff /2) (D < 0) 



(234) 



where D = B 2 -4ay is the discriminant of the beta func- 



tion in Eq. ( 233 1. The first case (D > 0) is found for sys- 



tem (ii). The coupling A3 then reaches an infrared fixed 
point. Intuitively, this can be understood from Pauli's 
principle, preventing two spin-1/2 fermions from com- 
ing too close to each other. However, for systems (i) 



5 The canonical power counting at the unitary point deviates from 
the canonical power counting for finite a , as can be inferred from 
the exact solution of the two-body problem, see e.g. 1421 . 



and (iii) it turns out that D < 0, and thus the solution 
is periodic in RG-time t = log(fe/A). We say that the 
flow approaches an infrared limit cycle. The periodicity 
of the solution with time is 



T = 2nl VlDi. 



(235) 



For each divergence of the three-body scattering ampli- 
tude at zero frequency and momentum approximated by 
A-sfr, a new Efimov bound state is expressed. This gives 
rise to an infinite number of exponentially spaced Efi- 
mov bound states, i.e. the binding energy of successive 
bound states vanishes exponentially fast. This consider- 
ation is exact at the resonance. Away from it, due to the 
scaling violations caused by cr l , only a finite number of 
Efimov states exists. Experiments with ultracold atoms 
have resolved the lowest (most deeply bound) Efimov 
states. 

In order to compute sq, we perform a scale identifi- 
cation, which relies on the fact that the flowing action 
can be seen as a theory at momentum scale k. However, 
the subtle association of k to physical scales like exter- 
nal momenta, bound state energies etc. is not unique. 
Here, we observe that from dimensional arguments we 
have E ~ k 2 . The prefactor drops out when calculating 
the ratio of two energies. We thus associate two neigh- 
boring dimer bound state energies E„ and E„ +1 with the 
corresponding scales k„ and k n+ \ in the limit cycle. We 
have t n +\ — t„ — T and accordingly 



k 2 

%+l _ -2T 

kl 



(236) 



We compare this to Eq. (232 1 and deduce Sq = n/T 



The Efimov parameter is related to the periodicity of 
the limit cycle. 

From this simple analysis including only the single 
parameter A3, one obtains Sq — 1.393. This is quantita- 
tively not yet very accurate, but we inferred the relevant 
qualitative physics. The FRG approach can be com- 
pleted by solving the flow equation for the momentum 
dependent vertex 73 ^ numerically. The result sq - 1 .0 is 
compatible with the value from other approaches [42]. 



Open challenges and further reading 

As discussed above, the truncation for the effective ac- 
tion given in Eq. ( 193 1 provides a consistent pic- 



ture of the whole crossover, interpolating smoothly be- 
tween the limiting cases of BCS superfluidity of atoms 
and Bose condensation of dimers, and in particular pro- 
viding a consistent treatment of the bosonic sector, with 
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Figure 29: The phase diagram of the BCS-BEC crossover can be di- 
vided into regions where different momentum scales are dominant. 
Including effects beyond mean field theory provides a challenge in 
each individual sector. In addition, a unifying description of the whole 
crossover has to match the correct limits. The zero crossing of the 
fermion chemical potential at ~ 0.5 separates a region where 

many-body effects are important from a region where the main fea- 
tures of the system can be captured by an effective theory of dimers. 
In the latter, we have |eg| a \ 3> k^, T such that microscopic scales 
are relevant. The situation is reversed to the left of the zero cross- 
ing. In the region around the critical line, long distance physics and 
infrared fluctuations become important. 



a number of quantitative improvements. Starting from 
this, open challenges remain. Some of them have been 
successfully addressed and others remain for the future. 
For a more detailed discussion of this subject, see the 
review article |36|. 

In order to identify these challenges, we structure the 
physics of different regimes of the BCS-BEC crossover 
more carefully according to the relevant length scales. 
This will help to find suitable improvements for our 
ansatz of the effective action. We visualize the situation 



in Fig. 29 which serves as a guide for the following 



discussion. We divide the phase diagram into three ma- 
jor regions, and we will argue that the main remaining 
challenges are determined by the few-body, many-body 
and long distance length scales, respectively. 

(i) BEC regime - On the BEC side of the crossover, 
we have a nonvanishing binding energy |s B | m 2 /a 2 . 
In particular, for a — > 0, its value exceeds all other 
scales of the problem, such as the many-body scales ep 
and T. We therefore have a clear separation of scales: 
The physics of the effective boson theory builds up at 
k ~ I /a, and will not be influenced noticeably by the 
presence of the many-body scales which become im- 
portant at k ~ kp, Vr. The many-body physics can 
then be described by a theory of effective bosons, cf. 



Fig 24 Nevertheless, neither the physics of the micro- 
scopic bound state nor the many-body physics of the 
effective bosons is trivial. The first challenge concerns 



the correct description of the bosonic self-interaction, 
i.e. the bosonic scattering length a^. On dimensional 
grounds, = Ba, but the determination of the dimen- 
sionless number B represents a genuinely nonperturba- 
tive problem with no small expansion parameter, despite 
a — > 0. The problem has been solved exactly from the 
four-body Schrodinger equation 11331 . with the result 
B - 0.6, see Tab. [3] This has also been derived from 
phenomenological two-loop self-consistency equations 
111851 11861. Recently, the result was also found in a trun- 
cation of FRG equations taking into account the feed- 
back from an atom-dimer vertex [187|. This trunca- 
tion could be used as a starting point for a many-body 
calculation, which would provide highly accurate low- 
temperature many-body results on the BEC side based 
on the separation of scales and the validity of Bogoli- 
ubov theory for many-body observables. These state- 



ments are illustrated in Fig. 30 where we plot the con- 
densate fraction at zero temperature as a function of 
(kpd) -1 in comparison with an effective Bogoliubov the- 
ory with a phenomenologically assumed molecular scat- 
tering length fl0 = 0.72a. Deviation from the effective 
theory of pointlike bosons occur near resonance, where 
|sb| becomes of comparable size to the thermodynam- 
ics scales. For further comparison, we plot the result 
from extended mean field theory approaches, in which 
= 2a. 

For higher temperatures, the infrared divergences be- 
come more severe, and Bogoliubov theory or its finite- 
temperature extension known as Popov approximation 
||69l fail to predict both the second order phase tran- 
sition as well as the shift in T c 11881 . These genuine 
many-particle effects, which can be traced back to the 
physics of the zero Matsubara mode only and thus the 
3D 0(2) model alone, are already quite well captured in 
the above truncation, see Tab. [3] Higher precision can 
be achieved in an FRG treatment within higher orders of 
the derivative expansion, e.g. 1 1891 1 1901 , or by treating 
the full momentum dependence, e.g. 11 15l [l91 1. 

(ii) BCS and unitary regime - No bound state is 
formed on the BCS side of the crossover and thus the 
many-body scales k F and T dominate the physics in this 
region. An ordering principle is at best given by the 
presence of a strongly expressed Fermi surface, such 
that only modes in its vicinity can contribute to the 
thermodynamics in a nontrivial way. In particular, the 
corrections to the equation of state must remain ex- 
ponentially small. An interesting beyond mean-field 
many-body effect is present in this regime nevertheless: 
While qualitatively the finite temperature transition is 
governed by the log T-divergence in the particle-particle 
channel in the RG flow of the four-fermion vertex, there 




Figure 31: Critical temperature of superfluidity in the BCS-BEC crossover. The limiting cases of the BCS and Gorkov critical temperature for 
negative fermion scattering length are indicated as well as the ideal gas condensation temperature on the BEC side, both with and without shift 
A7" c BEc/7c,id K on ' 3 . The basic truncation scheme for the effective action reproduces the correct limiting boson theory, but overestimates the 
critical temperature on the BCS side. By including particle-hole fluctuations within a rebosonization scheme 1 156], we find agreement with the 
Gorkov result on the BCS side. The corresponding many-body effect is only weakly expressed at resonance and vanishes at the zero crossing of 
the fermion chemical potential. We also plot the critical temperature obtained from extended mean field theory. 




Figure 30: Condensate fraction Clc at zero temperature. The FRG 
result is indicated by the red line. It agrees well with Bogoliubov 
theory of effective dimers with molecular scattering length a^/a = 
0.72. We also show the predictions from extended mean field theory. 



are additional contributions from the particle-hole chan- 
nel which remain regular for T — > 0. They can thus be 
treated perturbatively, and be taken into account via a 
shift of the dimensionless scattering length according 
to akp — > ak-p - a(akp) 2 , with a > a number of or- 
der unity. Taking this shift into account in the exponent 



governing the BCS critical temperature ( 145 1, and treat- 
ing it perturbatively in ak-p as appropriate in this regime, 
one obtains a multiplicative correction of the prefactor 
for the critical temperature e~ na l 2 = (2.2) -1 , resulting in 



! c.BCS 



' c, Gorkov 



2.2. 



(237) 



This suppression of the critical temperature due to 
screening via particle-hole fluctuations is know as 



Gorkov's effect 11921 . The particle-hole fluctuations are 
tied to the presence of a Fermi surface, around which 
these pairs are created. Consequently, the Gorkov cor- 
rection is absent in vacuum and constitutes a true many- 
body effect. 

One may now ask how to recover this effect in the 
FRG treatment |156|. Since the Hubbard-Stratonovich 
transformation of the original fermionic theory has to 
be performed in the particle-particle channel in order to 
capture the formation of the bound state - which is the 
essence of the BCS-BEC crossover - we may think that 
the particle-hole channel is lost. However, there is an el- 
egant trick to recover it with the help of the rebosoniza- 
tion technique j87l |193tf l95 1. More explicitly, we uti- 
lize the fact that, although we have chosen our initial 
values such that the four-fermion coupling vanishes, 



A 



I//.A 



= Abe = 0, 



(238) 



a nonzero value of this quantity is immediately gener- 
ated during the RG flow: d^A^A + 0. The corresponding 
contributions on the right hand side of this equation are 
precisely the particle-hole fluctuations. Thus, an exten- 
sion of our truncation according to 



AT* = 



(«AV) 2 



(239) 



captures the Gorkov correction. Spontaneous symme- 
try breaking during the flow is now signaled by a di- 
vergence of the overall four-fermion coupling A^k ~ 
fijl/m^. Since this is rather difficult to resolve numeri- 
cally, we conduct a Hubbard-Stratonovich transforma- 
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tion on each scale k to absorb the contribution from fluc- 
tuations to A^jc into the running of h^. This procedure is 
called rebosonization; it also allows to apply our stan- 
dard criteria for the determination of the phase diagram 
explained above. 



Eq. (237 1 suggests that our finding for T c /Tp at uni- 



tarity, which is above the experimental value, may be 
lowered substantially by the inclusion of particle-hole 
fluctuations. It is a particular strength of the FRG that 
it allows to answer the question whether particle-hole 
screening is a relevant effect at resonance. In particular, 
the feedback of the particle-hole fluctuations does not 
rely on whether the effect is perturbative (as in the deep 
BCS regime) or not. We simply have to improve our 
truncation accordingly and then solve the flow equation. 
While it is clear that the Gorkov correction is restricted 
to the presence of the Fermi surface /u > 0, and thus 
has to vanish on the BEC side after the zero crossing 
of the fermion chemical potential, it is still an impor- 
tant quantitative question to follow its evolution into the 
strongly interacting unitary regime. It has been studied 
in Ref. 111561 . with the result for the critical temperature 



shown in Fig. 31 From this figure we see that the rel- 
evance of the screening effect diminishes rather rapidly 
as we approach resonance, and in particular is too small 
to explain the large downshift in the critical temperature 
found from QMC simulations II196H199I . 

We now turn to the present status of the crossover 
truncation |200|, which in addition to our above im- 
provements takes into account the renormalization of 
the inverse fermion propagator generated by a diagram 
involving both a boson and a fermion propagator line. 
As a motivation, let us summarize our finding of the be- 
yond mean field effects due to fluctuation effects. The 
boson physics, namely the dimer-dimer interactions, 
drive the physics on the BEC side of the crossover and 
are responsible for the shift Ar c BEC of the critical tem- 
perature with respect to the ideal gas value. However, 
bosons are massive on the BCS side (except in an ex- 
ponentially narrow vicinity of the critical point), and 
thus their contribution is suppressed there. On the other 
hand, the Gorkov correction is bound to the presence of 
a Fermi surface and vanishes on the far BEC side: there, 
the fermions are gapped. In consequence, the diagram 
which renormalizes the inverse fermion propagator re- 
mains suppressed in both BCS and BEC regimes, but 
may be important in the unitary regime in between. The 
renormalization effect on the fermion propagator due 
to the diagram containing both bosons and fermions is 
most relevant at resonance. In each of the limiting cases 
of BEC and BCS, either the one or the other of them is 
massive. In the unitary regime, instead, no such simple 



(a hp)' 1 

Figure 32: Infrared value of the wave function renormalization = 
Z^jc=o introduced in 12401. We observe strong renormalization of the 
fermion propagator close to resonance. For further details see Ref. 
(2001 . 



ordering principle can be applied and a priori there is no 
suppression of this effect. Again, the influence of this 
effect can be investigated systematically by improving 
the truncation: The fermion propagator renormalization 
is taken into account by means of a derivative expansion 
according to 



iW(G) 



■ Z^kiito + q 2 



/4 



(240) 



see M200I . A strong renormalization of the propagator 
is indeed found in the region of anomalously large scat- 
tering length, but nevertheless, the effect on the critical 
temperature at unitarity is not strong enough to explain 
the discrepancy to numerical simulations. 

(iii) Critical domain and universal aspects - In Fig. 
29 we schematically indicate a third region in the phase 
diagram, which is dominated by the critical fluctuations 
on long distance scales k\a » kp, T 1 ^ 2 , IsbI 1 ^ 2 - While 
this regime is very difficult to address with techniques 
other than RG approaches, it is already well under con- 
trol in the simple truncation including boson fluctua- 
tions in our framework. We have seen that the infrared 
boson fluctuations drive the order parameter to zero in 
a continuous manner, cf. Fig. [27] Generically, nonper- 
turbative effects dominate the region around the phase 
boundary in the phase diagram. This regions is expo- 
nentially small in the BCS limit. From a systematic in- 
spection of the size of the universal region of long range 
fluctuations, one finds the largest extend of the critical 
domain at unitarity [35]. A quantitative test for the re- 
liability of the critical modeling is given by the critical 



-(d-2+r,) 



In our 



exponent 77 of the propagator, G(x) ~ r 
basic truncation, it is found to be 77 = 0.05 for the whole 
crossover. This has to be compared to the corresponding 
value T] = 0.038 of the three-dimensional 0(2) univer- 
sality class 111611 . 
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Further universal aspects in the crossover phase dia- 
gram relate to the width of the resonance. Here, we have 
exclusively focused on the broad resonances 11 1 3011 . An 
systematic FRG investigation has revealed the existence 
of two different fixed points governing the Fermi gas 
with divergent scattering length - a broad resonance, in- 
teracting (Wilson-Fisher) fixed point and a narrow reso- 
nance, Gaussian fixed point 111291 . which can be solved 
exactly II 1261 . While the narrow resonance fixed point 
is sensitive to microscopic details, the broad resonance 
fixed point is distinguished by a pronounced insensitiv- 
ity with respect to the precise microphysics. 

Finally, in order to foster comparison with QMC ap- 
proaches, a finite size study has been performed recently 
1 201 1, building on finite size studies in bosonic 1 202 1 
and fermionic theories [203, 204] . Since QMC calcula- 
tions are performed in finite volumes and lattices, such 
an analysis can provide valuable information on the do- 
main of lattice sizes where the extrapolation to infinite 
volume is justified, even if the absolute values for the 
observables were not fully accurate. Complementary, 
the effects of the lattice need to be negligible in order 
to be able to make reliable statements on the continuum 
limit from QMC simulations. For this purpose, the limit 
of vanishing filling has to be taken. This has been criti- 
cally examined in |205, 206 1. 

In summary, we have seen that the BCS-BEC 
crossover shows important nonperturbative effects on all 
scales, ranging form the microscopic scattering physics 
over genuine many-body effects down to the long dis- 
tance critical physics in the vicinity of the finite tem- 
perature phase transition. Full resolution of all of these 
effects requires a unified flexible framework, which is 
provided by the FRG. Beyond being free of intricate 
infrared divergence problems which represent a severe 
obstacle to alternative many-body approaches, this set- 
ting offers a high degree of flexibility for the inclu- 
sion of effects that are well-understood effects in the 
limiting cases, together with the possibility of follow- 
ing their impact when moving into the challenging uni- 
tary regime. This provides a systematic interpolation 
scheme between BCS and BEC regimes, as we have 
seen at the example of the Gorkov effect. It still remains 
to be seen if an effect can be identified which would 
be able to bridge the quantitative discrepancies between 
analytical approaches and Quantum Monte Carlo sim- 
ulations 11196111991 . as well as recent experiments ||43~1 - 
|45). 



5. Outlook 

In these lecture notes, we have given an introduction 
to many-body physics of ultracold atomic systems in 
a functional integral framework. We have worked-out 
the cornerstones of quantum condensation phenomena, 
Bose-Einstein condensation and the BCS mechanism in 
this language. We have also seen how these phenomena 
are connected in the presence of a Feshbach resonance. 

This was done by introducing and applying the con- 
cept of the Functional Renormalization Group, which 
already in a simple approximation to the full quantum 
theory of ultracold atoms allows to access the complete 
finite temperature phase diagram. On the technical side, 
we have seen on the basis of this example how the FRG 
concept can be applied to fermionic and bosonic sys- 
tems in the cold atoms, i.e. nonrelativistic context. In 
view of potential future applications, let us therefore 
come back to the discussion in the introduction, ask- 
ing to which of the challenges mentioned there the FRG 
framework could usefully contribute. 

Resolving physics at different scales and fostering 
comparison with experiment - Experimental tools such 
as Bragg or RF spectroscopy provide information be- 
yond thermodynamics, and in fact yield detailed knowl- 
edge of e.g. the full spectral function for a large regime 
of frequencies and momenta, including strongly inter- 
acting regimes. It is therefore a pressing issue to access 
such observables also theoretically with flexible tools 
beyond mean field theory with quantitative precision. 
Important steps in the direction of a full momentum res- 
olution are made by taking into account higher orders of 
the derivative expansion or performing vertex expansion 
schemes. Both steps have not yet been applied to the 
full quantum theory of ultracold atoms with fermionic 
and bosonic degrees of freedom, but have been tested 
in various general settings, for the higher order deriva- 
tive expansion see e.g. 11891119011 . for vertex expansions 
see e.g. 171 II 15111 171 . Direct calculation of (real time 
domain) spectral functions have been performed in the 
context of bosonic ll631l75Tf77Tl and fermionic 1 1 18 1 sys- 
tems. 

Lattice systems and exotic interactions - Lattice 
models have so far been investigated using the FRG 
framework mainly in the condensed matter context for 
fermionic systems Q. Optical lattices nowadays play 
a key role in the physics of cold atomic systems, and 
in part are crucially needed for a stable realization of 
many of the proposals involving long-range and mul- 
ticomponent interactions. In general, due to the pos- 
sibility of reaching high densities, lattice systems al- 
low to access regimes of strong correlations with rel- 
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ative ease, and therefore offer particularly rich quantum 
phase diagrams. While the study of such quantum phase 
transitions from a low energy viewpoint is interesting 
in its own right IT721 12071 , an outstanding challenge is 
clearly the quantitative assessment of the physics on 
various scales, such as the determination of the loca- 
tion of quantum phase boundaries, where the short dis- 
tance lattice physics needs to be accounted for explic- 
itly. First promising steps in this direction for the con- 
ceptually simplest problem - the Mott insulator to su- 
perfiuid phase transition in the Bose-Hubbard model - 
have been taken in [64 6 51 12081 . 

Non-equilibrium systems - The field of non- 
equilibrium physics with cold atoms is only Hedging 
but there is substantial potential for the discovery of in- 
triguing physics. The need for theoretical tools for the 
description of out-of-equilibrium many-body systems is 
thus pressing, enhanced by the fact that efficient numer- 
ical tools comparable to Monte-Carlo simulations are 
scarce or only applicable to specific circumstances - at 
least in dimensions larger than one. In the context of 
non-equilibrium, closed system dynamics, in addition to 
density-matrix-based approaches [209 -213 1 which have 
proven powerful in understanding aspects of thermal- 
ization, functional techniques 112 141 121511 based on the 
Keldysh real time path integral are most promising due 
to their high degree of flexibility in describing physics 
at different scales in one unified framework. FRG ap- 
proaches to this class of problems have been put for- 



ward in 12161 12171 for the investigation of the real- 
time evolution starting from a given initial state and 
in 1121 81 . applied to the analysis of strongly nonlinear 
wave-turbulent states far from equilibrium. Another in- 
teresting direction is provided by open systems, which 
in contrast to the dynamical phenomena above can ex- 
hibit stable non-equilibrium stationary states - they 
share exact time translation invariance with thermody- 
namic equilibrium, but are governed by different distri- 
bution functions. One can therefore hope that FRG ap- 
proaches based on the Keldysh formalism can still give 
analytical insights into such problems. Theoretical ap- 
proaches for classical non-equilibrium stationary states 
have been worked out in [219-222]. Non-equilibrium 
stationary states in cold atomic quantum systems are 
currently moving into the focus of research, being re- 
alized e.g. via the competition of particle loss and re- 
pumping in physical contexts as diverse as ensembles 
with optical Feshbach resonances II 1 241 . low dimen- 
sional systems of polar molecules [223 1, and metastable 
repulsive fermions [224 J, or in stationary states result- 
ing from tailored dissipation ||30l |3T1 [2251 . 
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Appendix A. Functional integral representation of 
the quantum partition function 

In this appendix, we derive the functional integral ex- 
pression for the partition function of a quantum many- 
body system. The construction utilizes coherent states. 
These are eigenstates of the annihilation operators a, 
and allow for a parametrization of Fock space, which 
is different from the occupation number representation. 
We show that bosonic or fermionic particles can be for- 
mulated in terms of a nonrelativistic field theory with 
euclidean time t. Whereas the former are expressed by 
a complex field <pi(j), the latter correspond to Grass- 
mann valued fields tA,( T )- We give a brief introduction 
to Grassmann numbers and the corresponding calculus. 

To begin our analysis, we note that the Hamiltonian 
of a many-body system can be expressed in terms of 
creation and annihilation operators, a ' and a, according 
to 



H = H(a\a) = Y 



tija'Mj + 



^ Vi jk ia]a\a k ai + 
ijki 



(A.l) 

As we have discussed in section |2.1[ this generic form 
of the Hamiltonian is important for cold atoms, where 
only two-body interactions (i.e. VW) are relevant. In 
Eq. (A.l i, all annihilation operators are to the right of 
the creation operators. If this is the case, we say that the 
operator is normal ordered. The labels i, j, . . . run over 
a particular choice of single particle states. Typically, 
we think of spin and momentum or points on a lattice. 
For a finite system, this enumerates a discrete and finite 
set. However, it is also common to write dp^ or a? with 
continuous variables p and x, keeping in mind that one 
has to go back to a discretized formulation if problems 
should occur. 
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The creation and annihilation operators satisfy 
commutation (anti-commutation) relations for bosons 
(fermions). We have 

[a,,a\] = a,at — a'a,- = 5n, [a,, a,] = (bosons), 

j J J J J 

{«,,«'.} = a,a'. + a'.fl; = 6a, {a;,a;} = (fermions). 

j j j J J 

(A.2) 

The particle number operator can be expressed as 

N = N{a\a) = Yaja,. (A.3) 

i 

We assume that there is a vacuum state |vac) which 
does not contain any excitations. Consequently, we have 
a,|vac) = 0. Since the states in occupation number rep- 
resentation form a basis of the many-body Fock state, 
we can write the partition function Z = Tre~^ H- '^ as 



1"2 ■ 



(A.4) 



l"l»2-> 



From the last formula we easily obtain Eq. ( 20 1 for 



the pressure of a noninteracting gas. However, for an in- 
teracting system, it will in general not be possible to ex- 
actly calculate the partition function. Therefore, we aim 
to rewrite Eq. (A.4i in terms of a functional integral, 



which is particularly well-suited for treating interaction 
effects in a systematic fashion. 



Coherent states for bosons 

For a functional integral representation, we would like 
to substitute the operators at and a, for corresponding 
classical fields (pi, which are then quantized. We first 
restrict ourselves to the bosonic case. The natural way 
to replace an operator by a number is to let it act on an 
eigenstate. We ask whether there are states \<p\<fi2 ■ ■ ■ ) 
with complex numbers ipt e C such that 



ai\(pi(f2 ■ ■ •) = <fi\<p\if2 ■ ■ ■) 



(A.5) 



for every i. We call such a state a coherent state. The 
(pi are not necessarily real, because d, is not selfadjoint. 
The creation operator a\ cannot have an eigenstate. In- 
deed, assume there was an eigenstate |*F) of ci\ . Then, as 
every Fock space state, \ x ¥) could be expressed in the oc- 
cupation number representation as a linear combination 
of several basis states \ri\nz . . .). In this linear combi- 
nation, one basis state has the smallest particle number 
Yji n,. Now, letting at act on I^F), the particle number of 
every basis state in the superposition gets increased by 



one. Due to the fact that this also increases the minimal 
particle number by one, ci\ \*¥) cannot be proportional to 
PP>. 

We claim that, for every collection {(/?,} of complex 
numbers, 



|^> = |^ 2 ...> = e 2-^|vac> 



(A.6) 



constitutes an eigenstate of {a,} with eigenvalues {<pi}. 
Applying a, to \ip), it will commute with all at for i + j. 
We can then - for fixed i - use the relation [a, = 
n(a')" -1 . This yields 

a\ip) = ae^lvac) = [a, e^]|vac> 

oo n 

= V ^-[^(a'nivac) 

= * X 7^T T (a 1 ")^ 1 |vac> = ip\tp), (A.7) 
(n - 1)! 



«=i 



which proves a,|<^) = <fi\(f). Via complex conjugation 
we find (ip\aj = (<p\<p*. 

Clearly, for every choice of complex numbers {i^,}, 
we can construct a coherent state. In contrast, the oc- 
cupation number basis |«i«2 • ■ ■ ) of the Fock space was 
limited to integers {«,}. This increase of degrees of free- 
dom results in the fact that coherent states are over- 
complete. They do not represent an orthonormal basis 
of the Fock space, but nevertheless any state can be ex- 
pressed as a superposition of coherent states. With the 
explicit formula in Eq. ( |A.6[ ) we find for the overlap of 
two coherent states 



(A.8) 



Functional integral 
One can show that the weighted sum 

D^*Dy>e" z ^>'»<^| = 1 (A.9) 



constitutes the unit operator in Fock space. The func- 
tional measure is defined as 



dipjdip* p-r d(Re^,)d(Im^,) 



2m 



n 



(A. 10) 



The integration over all values of ifj is similar to the 
summation over all occupation numbers, but due to the 
over-completeness we have to suppress the individual 
contributions by a weighting factor. 
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Since we assumed the Hamiltonian to be normal or- 
dered, we have 

{ip\ (H - fift) \<p) = {ip\(H{a\ a) - ^{a\a))\<p) 
= (H(tp'\(p)~fMN(<p'*,(pj)e I '>'e><. (A. 11) 

Obviously, we arrived at our goal to substitute the oper- 
ators for complex numbers. There are no operators ap- 
pearing in the functional integral. We could now insert 
the identity operator into Z = Tre - ^"^ and express 
the partition function as a path integral. However, there 
is one complication: The exponential of the normal or- 
dered operator H-fiN is no longer normal ordered. De- 
noting normal ordering by double dots we have 



e + . e 



(A. 12) 



The situation is less severe if we introduce a small pa- 
rameter e, because 



e .sA: = . e sA. + 0(c zj_ 



(A. 13) 



We therefore divide the inverse temperature /3 into M 
small intervals of length s such that ft = Ms. We then 
insert M— 1 resolutions of the identity. This introduces a 
discrete label s = 1, . . . , M according to pi — > {(f S i} s . It 
will later be associated with imaginary time. We arrive 
at 

ZQx, T) = Tre^"-^ 
= j D V lD Vo e-^^{ n \e-^%o) 

= J ,1V,, ;...IV (| IV„, 

x {^e-^%^) . . . (<pr\e-^-»% ) 

r M 
= I D^p*D^exp|-e^[^^ s i(— — ) 



+ 



(ffC^;,^!)-^;,^!))]} 

D^D^e -5 ^** 1 . (A. 14) 



In the second to last line we introduced a new complex 
field ip S i which depends on both i and s = 1 , . . . , M = 
Pie. The functional measure is naturally extended to be 



2n\ 



(A. 15) 



For these fields we introduced the action 

M _ 

+ 



,<P!-i)-frtf(<P*><Ps-i))]- 

(A. 16) 



The condition </) M( - = ^ for all z originates from the fact 
that the partition function is a trace. 

Sending M — > oo and e — > while keeping ft = Ms 
fixed, the discrete variable ss = r becomes continuous. 
The partition function then acquires the form 



D<p*D(pe 



-S [<p',<p] 



(A. 17) 



with microscopic action 

S [<p\ ip] = J o dr( ^(r)(<9 T - aiM(t) 

+ //(^(r),^(r))). (A. 18) 

Recall that z can be a continuous variable, too. Typi- 
cally we are interested in i = x and fields tp(T,x). If an 
expression happens to be not well-defined during a cal- 
culation, we can always go back to the discretized form 
of the action in Eq. (A. 16 1. The condition ipi(J3) = (fi(0) 



restricts the functions tp and tp* to be periodic functions 
in r. 

To summarize, we found a functional integral rep- 
resentation of a generic many -body system of bosons. 
The non-commutativity of operators introduced the time 
variable r. 



Coherent states for fermions 

How can these considerations be extended to include 
fermions? Most of the formulae from above remain 
valid or only receive corrections due to some signs. 
However, there is one important conceptual difference 
in the path integral representation for fermions. 

We already noticed that no operators appear in the 
functional integral formulation. In the bosonic case, 
complex numbers took the positions of the (normal or- 
dered) annihilation and creation operators. In order to 
satisfy the anti -commutation relations for fermionic op- 
erators instead, the eigenvalues of the {a,} cannot be 
complex numbers. Indeed, if \<p) = \<p\<p2...) was a 
coherent state of arbitrary complex numbers, we would 
arrive at the contradiction 

= \a h aj}\<p) = ((pupj + <pj<pd\<p) * 0. (A. 19) 
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Complex numbers thus cannot be applied for con- 
structing coherent states of fermions. However, there 
are objects which can formally be multiplied and obey 
anti-commutation relations. They are called Grassmann 
variables (or Grassmannians) and a calculus can be de- 
veloped for them. We will reduce our discussion of this 
issue to a minimum and refer the reader to the textbooks 
iBUl lM | for more details. Since any two Grassmannians 
iff and rj satisfy 

if/T] = -rpff, (A.20) 

the construction of functions of them is pretty sim- 
ple. Assume we have a polynomial expression f(X) - 
fiX' with e C. We then define 

M) = Jb+M, (A.21) 

which incorporates the condition iff 2 = iffiff = -iftiff => 
iff 2 = 0. Functions f(X u . ..,X n ) of more than one vari- 
able, and analytic functions which allow for a series 
expansion, can be extended to Grassmann-valued argu- 
ments as well. The associated series will always termi- 
nate at a finite order and we never have to bother about 
convergence problems. We will often have to deal with 
expression of the form 



1 + Clfj 



(A.22) 



for c being either complex or Grassmann valued. In the 
first (second) case, we have f(X) = e cX ( f(X, Y) = e XY ). 

The product 7r of an even number of Grassmannians 
satisfies mft - iffn for any Grassmannian iff, because 
the minus signs from anti-commuting through iff cancel 
each other. We refer to such products as being even. 
Complex numbers are trivially even. The remaining 
Grassmannians are then odd. 

We have just seen that functions of Grassmann vari- 
ables can be linearized and it is thus possible to define 
a differentiation rule, although there is no actual way in 
which a Grassmannian is either small or large. (In fact, 
they do not have particular values like complex num- 
bers, but will rather only serve for generating correla- 
tion functions.) We introduce a left- and right-derivative 
according to 



d_ 
dip 



(fo + <ffi)=fu 



(fO+M)gT =/l- 

Note that f x = f x (f x = -f{) if /, is even (odd). 



(A.23) 



Integration of Grassmann functions is defined to be 
linear and to satisfy 



J diffl = 0, j diffiff = 1. (A. 



24) 



Thus, it coincides with left-differentiation. Note how- 
ever that diffdn = -dndiff. 

We can now formulate the coherent state path integral 
for fermions in terms of Grassmann fields (^(t). It is 
given by 



Z(jx,T) 



J 



Diff*Difte 



-s[p.<i>\ 



(A.25) 



><ifM=-<i<m 

The action S [iff*, iff] is constructed from the Hamiltonian 



analogous to the bosonic case in Eq. (A. 18 i. The func- 
tional measure reads 

DrT>iff = Y]diff* si diff s<i . (A.26) 

s,i 

In contrast to the bosonic case, the partition function 
is now restricted to anti-periodic functions in time- 
direction, satisfying 



iffW = -<A,(0) 



(A.27) 



for all i. 

The (anti-)periodicity reduces the Fourier transfor- 
mation of the fields in time direction to a Fourier series 
with discrete frequencies. Indeed, for a function f(x) 
satisfying fix + L) = f(x) we have 



facwxjh 



(A.28) 



neZ 

Thus, in our case, we arrive at 



ip(T,x) = TY f Al^e^^ipM (A.29) 



with Matsubara frequencies 



2nnT (bosons), 
2n{n + 1/2)T (fermions). 



(A. 30) 



In the zero temperature limit, the stepsize Aoj = 2nT 
gets infinitesimally small and the sum is replaced by a 
Riemann integral over the continuous variable to. More 
explicitly, we have 

(A.31) 
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Appendix B. Lattice magnets and continuum limit 



periodic boundary conditions we have 



J^mitrij = -- ^(m,-m /+1 + m/m,^) 



In this appendix, we apply the concepts of section 3.1 
to the Ising model. The formulation of the system on 
a discrete and finite lattice allows for a transparent dis- 
cussion of spontaneous symmetry breaking and the con- 
struction of the effective action. We show, why sponta- 
neous symmetry breaking can only occur in an infinite 
system and in which sense the partition function and the 
effective action store the same physical information in 
a different manner. We then explain how our findings 
on the lattice can be translated to a continuum theory 
such as ultracold bosons and fermions. In particular, we 
give a brief overview onto functional differentiation and 
Gaussian functional integrals. 



Spontaneous symmetry breaking: Magnet case study 

We consider spins son a cZ-dimensional lattice. To each 
spin, we associate a magnetic moment m according to 
m/HB = ghs/h, where [i% and g^ are the Bohr magneton 
and the Lande factor, respectively. Choosing units such 
that m 2 = 1 we obtain the classical Heisenberg model 
of order «, where n is the number of components of in. 
The corresponding Hamiltonian is given by 



ijMi -nti + 



ntj 



(B.l) 



<U> 



with i labeling the individual lattice sites and the sum- 
mation (i, j) being restricted to nearest neighbors. For 
Ju > 0, the first term favors magnets on neighboring 
sites x; and xj to be aligned by lowering the energy of 
such a configuration. The magnetic field hi in the sec- 
ond term can be regarded as a source of magnetization 
at the site Xj. 

For the purpose of our analysis, it is sufficient to re- 
strict to the case of n — 1 such that the orientation of the 
magnets can be either up or down. The system is then 
known as the Ising model. The Hamiltonian function of 
the latter is given by 



H 



-J fflifflj 

<u> 



(B.2) 



with m, being either +1 or — 1. Moreover, we assumed 
Jij — J > to be an overall constant. The first term can 
then be interpreted as a kinetic term. Indeed, assuming 



Z m '( 



m, + i - 2m/ + rrii-\ 



Jo 2 



(B.3) 



with lattice spacing a, number of magnets N and dis- 
crete Laplacian V 2 . The constant shift is irrelevant here. 
For cold atoms it can be absorbed into the definition of 
the chemical potential. 

Assuming the system to be in equilibrium with a heat 
bath of temperature T, the orientation m, of the mag- 
nets on the individual lattice sites can be treated as a 
stochastic variable with a canonical probability distri- 
bution. We decompose m, into its mean and fluctuating 
part according to 



m: + 5m; 



(B.4) 



with (Smi) = 0. The mean magnetization is given by 



rhf = (ra,),, = i ^ m v 



1 d 
p Oil, 



Thus, taking a derivative of the partition function with 
respect to the ith component of the set {hi}, we obtain 
the magnetization at site i. 

If we set hi = after taking the derivative in Eq. 



(B.5 i, we can calculate the spontaneous magnetization 
at vanishing field, (m,)/, = o. Since the Hamiltonian for 
hi = is symmetric under m-, — > -nti (^-symmetry), 
we expect (m,)/, = o to be zero, because each configura- 
tion with a certain number of spins up has a correspond- 
ing configuration with all spins reversed and they should 
cancel when summing over all configurations. This rea- 
soning, however, can only be applied to finite systems 
(finite number of lattice sites). Indeed, the Ising model 
shows phases of nonvanishing spontaneous magnetiza- 
tion in the infinite volume limit, most prominently in 
the exact solution of the two-dimensional Ising model 
by Onsager. 

The magnetic field singles out a preferred direction 
at each site, because parallel alignment along the field 
minimizes the energy. If we first enlarge the system to 
infinite volume and then remove the field, we have 



I mo + ("broken phase") 
("symmetric phase"). 



(B.6) 
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Interchanging these limits we always get zero, because 
we can apply our reasoning from above. Note that 
Eq. ( |B.6| l already hints towards the answer to the deep 
question why there actually are phase transitions in na- 
ture. The calculation of the partition function in a finite 
volume consists of a summation of a finite number of 
terms, all of them being nice analytic functions. Even if 
the Hamiltonian contains a symmetry breaking term like 
Yji hjirii in the above example, the contribution of this 
terms to Z({/i,}) can be removed at the end of the calcu- 
lation, because Z is analytic in the /i,. However, if we 
calculate the corresponding intensive potential i log Z 
(e.g. free energy density or pressure) in the thermody- 
namic limit V, N — > oo we may introduce singularities 
and non-analytical behavior in Z({/;,}). In fact, the sum 
of infinitely many analytic functions will in general not 
be an analytic function again. Therefore, the thermody- 
namic potentials remember the field hi, even if we set it 
to zero at the end of the calculation. 

Leaving this aspect aside for the moment, we want 
to extract further information from Eq. (B.5i. Higher 



correlation functions can be obtained by taking higher 
derivatives with respect to the field. For instance, the 
connected two-point function 

((to,- — fhi)(nij — fhj)) — (mirrij) — fhifhj 

logZ({A,-}) (B.7) 



dhjdhj 



tells us how strongly deviations from the mean value at 
site i are correlated to deviations at site j. 

We now invert the magnetization problem. Given an 
arbitrary mean magnetization to,-, how do we have to 
choose {hi} in order to obtain exactly this magnetiza- 
tion? The corresponding generating function is obtained 
from the Legendre transform of logZ({/i,}) in the vari- 
ables {hi}. We build 

F({to,}) = sup(V rhihi - logZ({A/})) 

Vij) j 

= Y j m i h i -\ogZah i }). (B.8) 



In the second line, which is valid for continuously dif- 
ferentiable log Z, the magnetic field is defined implicitly 
through the equation m,- = diogzdfr}) Nctfe that the right 
hand side of Eq. (B.8 i does not depend on {hj}. We ar- 



be chosen for given mj, via 



dm 



dhi d log Z 
dfhj dfhj 



dlogZ dhi 



dfhj 
dfhj 




Functional differentiation 

Instead of studying fields which live on discrete lattice 
sites Xi, we now turn our attention to a continuum theory 
with space-time variable X = (t, x). This is of relevance 
for ultracold atoms considered in the main text. For a 
system of classical magnets as treated in this appendix, 
we do not have a time variable t. However, the func- 
tional methods developed here are not invalidated and 
we may replace X — > x in this case. The sets {/i,} and 
{to,} become functions h(X) and m(X). Sums 2; are re- 
placed by integrals JT and instead of partial derivatives 
d/dfhi we take functional derivatives 5/Sm(X). We for- 
mally write 

2->X, d-*6. (B.10) 
More explicitly, we have 



6m(X) 
8fh(X') 



= 6(X - X') = 5(t - t')6 w (x- x*), (B.ll) 



which mimics the well-known relation dfhi j dfhj = Sij. 

For a proper functional differentiation, we have to 
specify which variables are treated as independent of 
each other. Certainly, the field values on different space- 
time points (lattice sites) are distinct variables. In addi- 
tion, for a complex field fh(X) = (p(X), as it is the case 
for ultracold bosons, we have to specify whether (<p, <f>*) 
or (<p\,(f>2) are treated as independent variables. Here, 
the real fields 4>i and <p2 are derived from the represen- 
tation <p = ^(0i + i(p 2 ). 

Let us choose (f> and <p* as independent variables. As 
a generalization of 



a 



we have 



rive at the answer to our question, how the field has to 



6 
6(f> z 



- 4>* x A X y4>y = I i 

7. JXJ JY 



i Z Y<PY- 



(B.12) 



(B.13) 
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The expression Axy is an operator kernel and thus simi- 
lar to a matrix element. If the kernel contains derivative 
terms, we sometimes have to perform a partial integra- 
tion before taking the functional derivative, e.g. 

4- f f z (d T - V 2 )0 Z = -A- f M-dr ~ V 2 )^ 
' Jz 0(f>x Jz 



Given the above restrictions on A, we find both the left 
and the right hand side of this equation to be an analytic 
function in the complex coefficients Ay of A. Thus, by 



"x Jz 



Sep, 

(-d T -V 2 )cp x . (B.14) 



The second functional derivative of this expression is 
found to be 



virtue of analytic continuation, Eq. (B.19i is also valid 
in the case of complex A. 

For complex integration variables (f> = ((pi ,4>n)> 

the denominator of Eq. (B.19i appears without the 
square root. Introducing in addition complex source 
fields j = (ji, . . . , j N ) on each lattice site, we find from 
completing the square that 



Ko<px Jz 



1>z = (-d T -V 2 ) x 8(X-Y) 



1 



det(A) 



(B.15) 



(B.20) 



We emphasize that this is a kernel and not an operator 
acting to the right. This may be more transparent from 
writing 

f <f> m Y (d T - V 2 )0 y = f f Y 5(X - Y)(d T - V 2 ) x 4> x 
Jy Jx,y 

[(-d T -V 2 ) x 6(X-Y)]<p x , (B.16) 



Here, f = (j*)' has the usual meaning. 

In practical calculations, the matrix A (e.g. the in- 
verse propagator) is often given in terms of real basis 
fields. One might then either work with real variables 
as in Eq. (B.18 i or with complex ones as in Eq. ( B.20| >. 
To see this, we write 



JX.Y 



which allows for reading off the second functional 
derivative immediately. 



Gaussian functional integrals 

The partition function which corresponds to a Gaussian 
probability distribution 



with 



B = 



A' 
A 



(B.21) 
(B.22) 



The integral on the right hand side of Eq. (B.21 1 is of 
type ( |B7T9] > with m, = (0„ 0*) for i = I,..., 2N. We 
then find 



i r N C(T1 d W f \ -<b' 

(O) = -J(f]d^)O(WK^ (B.17) J\[[ 2m > 



can be computed for an arbitrary complex symmetric 
matrix A with nonvanishing eigenvalues Aj + and 
Re(A) > 0, i.e. m'Re(A)m' = m,Re(A), ;m'. > for 

J J 

all choices of m { and m\. We have 



d<p*d<pje 



^ni' VdetB 



1 



1 



1 



iJV [(-l)A'(detA) 2 ] 1 / 2 detA' 
(B.23) 



-^mjAjjitij 



(B.18) 



For real A, we can find an orthogonal transformation 
O such that OAO' = dmg(A u . . . , A N ) = B. Writing 
m'Am = (Om)'B(Om), we observe a transformation of 
the integration variable according tomn Om to factor- 
ize the integral into Gaussian integrals. Indeed, the 
Jacobian of the transformation is unity. We arrive at 



Thus, in both cases, the solution is found to be (detA) -1 
or (detB) -1 ^ 2 , respectively. In particular, applying the 
formula log detA = TrlogA, which is easily seen by 
writing both sides of the relation in terms of the eigen- 
values of A, we find up to an overall constant 



TrlogA = -TrlogZ?. 



Z n = 



(2;r)^ 2 (2^ 



i • • ■ an 



Vdet(A) 



(B.19) 



2m ' "2 

(B.24) 

For this reason, the factor of 1 /2 in the one-loop formula 



53]l or the flow equation ( 165 1 for the effective action 
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might be present or not, depending on the definition of 

S< 2) . 

The generalization of these lattice formulae to the 
continuum are straightforward. For bosonic atoms, rep- 
resented by the complex field <p(X), which are coupled 
to a complex source field j(X), we have 

J D<£*D</>exp(- J (ip* x A X Y<PY + fxfx + f* x jx)) 

= di) exp (/^ (A " w )- (B - 25) 

The functional measure can formally be written as 



Dip*Dip 



d<p*(X)d<p(X) 
2n\ 



(B.26) 



the precise definition being given by the discretized 
version. Note that infinite constant prefactors such as 
N = Ylx 2ti are not of relevance for our purposes, be- 
cause they always drop out in calculations of correlation 
functions. We can thus normalize them to be unity right 
from the beginning. 

For Grassmannians, although they do not have con- 
crete values, a complex conjugation can be defined. 
However, for all purposes of our interest, it is suffi- 
cient to work with two independent Af -vectors iff - 
(iff i, . . . , i[in) and if/* = (i//* v . . ., tf/L). We still write 
ijr = (iff*)', but i/r is in no way related to i/r. After intro- 
ducing two source terms rj and 77* we find 



D^Dif/e 



-^MW^n _ det(A)e'' t(A_1),7 . (B.27) 



Note that the determinant appears in the numerator. 
This peculiar feature of Gaussian integrals for Grass- 
mannians is related to the definition of the exponential 
function. We used the functional measure 



(B.28) 



Of course, formula (B.27i immediately applies to the 



continuous case of Grassmann fields i//(X) and if/*{X) 
with the usual replacements, for instance 



rf A 77 



JX,Y 



rfx( A 1 )xytjy- 



(B.29) 



Appendix C. One-loop effective potential for 
bosons 

This appendix provides computational details on the 
derivation of the one-loop effective potential for weakly 



interacting bosons. In particular, we derive the conden- 
sate depletion at zero temperature, which is a pure inter- 
action effect. 



Evaluation of the one-loop correction 

Here we perform the calculation of the effective action 
in the one-loop approximation. This computation pro- 
vides some useful formulae and serves as a basis for ex- 
tracting the condensate depletion present in weakly in- 
teracting Bose systems at zero temperature. It will also 
shed light on ultraviolet divergences, and how to cope 
with them. 
We start from 

r (i-ioo P ) M =sw + ATW =S[<f>] + iTrlogS (2) [0]. 

(CI) 

The field <p(r,x) is arbitrary. By taking functional 
derivatives with respect to <p, we obtain higher corre- 
lation functions F (n) to one-loop order. If we evaluate 
p(i-ioop)|-0j f or a constant field, we arrive at the effec- 
tive potential t/ (1 ~ loop) (0). Analogous to Eqs. (fssl and 



( 89 1, we then find the gap equation and equation of state 
according to (dU / d</>)(<po) = and (dU /d(i)(<po) = -n, 
respectively. 

We express the microscopic action in the basis of real 
fields^ = -^(ipi +i(p 2 ) as 



S[fi,<f2l = I U A (p x ) 
Jx 



-v 1 

2 M 



fix 
fix ) 



(C2) 



We introduced the microscopic or classical potential 



^ 8 2 



(C3) 



which depends on p — tp*<p = + <p$). It coincides 
with the classical effective potential from Eq. (85 1. The 
label A indicates that this expression for the effective 
potential is only valid at the energy scale of the UV cut- 
off, whereas it is changed due to the inclusion of fluctu- 
ations in the infrared regime. 

We calculate S (2) [<p* , <p] by expanding S [<f>* + dtp* , (p + 
S(p] to second order in dtp around a real, constant field 
<p = yfp. This allows for the implementation of spon- 
taneous symmetry breaking and to distinguish between 
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amplitude and phase fluctuations. We find 



matrix. We have 



+ 2 I ( S <PlJ> &<P2jd 



Pu,x id T 
-id T Prix t 



with 



Pujc = -^ + U' A (p) + 2pU^(p), 



8<pi,x 

5tfi2,x 

(C4) 



(C.5) 
(C.6) 



Here, a prime denotes differentiation with respect to 
p. Note that the minimum po of the full effective po- 
tential U(p) will usually not satisfy U' A (po) = 0, i.e. 
po(p, T) + pig, because the inclusion of quantum and 
thermal fluctuations shifts the position of the minimum 
in field space. Otherwise, no phase transitions would 
occur in nature. However, in the case of a weakly inter- 
acting Bose gas at low temperatures, the shift is small 
and can be treated perturbatively. 

In order to deal with the derivative terms in Eq. 
( |C.4| >, we transform our fields to momentum space Q — 
(co„,q). Note that the Fourier transformation of Eq. 
\CA) would be difficult, if the expansion point <p was 
not constant in space-time. We have 6ip*(Q) - Sipj(-Q), 
because 6ip\ and 6tp2 are real. This leads us to the ex- 
pression 



S[y/p + S<p u S<p2\^S[y/p,0] 
+ i£(^ 1 ,_ G ,^ 2 ,_ e )G A 1 (0 

with classical inverse propagator 
( 



S<f2,Q l 



(C7) 



s q + U' A (p) + IpU'^p) -a>„ 

0>n Sq + U' A {p) 



Sq-fi + gp t 



(C8) 



We also introduced the notation f = T^n J F° r 
the corresponding zero temperature limit we refer to Eq. 
COT) , 

The trace of the inverse propagator G A '(0 consists 
of an integration over Q and the usual trace of the 2x2- 



Tr log S (2) [ y/p, 0] = 5{P = 0) f (log G A ] (g)) fl 

Jo 



■■pv 



logdetG A I (0. 



We used here the matrix identity 

trlogA = logdetA, 



(C9) 



(CIO) 



which is easily seen by writing the left and right hand 
sides in terms of the eigenvalues of Ay. Moreover, we 
used that 



8{P 



= Q) = j e iPX \ p=0 = j dr J d 3 x=/3V (Gil) 



coincides with the volume of space-time. 

With these formulae, we find for the one-loop effec- 
tive potential 



U {Uoop) (p) 



^ 8 2 ^ 
-W + 2 P + 2 



2L 



logdetGxkfi) 



~W + 2 P + 2 



) 



C d 

= ^ + ip* + T J _ -L log sinh(£ 9 (p)/2r). 



(C12) 



The bosonic Matsubara summation was evaluated ac- 
cording to Y^n l°g(l + (^2) - 2 log sinh(x). We dropped 
an overall constant in line two, which is irrelevant for 
the thermodynamics, and introduced the abbreviation 



E q(P) = J( s g~f l + 3gp)(£ g - M + gP)- (C.13) 



Forp = p/g we find E q (p) = y/e q (e q + 2 gp~oj = E q . The 
integral on the right hand side of Eq. (C. 12 1 is diver- 
gent for large momenta. Regularization via introduction 
of a UV cutoff A yields terms proportional to A 3 and 
A. These unphysical divergences from perturbation the- 
ory are cured in a physically motivated ultraviolet (UV) 
renormalization scheme based on an investigation of the 
gap equation and the equation of state. 



Thermodynamics, condensate depletion, and UV 
renorma liza tion 



From the one-loop effective potential ( |C.12| i we obtain 
the phase structure and the equation of state. The order 
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parameter po(p, T) is found from the gap equation 



<9£/(l-l°op) 

0= (fio(M,T)) 



dp 

= -V- + gPo + 



-H + gp Q +g 



1 r d\ 

2 J Jin 

f — 

J (2n 



dE 



coth 



\2r) 



q s q -p + |gp 



(C14) 



The integrand tends to 1 for large momenta. The phys- 
ical origin of this divergence can be understood from 
the fact that we have assumed interactions which are 
local in coordinate space, and hence constant in mo- 
mentum space up to arbitrarily large momenta. It can 
be cured by a renormalization of the coupling constant 
gR - g + dg(A), where 6g(A) is an explicitly cutoff de- 
pendent term, which replaces all the microscopic details 
we left out during our calculation. (The prescribed UV 
renormalization procedure is performed for the case of 



fermions within BCS th eory in Eq. ( 143 1.) We then ob- 
serve from Eq. ( C. 14 1 that po = p/g + 0(g ) at zero 
temperature. 

Taking a //-derivative of Eq. ( |C.12[ ), we obtain 



8U 
dp 



(P) 



-p 



1 r d'q 

2 J (2nf 



d 3 ^ s q -p + 2gp jE q (p)\ 
coth - 



\ 27 / 
(C.15) 

At zero temperature, the coth can be replaced by unity. 
Inserting po - pig we then find 



n(p) 



"dp P =P ° ' "' 



1 r d'q 

2 J (271)1 



q Sq + gPO 



(C.16) 



which again diverges due to high momenta in the inte- 
gral. Writing the particle number as 

r d 3 

n = <p* <f>o + J ■^-^(dipldip^, (C.17) 

we obtain a recipe to cure this divergence: it relates to 
the fact that the functional integral works with fields 
rather than operators and thus does not contain infor- 
mation on the operator ordering which is present in a 
second quantized formulation. Indeed, we have the re- 
lation (6<p*6<p) = (SipStp*), We may contrast this with the 
particle number density defined in second quantization 
as n = (j>* Q <p~ + Tjqi=o(^^q)' where a^, are bosonic op- 
erators. We can reconcile these two approaches by iden- 
tifying permutation invariant combinations, which must 
be equal to each other: (dtp* dtp) = I (6<p* dip + 5tp5ip*) = 
Iffia+M') - (a'fl)+i for each momentum mode. We 



thus conclude that we overestimated the contribution of 
each mode by 1 /2 in the functional formulation, leading 
to a linear UV divergence in Eq. ( |C.16| l and a cubic one 
in the corresponding effective potential. Physically, we 
have thus identified this divergence as resulting from the 
quantum mechanical zero point shift for each oscillator 
mode <? 111271 . Subtracting this term we finally arrive at 



n(p) 



= <Pq<Po 



= Po + 



J (2tt) 3N i q/ 

1 r d 3 <? i B q +gp _ \ 

2 J (2tt) 3 1 E q J' 



(C.18) 



The resulting condensate depletion is not found in a 
noninteracting Bose gas. We can interpret this behav- 
ior as part of the bosons being kicked out of the con- 
densate due to the repulsive interactions. This is an ob- 
servable effect of quantum fluctuations, which occurs in 
the absence of thermal fluctuations at T = 0. In con- 
trast, superfluid density and particle density are equal 
at zero temperature. This is ensured by the Ward iden- 
tity which is related to Galilei symmetry of the effective 
action. Although the condensate spontaneously breaks 
Galilei symmetry, this does not invalidate the statement, 
because the Ward identity is a property the effective ac- 
tion. 

Appendix D. Symmetries of the effective action 

A global continuous symmetry of the classical action 
yields a (classical) conserved charge - this is Noether's 
theorem. In the absence of anomalies, there is also a 
conserved charge for the full quantum theory. Here, we 
briefly review the formalism for the construction of the 
conserved Noether charge for both the classical and the 
quantum case. 

The classical equations of motion are found from the 
action S[ip] and the corresponding Lagrangian the 
latter being defined by S — J x -C, from the variational 
principle 



= 



6S 
8<p(X) 



[<p< 



5Z 
S<p(X)- 



(DA) 



In particular, for a Lagrangian X. = -C(tp,d^tp) which 
only depends on the field and its first derivatives, we 
have 







ss 

dip 



8X 
dip 



d£ 



d(dn<p)L =i 



(D.2) 



f=fa 
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These are the standard Euler-Lagrange equations. In 
general, the Lagrangian may also contain higher deriva- 
tive terms of the field, which we indicate by writing 
£ = £{(f, djj<p, ...). 

Eq. ( |D. 1 [ > tells us that the solution (po is a stationary 
point of the action with respect to variations of the field. 
In particular, these variations <p — > ip + dtp may be gen- 
erated by a continuous transformation ip — > (pf a \ where 
a(X) is a continuous function quantifying the mapping 
such that a — corresponds to the identity map. We 
restrict the discussion to transformations, whose linear 
part in an expansion in powers of a is at most linear 
in the field <p, A typical example is provided by a U(l) 
transformation according to ip — » ip^ = e la <p ^ (1 +\a)<p 
with real or. Defining 

S (a) [<p] =S[ip (a> ] = f £ (a \ (D.3) 
Jx 

which is a functional of a(X), we apply the chain rule 
for functional differentiation and specialize Eq. ( D. 1 1 to 

SS (a) 



a local continuity equation d^J^ = for the classical 
Noether current 



J cl 



(D.8) 



It is related to the global conservation of the Noether 
charge 



Qc 



(D.9) 



Sa(X) 



= 0. 



(D.4) 



An important class of microscopic actions S [<p] are 
the ones possessing a global symmetry of the type de- 
scribed above, i.e. S*- a ' = S for constant a(X) = a. For 
a generic Lagrangian, we then find 



S W M= f L{ip (a) ,d^ a \...)= f £ (a) (a,d M c 
Jx Jx 



•) 



where 7° is the temporal component of the current. 

The global continuous symmetry of the microscopic 
action yields a conserved quantity Q. For ultracold 
atoms with global U(l) symmetry, for instance, Q - N 
is given by the particle number. Since the Noether 
charge is conserved globally, its local variations cannot 
be arbitrary, but rather must be such that the net change 
is locally zero. This fact is expressed by the continuity 
equation d^J^ = 0. 

We now turn our attention to a full quantum theory 
which is described by a microscopic action S[<p] and 
the partition function 

Z[j] = j Dipe- 51 *™* = J Yfyc-saW+i-** . (D.10) 

In the second line we performed a change of the integra- 
tion variable and assumed the functional measure to be 
invariant with respect to this mapping. This is the case 
for unitary transformations, which we want to consider 
here. However, there are cases (so-called anomalies), 
where the assumption D^ (cf) = Dtp is not valid. Since 



£ (a \a) = I £ = S[tp]. (D.5) the left hand side of Eq. ( |D~T0| does not depend on a 

X *JX \x/a finH frrvm thp> riaht hanH ciHp> tn firct nrHpr in rv thf 



Since this equation is valid to all orders in a, we find for 
an action with global symmetry the relation 



we find from the right hand side to first order in a the 
relation 



ia) 



da 



= 0. 



(D.6) 



/ 6S (a) 
' 6a(X) 



y- 

or=0 



6a(X) 



(D.ll) 



a=0 



a=0 



Whereas Eq. (D.4i is valid for any action evaluated 
at the solution to the equations of motion, Eq. ( |D.6| l is 
valid for all <p but only actions with global continuous 
symmetry. We can now combine both findings to obtain 



Switching to the effective action via a Legendre trans- 
formation in the variables (j,<p), we have j[<f>] = 
(dr/5<p)[4>] and deduce 



g S (a) 



6a(X) 



) =-• 

Jj 8<b 



6a(X) 



6F (a) 



6a(X) 



(D.12) 







6S 



(a) 



da 



a-=0,EoM 

d£ (a> 



da 



d£ (a) 
' did^a) 



did^a) 



u-=0,EoM 

(D.7) 



Analogous to the classical discussion, we introduced 

From the quantum action principle, i.e. 
(6T / S<P)[(/>q] = 0, we conclude that 



The leading term, which is relevant for most cases of 
interest, has a divergence structure. In fact, it constitutes 



6T (a) 



5a(X) 



= 0. 



(D.13) 



a-=0,EoM 
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This equation generalizes the classical relation (D.4i 



Together with the result in Eq. (D.12 1, we find for mi- 
croscopic actions with global continuous symmetry 







• 6a(X) 



=0,EoM 



d ( d£ia) ) 

' , \d(d„ar)/«=o,EoM 



(D.14) 



We again applied Eq. (D.6l. The ellipsis vanishes for a 
Lagrangian which only depends on the first derivatives 



of the field. Eq. (D.14i is a continuity equation for the 
full quantum Noether current 

with conserved charge Q - J d d xJ {) . 

Appendix E. Few-body physics in vacuum 

In this appendix, we show how the physics of a few par- 
ticles can be treated with the FRG. Special features of 
these vacuum calculations are diagrammatic simplifica- 
tions and a closed hierarchy of flow equations of the 
Af-body sector. For an introduction to this field see, for 
instance, Ref. [37|. Therein, emphasis is on the Efi- 
mov effect, but general features of vacuum physics with 
the FRG and its relation to the many-body problem are 
covered as well. 

The investigation of few-body problems with the 
FRG is motivated by several points. Firstly, the calcula- 
tion of vacuum observables from first principles allows 
for comparison with experiments or exact results from 
quantum mechanics and thus benchmarking the tech- 
nique. Moreover, the FRG approach reveals a different 
point of view on well-known results like the Efimov ef- 
fect. It also allows for the computation of nonuniversal 
features away from resonance. Finally, we already en- 
countered the importance of the vacuum problem for the 
many-body BCS-BEC crossover, where it provides the 
UV renormalization of the microscopic couplings and 
influences the physics on the BEC side. 

We obtain the effective action in vacuum by choosing 
the parameters such that n — T — with the tempera- 
ture always being above criticality. This implies con- 
densation to be absent. In an FRG setting, we then 
always remain in the symmetric phase with vanishing 
density. For the BCS-BEC crossover, we had found 
fi# = -mi < and m^>0 for the fermions and bosons, 
respectively. Since the ground state has to be stable, 
the propagators generically acquire a non-negative gap, 
r (2, (2 = 0) > 0. For the regularized propagator, we 



even have the stronger statement that their gaps are 
strictly positive. 



Simplification of the flow equations in vacuum 

All diagrams whose inner lines point into the same di- 
rection ( thereby forming a closed tour), do not con- 
tribute to the flow in vacuum^ This has important con- 
sequences for the FRG, because from Eq. ( |171| > we 
find the right hand side of the flow equation to be con- 
structed from one-loop diagrams. For this reason, sev- 
eral quantities do not get renormalized in vacuum (i.e. 
they do not receive contributions from quantum fluctu- 
ations), because there are simply no diagrams. To intu- 
itively understand the above claim, we consider such a 
cyclic diagram. It necessarily contains a line associated 
with a hole (or anti-particle). However, such excitations 
are not contained in the nonrelativistic quantum vacuum 
which conserves particle number. Therefore, these pro- 
cesses cannot occur. 

An example for a process which is described by 
such cyclic diagrams is provided by the particle-hole 
fluctuations around the Fermi surface which made up 
the Gorkov effect on the BCS side of the BCS-BEC 
crossover. The corresponding screening effect is not 
found in vacuum, because there is no Fermi surface. 

Cyclic diagrams vanish in vacuum, because all poles 
of the propagators lie in a definite half -plane of the com- 
plex w-plane. Therefore, the loop integral yields zero by 
virtue of the residue theorem. Indeed, consider a kinetic 
term § <p* (Q)Pb,f(Q)<p(Q) in the effective action for ei- 
ther bosons or fermions. In vacuum, there is no con- 
densation and the corresponding mean fields vanish. In 
the (<p, <^)*)-basis, we then find for the inverse propagator 
T (2) (2', = 5(Q + Q')G l (Q) the simple expression 



G b !f(0 = 



Pb,f(-0 
PbAQ) 



(El) 



Connected lines pointing into the same direction are 
represented by products of P B '(0 and P F '(0 with the 
same sign of the frequency and spatial momentum vari- 
able Q. Without loss of generality, we assume a> to be 
positive. As mentioned above, for a nonvanishing cutoff 
the masslike term in the propagator P(Q) will be strictly 
positive and the poles indeed lie in one distinct half- 
plane of the complex w-plane. We close the integration 



6 This holds under the mild assumption that the propagators are 
properly described in terms of a single pole. In other words, no addi- 
tional non-analyticities are generated during the RG flow. 
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contour in the opposite half-plane and the loop integral 
vanishes. Since the cutoff derivative d, only increases 
the multiplicity of the poles but not their location, the 
argument remains valid for the full flow equation. 

We give three examples of this diagrammatic sim- 
plification in vacuum, which are relevant for the BCS- 
BEC crossover. First, the fermion propagator does not 
get renormalized, because the mixed diagram contain- 
ing a fermion line and a boson line is cyclic. Indeed, 
this is a result of the form of the Yukawa coupling 
~ h((p*if/\if/2 + h.c). Second, the single box diagram 
which contributes to the flow of the four-fermion cou- 
pling A$ vanishes in vacuum. Thus we have A^^ = 
on all scales if A^^ = 0. The latter may be obtained by 
a suitable Hubbard-Stratonovich transformation of the 
microscopic action. This, in turn, implies the Yukawa 
coupling to be protected from renormalization, too, be- 
cause the only diagram contributing to its flow is propor- 
tional to Ay and thus zero. In particular, /u-(2i , Qi) = 
remains momentum independent. 

As a further consequence, we have a strict hierarchy 
of flow equations in vacuum. For instance, the two-body 
problem can be solved without knowledge of the three- 
(or more)-body problem. In general, the A^-body prob- 
lem only requires input only from processes involving 
M < N particles. This is easily understood in quantum 
mechanical terms, where the solution of the A^-body sys- 
tem is given by the Af -particle wave function. 

We are thus led to a very different truncation scheme 
for vacuum problems than we applied in the many-body 
sector. Indeed, it will turn out to be useful to expand the 
most general Tk[4>\ into monomials of the field, i.e. 



FuM +r*, 3 [0] +r kA [<p] + . 



(E.2) 



where Fyy is of A^-th order in the field <p. This a partic- 
ular case of a vertex expansion. 

Due to the translational invariance, we have for the 
momentum representation of a particular vertex 

r W (<2l, . . ., Q„) = JniQu- ■ ■ , Gn-lWGl + • • ■ + Qnl 

(E.3) 

The n-th momentum is determined by the remaining 
n - 1 ones due to momentum conservation. For ex- 
ample, we have T {1 \Q , Q) = 6(Q + Q')G \Q), which 
we already encountered several times in the earlier sec- 
tions. To get quantitative precision, one has to keep 
the full momentum dependence of the vertices y n . Be- 
cause of the simplifications which arise in vacuum, it 
might even in this case be possible to solve the system 
of equations. However, for a qualitative understanding, 
also the assumption of momentum independence can be 
employed. 



We dropped the zeroth order term in Eq. (E.2i, be- 
cause it represents the ground state energy of the vac- 
uum. The linear term vanishes for U(l)-symmetric sys- 
tem, which are the only relevant ones for our discussion. 
The higher order terms are determined such that they 
respect U(l)-symmetry and their total number of fields 
equals N. 



Two-body sectors and dimer binding energy for 
two-component fe unions 

From our above considerations on the non- 
renormalization of the fermion propagator, the 
four-fermion coupling and the Yukawa coupling, we 
see that the most general truncation to third order for a 
system of two-component fermions is given by 

T K2 [<p,ilf\ = J (^(<2)(iw - f-fx^MQ) 

- <P(Q^\Qx)^\Q2))s{Qi -Q.2- S 3 ). 

(E.4) 

The UV condition on the boson propagator reads 



^Va(G) = m \- Eq. iEAi contains everything that con- 



tributes to the two-body problem. The flow equation 

-h 2 



Jo 



q (P*{Q) + R^(Q))(P^(P -Q) + R*,k(P - Q)) 

(E.5) 

has been solved exactly for the cutoff = k 2 , and 
for the Litim cutoff discussed in the context of the BCS- 
BEC crossover. After UV renormalization, the physical 
dimer propagator is found to be 



=o(0 



lf_ 
8;r 



1 



id) q 1 
— + — 
2 4 



(E.6) 



The binding energy of the dimer is determined by the 
poles of the analytically continued propagator. We use 
our expression for the inverse propagator and compute 



^M=o(G = = £b/2) = 0. 



This yields 



e B = -- 



(E.7) 



(E.8) 



as expected. 
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